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ANNUAL  REPORT 


PCRP  060814 — “Transrectal  Near-Infrared  Optical  Tomography  for  Prostate  Imaging” 


INTRODUCTION 

The  objective  of  this  research  is  to  explore  the  technology  of  trans-rectal  near-infrared  (NIR) 
optical  tomography  for  accurate,  selective  prostate  biopsy.  Prostate  cancer  is  the  most  common 
non-dermatologic  cancer  in  American  men.  Prostate  cancer  suspicion  is  typically  based  on  an 
elevated  serum  prostate-specific  antigen  (PSA)  level  or  a  suspicious  nodule  found  during  a 
digital  rectal  exam  (DRE).  When  the  PSA  level  is  elevated  or  the  DRE  shows  abnormal,  there  is 
a  25  %  chance  that  cancer  is  present.  The  existence  of  prostate  cancer  can  only  be  confirmed 
by  a  needle  biopsy  that  is  guided  by  trans-rectal  ultrasound  (TRUS).  Since  there  are  no 
pathognomonic  findings  for  prostate  cancer  on  ultrasound  imaging,  biopsies  are  taken  following 
a  systematic  pattern  throughout  the  prostate  with  preference  given  to  the  peripheral  zone 
wherein  most  cancer  are  found.  The  accuracy  of  biopsy  is  questionable  and  many  men  undergo 
multiple  biopsies  due  to  the  lack  of  a  more  specific/sensitive  imaging  modality.  Pathologic 
studies  have  demonstrated  positive  correlation  between  increased  micro-vessel  density  and  the 
onset  of  the  disease.  Near-infrared  (NIR)  optical  tomography  is  known  of  sensitive  to  blood- 
based  contrast,  therefore  trans-rectally  implemented  NIR  optical  tomography  may  provide  a  new 
way  of  assessing  the  prostate  cancer.  One  of  the  outcomes  of  trans-rectal  NIR  tomography  of 
the  prostate  would  be  a  more  accurate  imaging  guidance  for  targeted  prostate  biopsy. 

BODY 


1.  Proposed  Specific  Aims: 

(1 )  To  demonstrate  that  endoscopic  NIR  tomography  at  a  probe  size  of  25mm  in  diameter  can 
be  achieved  by  use  of  spread-spectral-encoding  from  a  broad-band  light  source. 

(2)  To  demonstrate  that  trans-rectal  NIR  tomography  can  image  the  prostate  at  the  proximity  of 
the  rectum  with  significant  tumor-tissue  contrast. 

(3)  To  demonstrate  that  multi-spectral  trans-rectal  NIR  tomography  can  be  implemented  with  the 
single  trans-rectal  imaging  probe. 

(4)  To  demonstrate  that  trans-rectal  multi-spectral  NIR  tomography  can  quantify  the  hemoglobin 
concentration  and  oxygenation  saturation  in  phantom,  and  further  in  prostate  tumor  model  if  the 
time  of  research  allows. 

2.  Progresses  toward  trans-rectal  ultrasound-coupled  near-infrared  optical  tomography 

In  the  1st  year  of  this  research  (Mar,  2007 — Feb.  2008),  we  have  successfully  developed  an 
axial-imaging  20mm  diameter  trans-rectal  NIR  tomography  probe  which  is  equivalent  in  size  to  a 
standard  trans-rectal  ultrasound  (TRUS)  probe.  Calibrations  and  simulations  were  conducted  to 
evaluate  the  imaging  capabilities  of  this  axial-imaging  probe,  and  it  was  found  that  the  imaging- 
depth  of  this  axial  probe  may  be  limited  to  15mm  [Musgrove  2007].  In  vivo  axial  trans-rectal 
imaging  using  this  probe  has  been  conducted  on  healthy  canine  prostate  as  shown  in  Fig.  1 .  For 
this  study  a  planar  array  trans-rectal  US  probe  was  used  to  evaluate  the  location  and  ultrasonic 
structure  of  the  prostate  prior  to  imaging  with  this  stand-alone  axial-imaging  NIR  probe.  The 
depth  of  imaging  the  prostate  tissue  seemed  to  be  at  the  1cm  range.  The  tissue  being  imaged 
was  therefore  primary  the  peripheral  tissue.  The  limited  imaging  depth  may  be  caused  partially 
by  the  axial-imaging  geometry,  and  partially  by  the  strong  absorption  of  prostate  tissue. 
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We  have  since  investigated  the  approach  of  sagittal-imaging  trans-rectal  NIR  optical 
tomography  which  has  the  potential  of  interrogating  deeper  tissues,  and  developed  TRUS- 
coupled  NIR  tomography  for  concurrent  multimodality  trans-rectal  NIR/US  imaging. 


Figure  1.  Axial  trans-rectal  NIR  optical  tomography  of  normal  canine  prostate  in  vivo.  The 
absorption  value  of  the  prostate  tissues  agrees  with  that  in  Fig.  4. 


2.1  Principle  of  trans-rectal  NIR/US 

The  schematic  configuration  of  our  trans-rectal  NIR/US  prostate  imaging  is  illustrated  in  Fig.  2. 

In  order  to  interrogate  deep  prostatic  tissues,  it  is  necessary  to  have  longitudinally  distributed 
NIR  arrays  because  the  NIR  imaging  depth  is  approximately  !4  of  the  array  dimension.  By 
arranging  the  longitudinal  NIR  arrays  symmetrically  on  the  lateral  sides  of  the  sagittal  US 
transducer,  the  correlation  between  the  NIR  imaging  and  sagittal  US  can  be  made.  The  NIR 
arrays  could  have  a  number  of  optode  configurations.  However,  the  design  of  having  sources  or 
detectors  exclusively  on  one  of  the  two  NIR  arrays  is  adopted  after  evaluating  the  performance, 
the  endo-rectal  applicability,  and  the  prototyping  challenge  [Xu  et  al,  2008], 


NIR/US  imaging - Sagittal  view 


NIR  array 


NIR/US  imaging - Axial  view 


Prostate 


Tumor 


US  transducer 


Figure  2.  Illustration  of  trans-rectal  NIR/US  of  the  prostate.  Trans-retal  US  is  placed  in  the 
middle  of  the  trans-rectal  NIR  applicator  (optodes  distributed  longitudinally)  to  perform 
combined  NIR/US  imaging. 

The  completed  trans-rectal  NIR/US  probe  and  imager  are  shown  in  Fig.  3.  The  probe 
was  developed  based  on  the  ALOKA  Veterinary  US  unit.  The  NIR  applicator  was  integrated 
over  the  intact  7.5MHz  sagittal-imaging  transducer  of  the  bi-plane  ALOKA  trans-rectal  US 
probe.  The  ALOKA  US  probe  also  had  a  5MHz  axial-imaging  transducer.  The  NIR  probe 
consisted  of  one  source  and  one  detector  array  separated  laterally  20mm  and  placed  symmetric 
to  the  sagittal  US  transducer.  The  NIR  source  and  detector  array,  each  having  7  channels,  were 
60mm  in  longitudinal  dimension,  which  was  identical  to  the  longitudinal  imaging  view  of  the 
sagittal  US  transducer.  Each  NIR  optode  channel  had  a  micro  prism-lens  pair  for  coupling  the 
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light  to  and  from  the  probe  surface.  The  NIR  light  from  a  super-luminescent  diode  of  lOOmW  at 
840nm  was  focused  sequentially  onto  7  source  fibers  of  the  NIR  applicator  by  a  home-made 
translating  fiber  multiplexer.  The  NIR  remissions  collected  by  the  7  detection  fibers  were 
acquired  by  a  16-bit  intensified  CCD  camera  through  a  spectrometer.  Acquisition  of  signals  for 
all  7  source  channels  usually  took  5  seconds  after  the  prostate  was  localized  on  the  sagittal 
trans-rectal  US.  An  Aloka-900  (in  our  lab)  or  Aloka-3500  (belonging  to  OSU  Vet-Med)  US 
scanner  acquired  the  bi-plane  US  images,  of  which  the  sagittal  US  had  accurate  position 
correlation  with  the  trans-rectal  NIR.  The  absorption  and  reduced  scattering  images  of  the 
prostate  were  reconstructed  from  steady-state  measurements  [Farrell  1992]  using  non-linear 
optimization  method  [Xu,  2005,  Iftimia  2000,  Jiang  et  al,  2008], 

Using  this  probe/imager,  trans-rectal  NIR/US  has  been  performed  on  normal  and 
diseased  canine  prostate  in  vivo.  Dogs  have  been  previously  used  for  prostate  studies  because 
of  the  similarity  between  canine  and  human  prostate  glands  [Hazle  2002;  Forsberg  2002;  Rivera 
2005],  We  have  successfully  developed  a  canine  prostate  tumor  model  in  OSU  Center  for 
Veterinary  Health  Sciences  with  this  New  Investigator.  This  study  was  approved  by  the 
university's  Institutional  Animal  Care  and  Use  Committee.  The  protocol  was  also  approved  and 
underwent  an  on-site  inspection  by  the  US  Army  Medical  Research  and  Material  Command.  In 
the  following  sections  we  will  discuss  the  results  obtained  from  normal  canine  prostate  and 
canine  prostate  cancer  models.  It  should  be  noted  that  all  the  in  vivo  NIR  images  presented 
were  reconstructed  purely  from  NIR  signal,  without  any  US  priori  information. 


NIR/US  probe  based  on  Aloka  US  transducer 


Spectrometer 
&  CCD 


NIR/US  imaging  system  configuration 


Figure  3  (Left).  Photograph  of  the  Aloka  trans-rectal  US  transducer  and  the  completed 
NIR/US  probe  used  for  canine  imaging.  (Right)  Schematic  diagram  of  the  trans-rectal 
NIR/US  imaging  system  that  consists  of  a  custom-built  NIR  imager  and  a  commercial 
ALOKA  SSD-3500  US  scanner. 

2.2.  In  vivo  trans-rectal  NIR/US  of  normal  canine  prostate 

Figure  4  displays  the  trans-rectal  NIR/US  imaging  results  performed  on  normal  canine  prostate 
in  vivo.  The  sagittal  US  was  conducted  at  the  middle  of  the  right  lobe,  the  middle-line  of  the 
prostate  gland,  and  the  middle  of  the  left  lobe.  The  trans-rectal  NIR  images  were  in  3-dimension 
but  only  the  sagittal  plane  aligned  with  the  US  were  displayed  here.  The  NIR  images  were  for 
absorption  coefficient  and  reduced  scattering  coefficient,  respectively.  The  prostate  was  shown 
as  hyper-absorbing  (as  prostate  is  more  vascular),  and  slightly  hyper-scattering,  with  respect  to 
the  peripheral  tissues.  The  absolute  values  as  well  as  the  large  ranges  of  the  absorption  and 
reduced  scattering  coefficients  of  the  prostate  agree  with  the  values  reported  in  literatures 
[Svensson  et  al,  2007  reviewed  optical  properties  of  human  prostate]. 
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Figure  4.  Trans-rectal  NIR/US  of  normal  canine  prostate  in  vivo:  The  US  images  were 
taken  at  the  right  lobe,  middle-line,  and  the  left  lobe.  (Left),  US/NIR  absorption;  (right) 

US/NIR  reduced  scattering.  The  dimensions  of  the  images  are  50mmx30mm  (cranial- 
caudal  x  dorsal-ventral). 

2.3.  In  vivo  trans-rectal  NIR/US  of  transmissible  venereal  tumors  in  canine  pelvic  canal 

We  have  developed  a  canine  prostate  cancer  model  using  canine  transmissible  venereal  tumor 
(TVT).  The  unique  canine  TVT  is  a  round  cell  tumor  of  dogs  that  mainly  affects  the  external 
genitalia  and  can  be  transmitted  from  animal  to  animal  during  copulation,  regardless  of 
histocompatibility.  The  neoplastic  cells  are  generally  thought  to  be  of  histiocytic  origin  [Mozos  et 
al,  1996],  TVT  cells  can  be  propagated  in  immunocompromised  (SCID)  mice  and  transferred  to 
different  tissues  (such  as  the  prostate  gland)  of  the  dog  to  result  in  a  neoplastic  mass  effect 
useful  for  imaging  studies  [Rivera  2005].  For  these  studies,  an  existing  TVT  cell  line  was 
obtained  as  cryopreserved  tissue  from  MD  Anderson  Cancer  Center  (Houston,  TX).  The  tumor 
tissue  was  quick-thawed,  homogenized  and  inoculated  into  the  subcutis  of  a  NOD/SCID  mouse. 
After  approximately  12  weeks,  the  tumor  reached  an  appropriate  volume  (4. 0-4. 5cm3)  for 
recovery  and  were  processed  and  inoculated  into  a  second  NOD/SCID  mouse.  Following 
growth  in  the  second  mouse  (8  weeks),  the  neoplastic  cells  were  recovered  and  homogenized 
for  inoculation  into  the  canine  prostate  gland. 

A  12  kg  sexually  intact  adult  Beagle  dog  estimated  to  be  approximately  4  years  of  age 
was  housed  and  administered  using  general  anesthesia  for  TVT  cell  injection.  A  6-inch  16 
gauge  hypodermic  needle  was  inserted  transperineally  into  the  right  lobe  of  the  prostate  where 
~3  cc  of  the  prepared  TVT  cells  were  injected.  During  retraction  of  the  injection  needle,  it  was 
considered  unavoidable  that  TVT  cells  would  leak  from  the  prostate  injection  site  and  be 
"seeded"  along  the  needle  insertion  tract.  The  dog  was  examined  weekly  by  US  using  the  trans- 
rectal  NIR/US  probe  covered  by  condom  as  did  with  standard  trans-rectal  US.  During  the  first 
two  weeks  following  TVT  cell  injection  there  was  no  evidence  of  tumor  growth  on  US  and  rectal 
examination.  The  next  US  examination  occurred  5  weeks  post-injection  when  hypoechoic 
masses  were  observed  both  in  the  prostatic  parenchyma  and  peri-prostatically  around  the  right 
lobe  of  the  prostate.  The  dog  underwent  weekly  trans-rectal  US  monitoring  for  2  more  weeks, 
and  was  then  humanly  euthanized  8  weeks  after  the  initial  TVT  cell  injection  for  a  complete 
necropsy  and  histological  examination  of  prostate  and  peri-prostatic  tissues. 
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The  NIR  absorption  and  reduced  scattering  images  taken  at  5  weeks  post-injection  are 
displayed  in  Fig.  5(a)  and  (b),  respectively.  The  US  and  NIR  images  have  the  same  dimension 
of  60mmx30mm  (cranial-caudalxdorsal-ventral)  for  sagittal,  40mmx30mm  (right  lateral-left 
lateralxdorsal-ventral)  for  axial,  and  60mmx40mm  (cranial-caudalxright  lateral-left  lateral)  for 
coronal  views,  respectively.  The  sagittal  trans-rectal  US  images  were  taken  at  the  middle 
portion  of  the  right  lobe.  The  hypo-echoic  region  LI  indicates  an  intra-prostatic  mass.  The  large 
hypo-echoic  region  L2  indicates  a  mass  ventral  and  caudal  to  the  prostate  that  may  have 
connection  with  LI .  The  NT  on  the  right  lobe  US  image  denotes  the  needle  trajectory  for 
introducing  the  TVT  cells.  Longitudinal  hypo-echoic  regions  including  L3  are  seen  along  the  NT. 
The  hyper-absorptive  regions  on  sagittal  NIR  images  correspond  to  LI,  L2,  and  L3  on  sagittal 
US  image.  The  lOmm-medial  NIR  image  displays  the  absorptive  masses  at  reduced  NIR 
contrast,  and  the  lOmm-lateral  NIR  reveals  connected  strong  absorptive  masses.  The  bladder  is 
seen  as  slightly  hypo-echoic  at  other  medial  sagittal  US  images  (not  shown  here);  however,  no 
hyper-absorptive  mass  corresponding  to  the  bladder  is  presented  at  the  left-most  region  of  the 
sagittal  NIR  images.  The  three  axial  TRUS  images  were  taken  at  the  cranial  edge  of  the 
prostate  crossing  LI ,  the  caudal  edge  of  the  prostate  crossing  L2,  and  the  peri-rectal  region 
crossing  L3.  The  axial  US  images  show  a  small  hypo-echoic  intra-prostatic  mass  at  cranial  side 
of  the  prostate,  the  distortion  of  the  right  lobe  boundary  and  extending  of  the  L2  over  the 
prostate  mid-line  that  correlates  with  the  findings  on  mid-line  sagittal  US  (shown  in  Fig.  6),  and 
large  peri-rectal  hypo-echoic  mass  cranial  to  the  perineum.  The  locations  of  axial  NIR  images 
are  comparable  to  those  of  the  axial  trans-rectal  US.  The  axial  NIR  image  A2  should  contain  LI, 
and  the  axial  NIR  image  A5  should  contain  L3.  The  coronal  NIR  images  correspond  to  5mm, 
10mm,  and  15mm  ventral  to  the  rectum.  The  no-show  of  the  lesions  on  the  15mm  coronal 
image  is  due  to  the  decreased  NIR  sensitivity  of  side-firing  and  the  reconstruction  without  any 
spatial  prior.  The  hyper-absorptive  mass  indicative  of  LI  is  seen  medial  to  the  hyper-absorptive 
masses  indicative  of  L2  and  L3.  The  appearance  of  the  NIR  hyper-absorptive  masses  in  3-views 
agrees  with  that  of  the  US  hypo-echoic  regions  in  2-views. 

The  hyper-absorptive  masses  in  (a)  being  suspicious  of  TVTs  are  shown  in  (b)  with 
different  contrast  patterns  for  the  reduced  scattering  coefficient.  In  the  NIR  scattering  images, 

L3  has  much  higher  contrast  than  do  the  other  masses  corresponding  to  LI  and  L2.  The 
contrast  elsewhere  is  also  more  homogenous  compared  with  that  in  the  absorption  image. 

In  Fig.  6  we  present  the  retrospective  comparisons  of  the  trans-rectal  US  and  the  NIR 
images  obtained  at  the  middle  of  the  right  lobe,  the  middle-line  of  the  prostate,  and  the  middle  of 
the  left  lobe,  before  the  TVT  injection,  2  weeks  after  the  TVT  injection  when  the  US  and  rectal 
examination  showed  no  evidence  of  tumor  growth,  and  5  weeks  after  the  TVT  injection  when 
the  tumor  growth  was  evident  on  both  US  and  rectal  examination.  Both  the  absorption  and 
reduced  scattering  images,  correlating  to  the  same  set  of  US  images,  are  shown.  At  the  2-week 
NIR  absorption  images,  hyper-absorptive  regions  were  found  dorsal  to  the  prostate  in  the  right 
lobe  only,  and  dorsal  to  the  pelvic  bone  in  the  right  lobe  with  potential  extension  to  the  middle 
line.  The  locations  of  these  hyper-absorptive  regions  correlated  very  well  with  those  of  the  hypo- 
echoic  and  hyper-absorptive  regions  found  in  week-5.  The  NT  in  week-5  US  image  indicated  the 
needle  track  for  injecting  the  tumor  cells.  The  dimensions  and  contrasts  of  the  hyper-absorptive 
regions  implied  a  tumor  growing  in  right  lobe  and  extending  toward  the  middle-line.  The  growth 
of  the  tumor  was  indicated  earlier  by  the  absorption  images  of  trans-rectal  NIR  than  by  the 
trans-rectal  US.  The  growth  of  the  tumor  on  reduced  scattering  images,  although  not  as  evident 
as  in  the  absorption  images,  did  correlate  with  the  findings  in  trans-rectal  US.  It  is  evident  that 
combining  the  information  of  trans-rectal  US  and  trans-rectal  NIR  renders  earlier  and  much 
more  accurate  findings  of  the  tumor  growth. 
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(b) 

Figure  5  TVTs  developed  in  canine  pelvic  canal  assessed  by  trans-rectal  NIR  and  US  in 
vivo  at  5-weeks  post  TVT-injection.  The  NIR  images  are  in  sagittal,  axial,  and  coronal 
views.  The  US  images  are  in  sagittal  and  axial  views.  The  dimensions  of  the  images  are: 
60mmx30mm  (cranial-caudalxdorsal-ventral)  for  sagittal,  40mmx30mm  (right  lateral-left 
lateralxdorsal-ventral)  for  axial,  and  60mmx40mm  (cranial-caudalxright  lateral-left 
lateral)  for  coronal  views,  respectively. 
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Figure  6  Comparison  between  the  image-monitoring  of  TVT  development  in  canine  pelvic 
canal  assessed  by  trans-rectal  NIR  and  US  in  vivo,  over  a  period  of  5  weeks.  The  images 
were  taken  at  the  middle  of  the  right  lobe,  the  middle-line  of  the  prostate,  and  the  middle 
of  the  left  lobe,  before  the  TVT  injection,  2  weeks  after  the  TVT  injection  when  the  US 
and  rectal  examination  showed  no  evidence  of  tumor  growth,  and  5  weeks  after  the  TVT 
injection  when  the  tumor  growth  was  evident  on  both  US  and  rectal  examination. 
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The  gross  and  histological  findings  (taken  at  8-week  after  the  TVT  injection)  given  in  Fig. 
7  confirmed  intra-  and  peri-prostatic  neoplastic  infiltrates  with  masses  also  located  along  the 
urethra  and  peri-rectal  tissue;  the  latter  related  to  dissemination  along  the  needle  track  during 
TVT  inoculation.  All  masses  consist  of  diffuse  sheets  of  a  monomorphic  population  of  neoplastic 
round  cells  dissecting  through  pre-existing  fibrovascular  stroma.  The  neoplastic  cells  have  large 
hyperchromatic  nuclei,  single  conspicuous  nucleoli  and  moderate  amounts  of  featureless 
cytoplasm.  The  cytological  features  are  consistent  with  canine  transmissible  venereal  tumor 
(TVT). 


Figure  7.  (Top)  Postate  gland  taken  after  removal  of  the  urinary  bladder  and  intra-pelvic 
urethra  immediately  following  euthanasia.  The  prostate  gland  is  markedly  expanded  and 
misshapen  secondary  to  intra-prostatic  (LI)  and  peri-prostatic  (L2)  neoplastic  masses. 
Neoplastic  masses  were  also  distributed  along  the  distal  intra-pelvic  urethra  (L3)  and  in 
the  peri-rectal  subcutaneous  tissues  (not  shown).  Histologically,  expansile  masses  of 
transmissible  venereal  tumor  cells  (TVT,  bottom  left  panel)  displace  and  compress 
prostatic  tissue  (P).  Prostatic  capsule  (arrowheads),  Bar  =  230pm.  The  TVT  infiltrate 
(bottom  right  panel)  consists  of  sheets  of  neoplastic  round  cells  dissecting  through  pre¬ 
existing  fibrovascular  stroma  (arrowheads),  Bar  =  300pm. 

Among  the  US  hypo-echoic  masses  obtained  in  week-5,  the  LI  was  intra-prostatic,  the 
L2  was  largely  ventral  to  the  prostate,  and  the  L3  was  caudal  to  the  prostate.  These  3  masses 
were  shown  at  similar  depths  on  NIR  images,  which  was  partially  due  to  the  depth-dependent 
NIR  sensitivity  [Xu  et  al,  2008]  and  partially  due  to  the  large  volume  of  diffuse  photon 
interrogation  with  the  tumor.  Despite  the  depth  inaccuracy,  the  TVT  nodules  had  been  detected 
by  use  of  trans-recta  I  NIR  alone.  If  the  image  reconstruction  of  trans-rectal  NIR  could  be  guided 
by  the  US  structural  information  of  the  prostate  (and  the  tumor  if  the  tumor  is  visible  to  US),  the 
characterization  of  the  TVT  nodules  on  NIR  would  become  more  accurate. 

It  is  known  that  most  prostate  cancers  are  presented  as  multi-focal  [Wise  et  al,  2002].  It 
must  be  noted  that  the  intraprostatic  TVT  tumors  were  initiated  in  a  non-immunosuppressed 
canine  model  on  which  the  TVT  nodules  developed  at  multiple  sites  intra-prostatically  and  peri- 

8 


prostatically.  Although  not  all  of  the  TVT  tumors  were  confined  to  the  prostate,  the  successful 
imaging  of  multiple  TVT  nodules  by  trans-rectal  NIR  tomography  implies  the  utility  of  detecting 
multiple  intra-prostatic  tumors.  On  trans-rectal  US  all  of  the  TVT  tumor  nodules  developed  in 
this  subject  were  shown  strongly  hypo-echoic  that  correlated  with  the  strongly  hyper-absorptive 
and  moderately  hyper-scattering  findings  on  NIR  imaging.  The  bladder,  on  the  other  hand, 
showed  as  hypo-echoic  on  US;  however  it  was  shown  as  iso-absorptive  or  iso-scattering  on 
NIR.  This  suggests  that  trans-rectal  NIR  tomography  has  higher  specificity  in  imaging  the 
regions  hypo-echoic  to  TRUS. 

The  feature  of  TVT  as  strongly  hyper-absorptive  on  NIR  tomography  is  likely  due  to  the 
hyper-chromatic  nuclei  unique  to  TVT.  The  relatively  higher  scattering  of  the  TVT  may  be  due  to 
the  hyper-density  and  larger  nuclei  of  the  neoplastic  cells  [Wang  et  al,  2006].  As  the  neoplastic 
cells  are  arranged  into  microlobules  by  the  pre-existing  fibrovascular  stroma,  the  TVT  may  also 
present  certain  polarization-sensitivity.  Overall,  the  characteristic  features  of  TVT  on  NIR 
tomography  are  comparable  to  those  of  breast  cancer  angiogenesis  on  NIR.  The  micro¬ 
architecture  of  spontaneous  prostate  cancer  may  be  different  from  that  of  TVT.  So  the  intrinsic 
NIR  contrast  of  human  prostate  cancer  may  not  be  as  pronounced  as  those  of  TVT,  but  it  can 
be  expected  that  the  human  prostate  cancer  will  have  notable  contrast  on  trans-rectal  NIR 
tomography  if  the  contrast  is  taken  for  hemoglobin  concentration,  and  oxygenation  saturation. 
These  contrasts,  which  are  difficult  for  trans-rectal  US  to  acquire,  could  help  to  characterize  the 
malignancy  of  the  prostate  lesions. 

2.4.  Detecting  single  prostate  cancer  using  the  prototype  NIR/US  probe — simulation 
studies 

We  have  also  conducted  simulation  studies  on  the  performance  of  prostate  cancer  detection 
using  the  above-mentioned  prototype  NIR/US  probe  and  a  prior-guided  NIR  reconstruction 
method  developed  specifically  for  prostate  NIR  tomography  [Xu  et  al,  2008],  The  purpose  was 
to  evaluate  the  detectability  of  the  cancer  when  it  shows  conspicuous  on  US,  therefore  it  was 
assumed  that  the  morphological  information  of  the  prostate  and  lesion  was  available  from  trans- 
rectal  US.  This  is  not  entirely  representative  of  the  actual  clinical  cases  as  many  cancers  are  not 
visible  on  trans-rectal  US,  but  it  will  provide  valuable  insights.  Trans-rectal  US  images  of  human 
prostate  available  in  open  sources  were  used  for  the  simulation  study.  The  TRUS  image  was 
first  converted  to  COMSOL  Multiphysics  compatible  format.  The  prostate  tumor  is  mimicked 
using  a  spherical  shape  to  allow  for  flexibility  of  adjusting  its  size.  The  absorption  and  reduced 
scattering  coefficients  of  the  rectum,  peri-prostate  tissue,  prostate  and  the  tumor  are  assigned 
with  values  as  suggested  by  literature  [Li,  2007]. 


Figure  8  Finite-element  mesh  of  the  prostate  for  simulation 

Figure  8  illustrates  one  example  of  the  completed  FEM-mesh  for  trans-rectal  NIR  optical  j 
tomography  derived  from  a  trans-rectal  US  image,  where  x,  y  and  z  denote  the  longitudinal, 
lateral  and  the  depth  coordinates,  respectively.  A  typical  mesh  used  for  this  work  contains 
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approximately  4000  nodes  and  20000  linear  tetrahedral  elements.  The  mesh  in  the  figure 
corresponds  to  a  volume  of  80x80x80  mm3  and  the  ‘walnut’  shaped  prostate  has  a  dimension  of 
50x50x30  mm3.  The  rectum  wall  is  4mm  thick  with  a  curvature  radius  of  80  mm.  The  choice  of 
the  curvature  radius  is  due  to  the  fact  that  the  NIR  array  added  to  a  TRUS  probe  may  have  a  flat 
surface  that  would  transform  the  rectum  lumen  to  an  elliptical  shape.  The  source  modulation 
frequency  is  set  at  100MHz,  and  1%  Gaussian  noise  is  added  to  all  forward  calculations  to  form 
the  measurement  data.  Comparisons  of  detection  sensitivity  are  made  between  the  prototype 
probe  of  7sd  (indicating  7  source  or  detector  on  one  lateral  side  of  US)  and  an  NIR  probe  of 
26ssdd  (indicating  26  source  or  detector — separated  to  2  lines--on  each  lateral  side  of  US),  both 
having  6cm  span  of  NIR  array. 

2.4. 1  Reconstruction  accuracy  versus  target  longitudinal  location 

A  target  of  10mm  in  diameter  is  placed  at  the  middle-sagittal  plane  of  y=40mm,  z=26mm,  and 
varied  in  longitudinal  coordinates  from  x=25  to  55mm  with  a  step  of  5mm  (Fig.  9(a)).  The  optical 
properties  reconstructed  by  the  two  geometries  are  plotted  versus  the  true  values  in  Fig.  9(b) 
and  (c).  The  optical  properties  recovered  by  26ssdd  and  7sd  are  close  to  each  other  at  most  of 
the  longitudinal  locations,  but  the  7sd  design  shows  a  larger  variation  in  the  recovered 
absorption  contrast  at  x=30mm  and  x=50mm  compared  to  other  positions.  This  may  be  related 
to  fewer  source-detector  pairs  that  contribute  to  the  target  detection  when  close  to  the  boundary 
or  the  existence  of  any  irregular  elements  in  the  mesh. 
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Figure  9  Comparison  of  two  geometries  for  a  target  varying  in  longitudinal  location  in  the 
middle-sagittal  plane:  (a)  illustration  of  the  target  location  change  [Unit:  mm'1];  (b) 
comparison  of  absorption  coefficient  reconstruction;  (c)  comparison  of  reduced  scattering 
coefficient  reconstruction 

2.4.2  Reconstruction  accuracy  versus  target  depth 

A  target  of  10mm  in  diameter  is  placed  at  the  middle-sagittal  plane  at  x=40mm,  y=40mm,  and 
the  depth  is  varied  from  z=25  to  40mm  at  a  step  of  2.5mm  (the  last  data  point  is  simulated  at 
z=39mm  as  at  40mm  the  target  is  out  of  the  prostate)  (Fig.  10(a)).  The  reconstructed  optical 
properties  are  plotted  in  Fig.  10(b)  and  (c).  The  26ssdd  configuration  outperforms  the  7sd  one 
again,  however  beyond  z=30mm,  both  designs  are  incapable  of  recovering  the  absorption 
coefficient  of  the  target  from  the  prostate  background.  This  depth  limitation  is  related  to  the 
maximum  span  of  the  NIR  array,  the  large  absorption  coefficient  (0.06mm'1)  of  the  prostate,  and 
the  size  of  the  target.  For  a  larger  target  with  a  diameter  of  14mm,  it  is  verified  that  the  target 
may  be  resolved  up  to  36mm  depth  from  the  NIR  array  in  comparison  to  30mm  for  a  target  of 
10mm  diameter.  A  potentially  smaller  absorption  coefficient  of  intact  prostate  may  also  increase 
the  depth  limit  of  target  detection  due  to  an  increase  in  the  depth  sensitivity. 
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Fig.  10  Comparison  of  two  geometries  for  a  target  varying  in  depth  in  the  middle-sagittal 
plane:  (a)  illustration  of  the  target  location  change  [Unit:  mm"1];  (b)  comparison  of 
absorption  coefficient  reconstruction;  (c)  comparison  of  reduced  scattering  coefficient 
reconstruction 

2.4.3  Reconstruction  accuracy  versus  target  size 

A  target  is  placed  at  middle-sagittal  plane  of  x=40mm,  y=40mm  and  z=26mm,  and  the  diameter 
of  the  target  is  varied  from  4mm  to  14mm  with  a  step  of  1  mm.  The  change  of  the  target 
diameter  is  illustrated  in  Fig.  11(a).  The  reconstructed  optical  properties  are  shown  in  Fig.  11(b) 
and  (c).  It  is  clear  that  the  larger  the  target,  the  better  the  accuracy  of  reconstruction.  The 
26ssdd  can  recover  the  absorption  contrast  of  the  target  when  the  diameter  is  greater  than  6mm 
and  the  7sd  can  recover  the  target  for  target  diameter  greater  than  8mm. 

These  comparisons  suggest  that  the  accuracy  of  recovering  the  absorption  properties  is 
strongly  dependent  on  the  number  of  NIR  optodes,  which  weakly  affects  that  of  scattering 
properties. 
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Fig.  1 1  Comparison  of  two  geometries  for  a  target  varying  in  size  in  the  middle-sagittal 
plane:  (a)  illustration  of  the  target  size  change  [Unit:  mm"1];  (b)  comparison  of  absorption 
coefficient  reconstruction;  (c)  comparison  of  reduced  scattering  coefficient  reconstruction 

2.5.  Detecting  dual-center  prostate  cancer  using  the  prototype  NIR/US  probe — simulation 
studies 

Capability  of  differentiating  two  targets  is  of  particular  relevance  to  prostate  cancer  imaging 
owing  to  the  existence  of  secondary  or  multifocal  tumors.  In  this  section,  we  investigate  only  the 
case  of  having  two  targets  in  the  same  sagittal  plane  and  conduct  the  simulations  with  the  7sd 
probe  design  only. 

2. 5. 1  Reconstruction  of  two  targets  located  at  the  same  depth  in  sagittal  plane 
Figure  12  shows  two  lOmm-diameter  regions  being  added  to  the  prostate,  at  coordinates  (25, 
40,  26)  and  (55,  40,  26),  respectively.  In  (a)  only  one  region  has  optical  contrast,  and  in  (b)  both 
regions  have  optical  contrasts.  In  both  cases  the  optical  contrast  can  be  reconstructed  with 
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good  accuracy.  The  /ua  of  the  target  with  true  optical  contrast  is  reconstructed  within  ±20%  of 

the  set  value  and  jus  can  be  reconstructed  within  ±23%  of  the  set  values.  The  target  with  no 

optical  contrast  is  reconstructed  with  some  artifacts,  nevertheless,  the  target  with  optical 
contrast  can  be  easily  differentiated  from  the  one  without. 


Fig.  12  Two  suspicious  regions  at  the  same  depth  [Units:  mm  and  mm'1] 

2.5.2  Reconstruction  of  two  targets  located  at  different  depth  in  sagittal  plane 
Two  targets  of  10mm  diameter  are  added  in  the  prostate  at  coordinates  of  (25,  40,  28)  and  (55, 
40,  24),  respectively.  Figure  13  shows  the  reconstructed  images  for  the  case  of  both  targets 
having  negative  contrast. 
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Fig.  13  Two  targets  at  different  depths:  negative  tumor  contrast  cases 

[Unit:  mm  and  mm"1] 

When  the  target  is  at  a  depth  of  24mm,  the  ^  and  can  be  reconstructed  within  ±20% 

and  ±25%  of  the  set  values,  respectively.  However,  a  target  at  28mm  depth  cannot  be 
reconstructed.  This  is  due  to  the  high  absorption  coefficient  of  the  prostate  set  at  0.06mm-1 . 
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When  the  prostate  absoption  coefficient  is  reduced  to  0.006mm-1  which  will  provide  positive 
optical  contrast  in  the  two  target  regions,  both  tartgets  can  be  recovered  as  shown  in  Fig.  14. 
The  jua  of  the  target  can  be  reconstructed  within  ±5%  of  the  set  value,  while  the  is  still 

reconstructed  within  ±23%  of  the  expected  values. 


Fig.  14  Two  targets  at  different  depths:  positive  tumor  contrast  cases  [Unit:  mm  and  mm'1] 


KEY  RESEARCH  ACCOMPLISHMENTS 

The  following  research  accomplishments  have  been  made  during  the  2nd  year  of  this  project: 

1.  The  20mm  diameter  axial-imaging  trans-rectal  NIR  tomography  probe  has  been  tested  on 
imaging  healthy  canine  prostate  in  vivo.  Prostate  tissues  in  the  peripheral  zone  were  imaged 
and  reconstructed  with  optical  contrasts  agree  with  literature  values. 

2.  A  TRUS-coupled  trans-rectal  NIR  tomography  probe  has  been  developed.  The  combined 
probe  performs  sagittal  US  and  NIR  imaging  concurrently.  The  NIR  imaging  at  the  middle- 
sagittal  plane  correlates  with  the  TRUS  image. 

3.  In  vivo  trans-rectal  NIR/US  imaging  of  normal  canine  prostate  has  been  demonstrated. 

4.  A  canine  prostate  tumor  model  has  been  developed  using  canine  TVT  cells.  We  now  have  a 
reproducible  large  animal  prostate  tumor  model. 

5.  In  vivo  trans-rectal  NIR/US  imaging  of  canine  TVT  development  has  been  demonstrated. 
There  is  evidence  that  when  TRUS  is  aided  by  trans-rectal  NIR,  the  tumor  development  can  be 
detected  earlier  and  more  accurately. 

6.  Quantitative  evaluation  of  the  imaging  performance  of  the  combined  trans-rectal  NIR/US 
imaging  has  been  conducted  by  simulation.  The  evaluation  includes  both  detecting  single  tumor 
and  multi-centric  tumors. 
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REPORTABLE  OUTCOMES 


The  progress  in  the  2nd  year  of  this  project  has  results  in  following  publications  or  manuscripts: 

Provisional  Patent  Applications 

1.  Daqing  Piao,  Zhen  Jiang,  Guan  Xu,  Charles  F.  Bunting,  “Method  and  Apparatus  for 
combined/concurrent  trans-channel  near-infrared  optical  tomography  and  ultrasound  for 
imaging  internal  organs”,  US  provisional  patent  application,  #61/144,556,  filed  on  Jan.  16, 
2009. 

Journal  Papers 

1.  Jiang  Z,  Piao  D,  Xu  G,  Ritchey  JW,  Holyoak  GR,  Bartels  KE,  Bunting  CF,  Slobodov  G, 
Krasinski  JS,  “Trans-rectal  ultrasound-coupled  near-infrared  optical  tomography  of  the 
prostate  Part  II:  Experimental  demonstration,”  Optics  Express,  16(22):  17505-17520  (2008). 
[PubMed  index:  18958031], 

2.  Xu  G,  Piao  D,  Musgrove  CH,  Bunting  CF,  Dehghani  H,  “Trans-rectal  ultrasound-coupled 
near-infrared  optical  tomography  of  the  prostate  Part  I:  Simulation,”  Optics  Express.  16(22): 

1 7484-1 7504  (2008).  [PubMed  index:  1 8958030], 

Proceeding  Papers  and  Conference  Presentations 

3.  Jiang  Z,  Ritchey  JW,  Holyoak  GR,  Bartels  KE,  Xu  G,  Bunting  CF,  Slobodov  G,  Krasinski  JS, 

Piao  D,  “In  vivo  trans-rectal  ultrasound  coupled  trans-rectal  near-infrared  optical  tomography 
of  canine  prostate  bearing  transmissible  venereal  tumor,”  Proceedings  of  SPIE,  Vol.  7174, 
Paper  #71741 U.. 

4.  Xu  G,  Bunting  CF,  Dehghani  H,  Piao  D,  “A  hierarchical  spatial  prior  approach  for  prostate 

image  reconstruction  in  trans-rectal  optical  tomography,”  Proceedings  of  SPIE,  Vol.  7171, 
Paper #71 71  OS. 

5.  Jiang  Z,  Piao  D,  Krasinski  JS,  “Development  of  a  continuous-wave  dual-band  trans-rectal 
optical  tomography  system  for  concurrent  sagittal  imaging  with  trans-rectal  ultrasound,” 
Proceedings  ofSPIE,  Vol.  71 71 ,  Paper  #71 71 0G. 

6.  Piao  D,  Holyoak  GR,  Bartels  KE,  Ritchey  JW,  Jiang  Z,  Xu  G,  Bunting  CF,  “In  vivo  trans-rectal 

optical  tomography  of  normal  canine  prostate — demonstration  of  optical  contrast  of  intact 
prostate  over  its  peripheral  tissue,”  Saratov  Fa//  Meetings  2008,  Saratov,  Russia,  Sep.  23- 
26,  2008  (Internet  session  invited  lecture). 

Students  Fully  or  Partially  Supported  During  The  2nd  Year: 


1.  JIANG,  Zhen 

PhD  student 

2.  XU,  Guan 

PhD  student 

3.  Mukherjee,  Sovanlal 

PhD  student 

4.  XIE,  Hao 

MS  degree  in  May,  2008 

Grants  Being  Awarded  That  Can  Be  Attributed  To  This  Research: 

1.  Big  XII  Faculty  Fellowship  Program  —  Faculty  Fellow  Award 

Piao,  D  (PI),  02/2009— -08/2009,  $2,500 

Title:  Faculty  Fellowship  in  University  of  Missouri,  Columbia 

Objective:  This  is  a  faculty  fellowship  award  for  visiting  Big  XII  universities  to 

enhance  research  collaborations. 

2.  Oklahoma  EPSCoR  —  Research  Opportunity  Award 

Piao,  D  (Supporting  PI),  06/2009— -08/2009,  OSU-  $3,500 
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Jiang,  Y  (PI)  University  of  Central  Oklahoma  (Total  $10,000) 

Title:  Advanced  Reconstruction  Method  for  Optical  and  Ultrasonic  Imaging 

Objective:  This  is  a  summer  research  opportunity  award  grant  for  faculties  in 

Oklahoma  regional  universities  to  conduct  research  in  Oklahoma 
research-intensive  universities,  including  Oklahoma  State  University. 

3.  Oklahoma  State  University — Technology  Business  Assessment  Group  (TBAG) 

Piao,  D  (PI),  12/01/2008— -11/30/2009,  $25,000 

Title:  Clinical  Trans-rectal  Optical  Imaging  Applicator  for  Prostate  Biopsy 

Objective:  To  develop  a  B&K  ultrasound-coupled  near-infrared  optical  imaging 

applicator  for  trans-rectal  near-infrared  tomography  of  human  prostate  for 
future  clinical  trials. 


CONCLUSIONS 

The  objective  of  this  research  is  to  explore  the  technology  of  trans-rectal  near-infrared  (NIR) 
optical  tomography  that  may  benefit  accurate,  selective  prostate  biopsy.  The  research  in  the  2nd 
year  has  made  important  advancements.  The  advancements  include  the  invention  of  TRUS- 
coupled  NIR  tomography,  the  development  of  a  canine  prostate  tumor  model,  the  demonstration 
of  TRUS-coupled  NIR  tomography  of  canine  prostate  tumor  in  vivo,  and  the  establishment  of 
preliminary  evidence  that  augmenting  the  NIR  imaging  to  the  standard  TRUS  detects  the 
prostate  tumor  earlier  and  with  better  accuracy.  The  proposed  technology  has  advanced  an 
important  step,  from  trans-rectal  NIR  optical  tomography  to  TRUS-coupled  NIR  optical 
tomography,  toward  clinical  validations. 
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Abstract:  We  investigate  the  feasibility  of  trans-rectal  optical  tomography 
of  the  prostate  using  an  endo-rectal  near-infrared  (NIR)  applicator  that  is  to 
be  integrated  with  a  trans-rectal  ultrasound  (TRUS)  probe.  Integration  with 
TRUS  ensures  accurate  endo-rectal  positioning  of  the  NIR  applicator  and 
the  utility  of  using  TRUS  spatial  prior  information  to  guide  NIR  image 
reconstruction.  The  prostate  NIR  image  reconstruction  is  challenging  even 
with  the  use  of  spatial  prior  owing  to  the  anatomic  complexity  of  the 
imaging  domain.  A  hierarchical  reconstruction  algorithm  is  developed  that 
implements  cascaded  initial-guesses  for  nested  domains.  This  hierarchical 
image  reconstruction  method  is  then  applied  to  evaluating  a  number  of  NIR 
applicator  designs  for  integration  with  a  sagittal  TRUS  transducer.  A  NIR 
applicator  configuration  feasible  for  instrumentation  development  is 
proposed  that  contains  one  linear  array  of  optodes  on  each  lateral  side  of  the 
sagittal  TRUS  transducer.  The  performance  of  this  NIR  applicator  is 
characterized  for  the  recovery  of  single  tumor  mimicking  lesion  as  well  as 
dual  targets  in  the  prostate.  The  results  suggest  a  strong  feasibility  of  trans- 
rectal  prostate  imaging  by  use  of  the  endo-rectal  NIR/US  probe. 

©2008  Optical  Society  of  America 
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1.  Introduction 

Prostate  cancer  is  the  second  most  commonly  diagnosed  cancer  and  the  second  leading  cause 
of  cancer  deaths  in  American  men  [1].  Prostate  cancer  screening  is  performed  by 
measurement  of  serum  prostate-specific  antigen  (PSA)  [2],  digital  rectal  examination  (DRE), 
and  in  many  cases  a  combination  of  both  tests  [3].  The  introduction  of  PSA  test  contributed  to 
substantially  increased  detection  rate  of  organ-confined  prostate  cancer  or  considerable  stage 
migration  [4].  However,  PSA  is  not  a  specific  indicator  of  prostate  malignancy  and  post¬ 
treatment  tumor  recurrence,  except  after  radical  prostatectomy  [5].  A  clearly  increased  serum 
PSA  value  (>20  ng/ml)  may  indicate  the  presence  of  a  prostate  carcinoma  at  a  very  high 
probability  [6].  In  the  gray  zone  between  4  and  10  ng/ml  the  tissue  marker  PSA  is  frequently 
influenced  by  benign  alterations,  so  that  it  is  not  possible,  on  the  basis  of  the  PSA  value  alone, 
to  differentiate  between  benign  and  malignant  cases  [6,  7].  DRE  can  often  distinguish  between 
prostate  cancer  and  non-cancerous  conditions;  it  may  also  detect  prostate  cancers  having 
normal  PSA  levels.  However,  palpation  during  a  DRE  is  subjective,  insensitive,  and  more 
than  half  of  all  prostate  cancers  detected  are  not  palpable  [3].  When  the  suspicion  of  prostate 
cancer  is  raised  by  abnormal  PSA  and/or  DRE,  the  diagnosis  is  made  by  biopsy.  The 
technique  of  trans-rectal  ultrasound  (TRUS)  based  trans-rectal  prostate  biopsy,  carried  out 
with  a  semi-automatic  coil  spring  device  and  an  18-gauge  needle,  is  to  date  considered  as  the 
gold  standard  [6] . 

Prostate  neoplastic  lesions  may  be  identified  on  TRUS  as  being  hypoechoic  [8].  However 
at  most  60%  neoplastic  lesions  appear  hypoechoic  on  TRUS  while  most  of  the  remaining 
neoplastic  lesions  appear  isoechoic  [9].  The  hypoechoic,  cancer-suspicious  areas  may  be 
histologically  either  benign  or  malignant  [9].  The  lack  of  TRUS  specificity  thereby  prompts 
the  practice  of  “systematic  biopsy”  of  the  prostate.  The  current  trend  is  to  use  10-  to  12-core 
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biopsy  with  a  preference  in  the  peripheral  zone,  where  most  neoplastic  lesions  are  found,  as 
the  initial  biopsy  strategy.  It  should  be  noted  that  the  majority  of  biopsies  are  found  to  be 
negative,  and  in  men  with  persistent  suspicion  of  prostate  cancer  after  several  negative 
biopsies,  more  extensive  protocols  (>12  cores)  up  to  saturation  biopsy  (24  cores)  represent  a 
necessary  diagnostic  procedure  [10].  However,  despite  years  of  research,  the  exact  number  of 
biopsies  to  be  taken  is  still  largely  unknown  [11]. 

The  need  of  having  many  biopsy-cores  for  systematic,  yet  random,  tissue  sampling  of  the 
prostate  may  be  alleviated  if  the  acoustic  contrast  that  TRUS  relies  on  is  augmented  with 
functional  or  “surrogate”  markers  of  the  prostate  tumor  such  that  the  biopsy  is  directed  to  the 
malignant  lesions.  A  functional  imaging  modality  augmenting  TRUS  is  certainly  more 
desirable  if  it  is  non-ionizing  and  minimally-invasive  as  is  TRUS.  Optical  tomography  based 
on  near-infrared  (NIR)  light  could  emerge  as  such  a  modality. 

Near-infrared  measurements  of  attenuation  through  tissue  have  demonstrated  significant 
contrast  gradients  between  blood  and  parenchymal  tissue  that  is  otherwise  difficult  to  obtain 
[12-17].  The  alteration  of  vascularity  or  the  hemoglobin  content  in  the  tumor  provides  high 
intrinsic  optical  contrast  between  the  tumor  and  benign  tissues  which  has  been  well- 
demonstrated  in  breast  cancer  imaging  [12-17].  When  multi-spectral  detection  is  engaged, 
NIR  imaging  is  also  capable  of  directly  quantifying  the  chromorphore  concentrations 
important  for  characterization  of  the  malignancy  [12-17].  In  prostate,  studies  have  shown 
vascular  density  gradient  in  malignant  versus  benign  tissue  specimens  [18],  and  different 
water  concentrations  in  cancerous  and  benign  tissues  in  vitro  [19].  Invasive  NIR 
measurements  of  prostate  have  been  conducted  for  experimental  prostate  tumors  [20]  and 
human  prostate  [21,  22].  Surface  measurements  of  implanted  prostate  tumor  have  also  been 
reported  [23,  24].  All  these  studies  demonstrate  the  potential  of  using  NIR  to  detect  and 
characterize  prostate  cancer.  NIR  diffuse  optical  measurement,  performed  interstitially,  is  also 
becoming  an  important  tool  for  monitoring  photodynamic  therapy  in  prostate  [21,  22,  25]. 
Prostate  NIR  imaging  via  trans-urethral  probing  had  been  analyzed  and  tested  [26].  Recently, 
trans-rectal  prostate  NIR  imaging  has  been  investigated  in  simulation  in  the  context  of 
assisting  MRI  for  treatment  decision  [27]. 

To  our  knowledge,  experimental  work  on  trans-rectal  NIR  tomography  has  not  been 
performed  except  for  our  recent  attempts  [28,  29]  which  may  be  largely  due  to  the  challenge 
of  fabricating  a  suitable  trans-rectal  applicator.  Optical  tomography  typically  needs  10s  of 
channels  of  NIR  optodes  in  order  to  achieve  reasonable  spatial  resolution  as  a  large  tissue 
volume  is  being  interrogated.  The  NIR  illumination  can  be  delivered  by  small  diameter  fibers, 
but  the  detection  of  weak  scattered  light  is  in  favor  of  fibers  of  larger  diameters  and/or 
suitably  larger  numerical  apertures.  Unlike  in  breast  NIR  tomography  where  there  is  minimum 
spatial  restriction  for  the  optode  configuration,  trans-rectal  applicator  for  NIR  tomography  has 
to  deploy  many  optodes  in  a  very  compact  space.  This  restriction  could  become  more 
pronounced  when  trans-rectal  NIR  applicator  is  also  to  be  combined  with  TRUS  transducer. 
Since  the  depth  of  tissue  interrogation  by  diffuse  NIR  light  is  roughly  one-half  of  the  source- 
detector  separation  for  typical  scattering-dominant  biological  tissue,  reaching  targets 
centimeters  deep  in  prostate  implies  a  NIR  array  of  several  centimeters  in  size.  Such  an  NIR 
array  is  feasible  for  trans-rectal  application  if  the  optodes  are  arranged  longitudinally,  which 
should  provide  a  sagittal-view  in  trans-rectal  NIR  imaging.  The  images  obtained  by  trans- 
rectal  NIR  tomography  alone  would,  however,  be  difficult  to  correlate  with  the  anatomy. 
Unlike  in  breast  imaging  where  the  NIR  applicator  could  be  accurately  positioned,  accurate 
positioning  of  a  trans-rectal  NIR  applicator  with  respect  to  the  prostate  is  difficult  which  is 
due  to  the  “blind”  location  of  the  prostate,  slow  NIR  image  reconstruction,  and  lack  of 
anatomic  details  in  NIR  tomography  images.  It  is  thereby  imperative  to  use  a  real-time 
morphological  imaging  modality  concurrently  with  trans-rectal  NIR  tomography  to  provide  a 
positioning  guidance  for  NIR  applicator  in  order  to  correlate  the  NIR  tomography  findings 
with  the  prostate  anatomy.  The  structural  information  of  the  prostate  can  further  be  utilized  as 
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the  spatial  prior  [30]  to  improve  the  accuracy  of  NIR  image  reconstruction.  Among  the 
prostate  imaging  techniques,  TRUS  is  perhaps  the  best  modality  for  trans -rectal  NIR 
tomography  to  combine  owing  to  the  operational  similarity  between  these  two  modalities. 

The  diagnostic  benefit  of  augmenting  NIR  contrast  to  US  has  been  demonstrated  in  breast 
cancer  detection  [31,  32].  The  methodology  of  combining  NIR  &  US  can  certainly  be 
extended  from  breast  imaging  to  prostate  imaging;  nevertheless,  the  technique  cannot  be 
extended  from  imaging  the  breast  to  imaging  the  prostate  without  an  NIR/US  applicator 
suitable  for  trans-rectal  manipulation.  In  this  work,  we  demonstrate  TRUS  coupled  trans- 
rectal  optical  tomography  of  the  prostate.  This  work  is  reported  in  two  consecutive  papers. 
The  Part-I  paper,  based  on  simulation,  investigates  designs  of  NIR  tomography  applicator 
suitable  for  integrating  with  a  commercial  TRUS  transducer.  A  hierarchical  NIR  image 
reconstruction  algorithm  is  developed  for  utilizing  the  TRUS  structural  a  priori  information 
which  is  then  used  to  evaluate  several  NIR  applicator  configurations  for  integrating  with  a 
TRUS  transducer.  The  Part-II  paper  implements  the  probe  design  suggested  in  Part-I,  presents 
the  instrumentation  details  of  the  TRUS  coupled  trans-rectal  NIR  tomography  probe  & 
system,  and  demonstrates  TRUS-coupled  trans-rectal  NIR  tomography  of  a  canine  prostate. 

2.  The  geometry  of  trans-rectal  NIR  imaging  and  the  utility  of  TRUS  information 

2.7  Configuration  of  sagittal  trans-rectal  NIR  array  for  coupling  with  sagittal  TRUS 

The  TRUS  is  typically  performed  in  bi-plane  (sagittal  and  transverse)  for  prostate  imaging. 
The  sagittal  and  transverse  views  are  switched  during  prostate  imaging,  but  for  biopsy 
procedure,  the  firing  of  the  spring-loaded  needle  is  monitored  in  the  sagittal  plane  wherein  the 
needle  trajectory  may  be  accurately  marked.  We  have  acquired  a  bi-plane  TRUS  probe  (Aloka 
UST-672-5/7.5)  as  is  shown  in  Fig.  1(a),  which  has  a  proximal  sagittal  transducer  window  of 
60mmxl0mm  and  a  distal  transverse  transducer  window  of  120°xl0mm.  The  cylindrical 
TRUS  probe  has  a  maximum  diameter  of  20mm.  Integrating  NIR  applicator  to  a  TRUS  probe 
implies  that  the  dimension,  particularly  the  radial  one,  of  the  NIR  array  is  quite  restricted.  As 
discussed  earlier,  longer  source-detector  separation  is  needed  to  interrogate  deeper  targets; 
therefore  the  most  feasible  configuration  of  an  NIR  array  may  be  through  distributing  a  linear 
array  of  optode  on  each  lateral  side  of  the  sagittal  TRUS  transducer,  as  illustrated  in  Fig.  1(b). 
Taking  into  account  the  mechanical  structure  necessary  to  support  the  optodes,  the  combined 


Sagittal  transducer 


Transverse 

transducer 


O  x  are  optodes 


(b) 


(c) 


(d) 


Fig.  1.  (a)  Photograph  of  a  bi-plane  TRUS  transducer,  (b)  Configuration  of  an  NIR  array 
feasible  for  coupling  with  sagittal  TRUS.  (c)  The  NIR  imaging  geometry  for  the  one  depicted 
in  (b).  (d)  Illustration  of  a  trans-rectal  probe  that  integrates  NIR  and  US  synergistically. 
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probe  will  likely  have  a  radial  size  of  at  least  a  few  millimeters  larger  than  that  of  the  original 
TRUS  probe.  NIR  tomography  has  to  use  multimode  fibers  to  detect  weak  diffuse  light,  and 
these  multimode  fibers  must  be  delivered  longitudinally  before  being  side-fired.  Bending  the 
fiber  is  not  a  viable  solution  here  for  side-firing  unless  the  fibers  are  passed  inside  the  TRUS 
probe.  The  side-firing  alternatively  may  however  be  realized  by  implementing  micro-optical 
components,  but  the  compact  space  inside  or  surrounding  the  probe  may  not  accommodate  a 
large  number  of  such  channels.  A  NIR  array,  which  could  couple  with  TRUS,  may  be 
configured  by  fabricating  7  fiber  channels  on  each  lateral  side  of  the  TRUS  to  span  60mm 
longitudinally  as  the  TRUS  does,  and  placing  the  optical  channels  10mm  apart.  The  NIR  array 
must  leave  the  lOmm-wide  sagittal  TRUS  transducer  unblocked;  therefore  a  20mm  spacing  of 
the  NIR  optodes  from  one  lateral  side  to  the  other  is  perhaps  needed.  These  considerations 
lead  to  the  NIR  array  geometry  shown  in  Fig.  1(c)  where  14  optodes  are  spaced  10mm 
longitudinally  and  20mm  laterally.  Fig.  1(d)  illustrates  an  NIR/US  probe  if  the  sagittal  NIR 
array  can  be  fabricated  synergistically  with  the  TRUS  probe. 

2.2  Forward  and  inverse  methods  for  sagittal  trans -rectal  optical  tomography 

The  prostate  and  peripheral  tissues  are  scattering-dominant  in  NIR  [19-22].  We  use  the 
diffusion  approximation  to  the  radiative  transport  equation  in  frequency-domain  [33]: 

V  •  D(r)VU(r,  a>)  -  (jua  +  —)U(r,a)  =  -S(f,aj)  0) 

c 

where  U(r,co)  is  the  photon  fluence  rate  at  position  f ,  S(r,cd)  is  the  source,  CO  is  the 
source  modulation  frequency,  c  is  the  speed  of  light  in  the  medium,  jia  is  the  absorption 
coefficient,  and  d  =  [ 3(p  +//’)]  1  *s  the  diffusion  coefficient  with  p  being  the  reduced  or 

transport  scattering  coefficient.  Finite-element  method  [34]  is  used  to  solve  Equ.  (1)  under  the 
Robin-type  boundary  condition  [33]: 

U(r0,co)  +  2DAh0  •VU(r0,a))  =  0  (2) 

where  rQ  is  the  boundary,  h{)  is  the  outward  normal  vector  of  the  boundary  and  A  is  the 

refractive  index  mismatch  coefficient.  The  refractive  indices  used  for  air  and  tissue  are  1  and 
1.33  respectively,  leading  to  A=2.82  as  in  [33]. 

The  imaging  volume  is  divided  to  4  regions-of-interest  (ROIs):  the  rectum  wall,  the  peri- 
prostate  tissue,  the  prostate,  and  the  prostate  tumor.  The  Jacobian  (sensitivity)  values  are 
calculated  for  each  ROI  rather  than  each  node  which  has  the  form  of: 
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where  I..  (i,j  =  1,2,. ..,7)  and  (/>..  (i,j  =  1,2,. ..,7)  are  the  intensity  and  phase  terms  of 
U(r,C0 ),  respectively.  In  Equ.  3,  “rect”,  “peri”,  “pros”,  and  “lesi”  denote  “rectum  wall”, 
“peri-prostate  tissue”,  “prostate”,  and  “prostate  lesion”,  respectively. 

The  Levenberg-Marquart  (LM)  algorithm  [35]  governs  the  iterative  recovery  of  the  optical 
properties  by  updating  the  ROI-specific  values  of  jua  and  jfs  according  to 

xk+ 1  =  xk  +  a  -  [JT (xk)J(xk)  +  Aiy1  JT (xk)Av(xk)  (4) 

where  x  is  the  array  of  parameters  to  be  optimized,  Av  is  the  forward  projection  error  and  X 
is  a  penalty  or  regularization  term.  A  small  damping  factor  (X  in  the  range  of  (0,  1)  is 
introduced  in  Equ.  4  to  stabilize  the  convergence.  It  is  shown  that  an  empirically  chosen  (X 
could  make  the  LM  algorithm  more  reliable  and  computationally  more  efficient  [36]. 


2.3  TRUS  prior  assisted  finite- element  mesh  for  trans -rectal  NIR  tomography  reconstruction 


The  TRUS  prostate  images  which  are  available  in  open  sources  [37]  are  used  for  the 
simulation  study.  The  TRUS  image  was  first  imported  into  a  pre-processing  software  3ds- 
MAX  [Autodesk  Inc]  (shown  in  Fig.  2(a)).  The  3ds-MAX  provides  very  flexible  geometry¬ 
deforming  functions,  with  which  a  basic  3-D  geometry  of  the  prostate  can  be  outlined 
manually.  The  finalized  3-D  mesh  of  the  prostate  is  then  converted  to  COMSOL  Multiphysics 
[COMSOL  AB]  compatible  format  (shown  in  Fig.  2(b))  using  MeshToSolid  [Syncode  Inc]. 
The  prostate  tumor  is  then  mimicked  using  a  spherical  shape  to  allow  for  flexibility  of 
adjusting  its  size.  The  absorption  and  reduced  scattering  coefficients  of  the  rectum,  peri- 
prostate  tissue,  prostate  and  the  tumor  are  assigned  with  values  as  suggested  by  literature  [27]. 
Figure  3  illustrates  one  example  of  the  completed  FEM-mesh  for  trans-rectal  optical 
tomography  derived  from  a  TRUS  image,  where  x,  y  and  z  denoting  the  longitudinal,  lateral 
and  the  depth  coordinates,  respectively.  A  typical  mesh  used  for  this  work  contains 
approximately  4000  nodes  and  20000  linear  tetrahedral  elements. 


Mesh  To  Solid  3.0  -  [ 


&  Fie  Edit  View  Measure  Tools 


Cfck  and  drag  the  mouse  to  pan  over  the  mode  ^ 


(a) 


(b) 


Fig.  2.  (a)  3ds-MAX  Interface,  (b)  Mesh-to-Solid  Interface 


The  mesh  in  Fig.  3  corresponds  to  a  volume  of  80x80x80  mm3  and  the  ‘walnut’  shaped 
prostate  has  a  dimension  of  50x50x30  mm3.  The  rectum  wall  is  4mm  thick  with  a  curvature 
radius  of  80  mm.  The  choice  of  the  curvature  radius  is  due  to  the  fact  that  the  NIR  array  added 
to  a  TRUS  probe  may  have  a  flat  surface  that  would  transform  the  rectum  lumen  to  an 
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elliptical  shape.  A  larger  radius  also  gives  more  flexibility  in  handling  the  posterior  prostate 
region  within  the  mesh. 

Although  only  the  rectum  wall  is  a  physical  boundary,  treating  the  other  5  surfaces  as 
physical  boundaries  (Robin  type,  Equ.  2)  should  have  negligible  effect  upon  the  results  as  the 
lateral-medial  and  ventral-dorsal  dimensions  well  exceed  the  potential  path  of  photon 
propagation  for  the  NIR  array  given  in  Fig.  1(c).  The  modulation  frequency  CO  of  the  source 
in  Equ.  1  is  set  at  100MHz,  and  1%  Gaussian  noise  is  added  to  all  forward  calculations  to 
form  the  measurement  data. 


(a)  Transverse  View  (b)  Sagittal  View  (c)  Coronal  View  (d)  Perspective  View 

Fig.  3.  FEM  mesh  generated  based  on  approaches  in  Fig.  2.  [Unit:  mm]. 

3.  A  hierarchical  spatial  prior  approach  for  trans-rectal  NIR  tomography  reconstruction 

3.1  Sensitivity  of  the  sagittal  trans-rectal  NIR  array 

The  NIR  array  proposed  in  Fig.  1(c)  has  7  source  channels  occupying  one  linear  array  and  7 
detection  channels  occupying  the  other  array.  The  sensitivity  with  respect  to  a  perturbation  of 
a  specific  optical  property  is  determined  by  the  corresponding  Jacobian  values  in  Equ.  3. 
Figure  4  plots  the  sensitivity  specific  to  absorption,  or  d  in  I..  jdju  ,  for  a  medium  with  optical 

properties  of  jua  =  0.01mm-1  and  jus  =  1.0mm-1  ,  calculated  by  projecting  the  Jacobian  values 

along  a  line  in  the  imaging  volume.  Figure  4(a)  is  the  longitudinal  sensitivity  in  the  mid- 
sagittal  plane  for  a  line  from  (0,  40,  30)  to  (80,  40,  30),  Fig.  4(b)  is  the  lateral  sensitivity  in  the 
mid-transverse  plan  for  a  line  from  (40,  0,  30)  to  (40,  80,  30),  and  Fig.  4(c)  is  the  depth 
sensitivity  in  the  mid-sagittal  plane  for  a  line  from  (40,  40,  15.1)  (here  15.1  is  the  z 
coordinate,  but  the  actual  depth  from  the  rectum  surface  is  0mm  owing  to  the  curvature  of  the 
rectum)  to  (40,  40,  80),  respectively.  The  dimension  or  the  locations  of  the  source  &  detector 
array  is  marked  on  the  abscissa  of  all  three  plots.  The  TRUS  sagittal  plane  is  located  at 
y=40mm,  which  is  the  mid-sagittal  plane  within  the  NIR  imaging  volume. 


Fig.  4.  Sensitivity  profile,  (a)  Mid-sagittal  plane,  longitudinal  sensitivity;  (b)  mid-transverse 
plane,  lateral  sensitivity;  (c)  mid-saggital  plane,  depth  sensitivity.  The  marks  on  abscissa  or  the 
origin  show  the  positions  of  optodes. 


Figure  4  indicates  that  the  longitudinal  sensitivity  has  ~6dB  variation  in  the  middle  75% 
range  of  the  array,  and  the  lateral  sensitivity  peaks  at  the  mid-saggital  plane.  In  the  middle- 
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sagittal  plane  the  sensitivity  degrades  ~ldB/mm  as  z-coordinate  increases  from  20mm,  which 
is  apparently  due  to  the  side-way  placement  of  the  NIR  array.  The  depth-degrading  sensitivity 
will  cause  deeper  targets  to  be  reconstructed  at  a  shallower  position  [28]  if  no  spatial  prior  is 
incorporated. 

3.2  Trans-rectal  NIR  image  reconstruction  without  a  priori  information 

The  performance  of  recovering  tumor-mimicking  target  by  trans-rectal  NIR  tomography 
without  any  structural  prior  is  examined.  Figure  5  lists  the  results  for  the  tumor  target  being 
placed  at  left,  middle,  and  right  within  the  prostate.  The  top  row  in  Fig.  5  lists  the  target 
images  of  jua  and  jus  generated  by  the  TRUS-defined  geometry  as  shown  in  Fig.  3.  The 

optical  properties  of  the  10mm  diameter  tumor  target  are  pa  =0.02 mm~l  and  ju'  =1.6 mm~\ 

with  the  parameters  of  other  regions  as  listed  in  Table  1.  These  target  images  are  used  to 
generate  the  noise-added  simulated  measurement  data.  The  image  reconstruction  is  then 
conducted  using  a  mesh  of  homogenous  element  density  throughout  the  entire  volume  and 
updated  element  by  element.  The  results  are  given  in  the  bottom  row  of  Fig.  5.  As  seen,  the 
tumor  target  may  be  localized,  but  the  recovered  resolution  or  spatial  information  is  poor.  The 
accuracy  of  optical  property  recovery  is  also  low.  Further,  it  is  found  that  a  tumor  target  with 
negative  contrast  in  absorption  cannot  be  accurately  recovered  using  similar  settings.  It  is 
however  expected  that  the  detection  and  characterization  of  the  tumor  target  will  be  improved 
when  the  spatial  information  of  prostate  and  the  tumor  is  available. 


(a)  (b) 


Fig.  5.  NIR-only  element-based  reconstruction  of  a  tumor  target  in  various  longitudinal 
locations.  Row  1:  target  image  for  calculating  the  forward  data;  Row  2:  Images  reconstructed 
without  any  spatial  prior,  (a):  na  images;  (b):  jus  images.  [Unit:  mm'1] 

3.3  A  hierarchical  spatial  prior  method  for  TRUS  guided  trans-rectal  NIR  reconstruction 

The  tissue  volume  interrogated  by  trans-rectal  NIR  imaging  constitutes  a  nested-domain 
including  a  thin  layer  of  rectum  wall,  a  large  volume  of  peri-prostate  tissue,  a  relatively 
absorbing  prostate,  and  the  lesion  within  the  prostate.  These  nested  imaging  domains  may  be 
further  complicated  by  the  pelvic  bone  that  could  interfere  with  the  light  propagation. 
Schweiger,  et  al.  [38],  Kolehmainen  et  al.  [39-41],  and  Srinivasan  et  al.  [42]  have  previously 
investigated  the  issue  of  recovering  the  shapes  and  optical  properties  of  regions  with  optical 
contrast  inside  a  non-nested  or  nested  domains,  where  the  shapes  of  the  ROIs  were  derived 
from  optical  information  when  no  spatial  prior  is  available  from  other  complementary 
imaging  modalities.  These  methods  have  shown  sufficient  robustness  in  recovering  the  shapes 
and  optical  properties  of  the  ROIs,  yet  the  problem  of  stability  and/or  slow  convergence  was 
noticed  in  such  approaches  dealing  with  nested-domains.  In  trans-rectal  NIR  tomography 
reconstruction  the  spatial  information  from  TRUS  may  be  implemented  by  assigning 
homogenous  optical  properties  within  each  ROIs  of  the  imaging  domain.  However  the 
convergence  and  the  accuracy  of  reconstruction  will  still  depend  upon  the  initial  guess  in 
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addition  to  the  accuracy  of  the  prior  information.  The  dependence  on  initial  guess  in  a 
gradient  based  solver  is  due  to  the  local  minimum  feature  [36],  as  indicated  in  Fig.  6,  which 
could  be  exaggerated  in  prostate  imaging  due  to  the  possible  multiple  combinations  of  optical 
properties  in  the  nested-structures.  The  image  reconstruction  in  trans-rectal  optical 
tomography  is  further  complicated  by  the  discrepancy  regarding  the  optical  contrast  that  the 
prostate  tumor  could  have,  namely  positive  or  negative  [19-22,  27]. 


Fig.  6.  Local-minimum  issue  in  reconstruction.  The  forward  calculation  is  based  on  Fig.  3(a) 
with  assumption  of  homogeneous  imaging  volume.  The  projection  error  is  calculated  by  using 
reduced  scattering  coefficient  of  the  true  value  (0.008mm1),  and  using  an  absorption 
coefficient  value  from  0  to  0.15mm'1  at  a  step  of  0.002mm'1,  with  respect  to  the  forward  data. 
Other  than  the  global  minimum,  three  local  minimums  can  be  observed  where  the  iteration  can 
stop  incorrectly.  This  is  the  effect  of  varying  only  one  parameter.  More  local  minimums  may 
occur  when  reconstructing  more  parameters. 


When  the  TRUS  is  available,  a  conventional  method  of  utilizing  the  spatial  information 
would  be  the  having  the  optical  property  of  each  ROI  set  as  homogenous  and  updated 
simultaneously  at  each  iteration.  However,  we  have  found  that  this  conventional  approach 
may  not  lead  to  reliable  convergence  for  prostate  imaging,  which  is  attributed  to  the  local- 
minimum  problem.  One  example  is  given  in  Table  1  for  the  NIR  array  shown  in  Fig.  1(c).  The 
prostate  model  is  generated  according  to  a  previous  work  [27]  (details  of  which  are  given  later 
in  Fig.  8),  and  a  target  of  10mm  diameter  is  located  at  (40,  50,  15)  which  is  15  mm  from  the 
rectal  surface.  When  the  four  ROIs  including  the  rectum,  the  peri-prostate  tissue,  the  prostate, 
and  the  prostate  tumor  are  updated  simultaneously  from  the  same  initial  guess  of  jua 

=0.0 1mm"1  and  p  =  1.0mm"1,  the  iteration  stops  after  1  update  due  to  the  negative  pa  value 

obtained  for  the  rectum  wall.  The  iteration  fails  to  continue  apparently  due  to  the  local 
minimum  issue. 


Table  1 .  Results  of  simultaneously  updating  the  4  ROIs  from  the  same  initial  guess 


(mm1) 

t*s  (mm’1) 

Regions 

Surrounding 

Tissue 

Rectum 

Wall 

Prostate 

Tumor 

Surrounding 

Tissue 

Rectum 

Wall 

Prostate 

Tumor 

Set  value 

0.002 

0.01 

0.06 

0.02 

0.8 

1 

1.27 

1.6 

Simultaneous 

Update 

0.1216 

-0.008 

0.026 

0.0215 

1.1482 

2.3602 

0.6173 

0.7073 

The  local  minimum  problem  may  be  mitigated  by  a  cascaded  initial-guess  approach  or  a 
hierarchical  spatial  prior  method.  The  principle  of  this  method  is  to  first  reconstruct  the  global 
optical  properties  of  the  entire  volume,  then  to  reconstruct  the  optical  properties  of  prostate 
and  rectum  wall,  and  last  to  reconstruct  the  tumor  lesion  area.  The  2nd  and  3rd  steps  use  the 
value  obtained  in  the  previous  step  as  the  initial  guess  of  that  specific  ROI.  Therefore  at  each 
step,  the  perturbation  by  a  relatively  smaller  region  is  less  influential  and  convergence  of  the 
iteration  is  better  achieved.  The  detailed  steps  are  shown  in  Fig.  7  and  described  in  below: 
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(a)  The  first  iterations  assign  an  entirely  homogenous  imaging  volume.  In  this  round  the 
initial  projection  error  will  be  large  and  the  convergence  is  most  likely  dominated  by  the 
global  minimum.  A  single  set  of  jua  and  ps  are  determined  using  LM  algorithm  (Equ.  4)  and 

will  be  used  as  the  initial  guess  for  the  second  step. 

(b)  The  second  iterations  consider  three  regions,  the  rectum  wall,  peri-prostate  tissue,  and 
prostate,  within  the  imaging  volume.  The  calculation  of  the  optical  properties  of  these  three 
ROIs  start  at  the  same  initial  guess  as  provided  in  step  (a)  but  converges  at  different  values. 

(c)  The  values  obtained  from  step  (b)  are  used  as  the  initial  guess  for  the  three  ROIs  but 
with  a  tumor  added  to  the  prostate.  The  tumor  and  the  prostate  take  the  same  initial  values  as 
determined  by  the  previous  step.  Each  of  the  four  ROIs  (rectum  wall,  peri-prostate  tissue, 
prostate,  and  tumor)  converge  to  different  end  value. 

The  change  of  the  overall  projection  error  for  the  three  steps  is  plotted  in  Fig.  7(d),  and 
as  evident,  rapid  and  reliable  convergence  is  observed.  The  hierarchy  of  the  implementation  of 
initial  values  for  iteration  is  illustrated  alternatively  in  Fig.  7(e). 


0  0.02  0.04  0.06  0.8  1  1.2  1.4  1.6 


(b) 


Step  1  Step  2  Step  3 


(e) 


Fig.  7.  The  3-step  hierarchical  reconstruction  method  (a)  Step  1 — one  ROI  for  the  entire 
volume;  (b)  step  2 — three  ROIs  representing  rectum  wall,  peri-prostate  tissue,  prostate;  (c)  step 
3 — four  ROIs  representing  rectum  wall,  peri-prostate  tissue,  prostate,  tumor;  (d)  change  of  the 
overall  projection  error,  where  the  dash  lines  separate  the  converging  of  the  three  steps  in  (a) — 
(c);  (e)  block  chart  of  the  hierarchical  initial  guess  assignment.  [Unit  in  (a)— (c):  mm'1] 


3.4  Validation  of  the  hierarchical  spatial  prior  method 

Recently  Li  et  al.  reported  simulation  results  for  trans-rectal  optical  tomography  in  the  context 
of  using  MRI  anatomic  information  [27],  which  is  referred  to  as  “NIR/MRI”  in  the  following 
text.  The  proposed  hierarchical  spatial  prior  method  is  evaluated  using  the  same  probe 
geometry,  prostate  geometry  and  the  optical  properties  (also  in  Table  1)  presented  in  the 
NIR/MRI  work.  The  size  and  depth  of  the  tumor  for  simulation  were  not  specified  in  the 
NIR/MRI  work,  but  a  tumor  with  diameter  of  10  mm  and  a  depth  of  15mm  from  the  planar 
probe  surface  is  considerably  close  to  the  one  presented  in  the  NIR/MRI  paper.  The  NIR/MRI 
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work  also  indicated  the  challenge  of  reconstruction  that  may  be  due  to  the  local  minimum. 
The  authors  set  an  arbitrary  searching  range  for  the  optical  properties  (jua  :0-0.1  mm~\jus  :0- 

2  mm-1),  and  in  4  sets  of  the  results,  3  of  the  tumor  absorption  value  reached  the  limits  and 
were  stopped  from  further  iteration,  whereas  the  4th  value  converged  at  a  value  more  than  2 
folds  of  the  set  value. 


Fig.  8.  The  FEM  mesh  generated  by  following  the  geometry  in  NIR/MRI  paper:  (a)  the 
geometry  of  the  optodes,  where  the  dash  rectangle  delineates  the  dimension  of  the  NIR  array 
that  will  be  evaluated  later  for  integration  with  US;  (b)  the  3-d  view  of  the  optodes  and  the 
imaging  volume;  (c)  FEM  containing  the  prostate  and  the  tumor.  [Unit:  mm] 


In  the  NIR/MRI  work,  a  stand-alone  trans-rectal  NIR  probe  is  simulated.  The  stand-alone 
NIR  probe  could  have  contained  more  optode  channels  incorporated  than  a  TRUS-coupled 
NIR  probe.  In  the  NIR/MRI  work,  the  best  result  is  deducted  for  10  sources  and  28  detectors, 
which  is  used  to  evaluate  the  hierarchical  method.  The  NIR  probe  geometry  and  the  imaging 
domain  of  the  NIR/MRI  work  are  re-plotted  in  Fig.  8  for  clarification  and  the  hierarchical 
method  is  also  preformed  in  transverse-view  as  did  the  NIR/MRI  work.  A  1%  noise  is  also 
used  in  both  methods  assuring  the  consistency  of  measurement  data. 


Table  2.  Reconstruction  of  a  prostate  tumor  of  negative  contrast  with  respect  to  the  prostate 


(mm1) 

Vs  (rnm1) 

Regions 

Surrounding 

Tissue 

Rectum 

Wall 

Prostate 

Tumor 

Surrounding 

Tissue 

Rectum 

Wall 

Prostate 

Tumor 

Set  value 

0.002 

0.01 

0.06 

0.02 

0.8 

1 

1.27 

1.6 

NIR/MRI 

0.0025 

0.01 

0.0575 

0.0448 

0.8324 

1 

1.339 

1.075 

3- step 

0.002 

0.0099 

0.06 

0.0208 

0.8012 

1.0028 

1.2824 

1.3495 

Table  3.  Reconstruction  of  a  prostate  tumor  of  positive  contrast  with  respect  to  the  prostate 


Va  (mm1) 

Vs  (mm1) 

Regions 

Surrounding 

Tissue 

Rectum 

Wall 

Prostate 

Tumor 

Surrounding 

Tissue 

Rectum 

Wall 

Prostate 

Tumor 

Set  value 

0.002 

0.01 

0.006 

0.02 

0.8 

1 

1.27 

1.6 

3- step 

0.0020 

0.0100 

0.0061 

0.0163 

0.7998 

0.9997 

1.2863 

1.2434 

Table  2  lists  the  results  of  the  hierarchical  method  in  comparison  with  those  given  in 
NIR/MRI  paper.  The  hierarchical  method  (listed  as  “3-step”  in  the  table),  as  expected,  slightly 
outperforms  the  NIR/MRI  method  in  terms  of  the  accuracy  of  recovering  optical  properties. 
The  results  of  recovering  a  target  of  positive  absorption  contrast  are  listed  in  Table  3  using  the 
NIR/MRI  probe  geometry  and  our  proposed  3-step  method.  The  case  of  reconstructing  a 
target  with  positive  absorption  contrast  is  not  presented  in  the  NIR/MRI  paper,  therefore  only 
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the  3 -step  method  is  presented  in  Table  3  for  comparison  with  the  set  values.  In  Table  3  the 
absorption  coefficient  of  prostate  is  set  much  lower  than  that  in  Table  2  but  the  tumor  optical 
properties  are  kept  the  same  as  those  in  Table  2.  It  is  found  that  if  the  absorption  of  prostate  in 
Table  3  is  kept  the  same  as  in  Table  2,  the  positive  absorption  target  can  hardly  be 
reconstructed.  The  choice  of  lower  prostate  absorption  is  for  testing  our  hierarchical  method 
and  the  values  may  be  much  lower  than  the  values  reported  of  prostate  [19-22].  It  is  noted  that 
the  reported  values  of  prostate  absorption  coefficient  vary  in  literatures  where  the 
measurements  were  taken  either  from  in  vitro  tissue  or  from  in  vivo  tissue  using  invasive 
methods.  The  absorption  coefficient  of  intact  or  native  prostate  is  in  fact  unknown,  and  is 
likely  to  be  lower  than  the  values  reported  in  literatures. 

The  reconstructed  images  for  targets  of  both  negative  and  positive  contrasts  are  listed  in 
Fig.  9.  These  results  demonstrate  the  capability  of  our  hierarchical  spatial  prior  method  in 
reconstructing  prostate  lesion  with  either  negative  or  positive  absorption  contrast. 


0.02  0.04  0.06  1  1.2  1.4  1.6  0.005  0.0125  0.02  1  1.2  1.4  1.6 


0.02  0.04  0.06  1  1.2  1.4  1.6  0.005  0.0125  0.02  1  1.2  1.4  1.6 

(a)  (b) 


Fig.  9.  Reconstructed  images  for  a  target  with  absorption  contrast:  (a)  negative  contrast  (b) 
positive  contrast.  Top  row:  target  setting;  Bottom  row:  reconstructed  image.  [Unit:  mm'1] 
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Fig.  10.  NIR  array  designs.  The  dimension  of  the  sagittal  TRUS  is  shown  in  (a).  The  number 
of  the  Fig.  caption  denotes  the  number  of  opdotes  on  each  lateral  side  of  the  TRUS  as  is 
depicted  in  (a),  “sd”  denotes  one  line  of  source  and  one  line  of  detector;  “ssdd”  denotes  two 
lines  of  source  and  two  lines  of  detector;  “sdsd”  denotes  mixed  source/detector  in  one  line; 
“sym”  denotes  symmetric  distribution  of  the  optodes  with  respect  to  the  sagittal  TRUS.  [Unit: 
mm] 
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4.  Assessment  of  NIR  applicator  designs  for  coupling  NIR  with  TRUS 

It  is  stated  previously  that  an  NIR  array  of  dual-line  geometry  is  feasible  for  concurrent  trans- 
rectal  NIR/US  imaging  considering  the  space  limitation  when  coupling  NIR  with  TRUS 
transducer  for  endo-rectal  application.  Based  on  the  fabrication  constraints,  we  have  also 
suggested  that  each  line  array  consist  of  7  channels.  The  7  channels  could  be  exclusively 
source  or  detector  as  shown  in  Fig.  10(a),  or  interspersed  source/detector  as  in  Fig.  10(b). 
There  are  certainly  a  number  of  NIR  geometries  that  can  be  coupled  to  TRUS  sagittal 
transducer  if  not-limited  by  difficulties  in  fabrication  or  endo-rectal  use.  Compared  with  the 
geometry  in  Fig.  10(a),  more  channels  could  be  added  to  each  line-array  as  shown  in  Fig. 
10(f),  more  lines  can  be  added  as  shown  in  Fig.  10(c),  or  more  lines  and  more  channels  added 
as  in  Fig.  10(h).  More  options  are  also  listed  in  Fig.  10. 

The  array  in  Fig.  10(a)  is  the  most  desirable  in  terms  of  the  fabrication  easiness  and  endo¬ 
rectal  applicability.  The  designs  in  Fig.  10(a),  (b),  (f)  and  (g)  correspond  to  an  NIR  probe  with 
a  minimum  lateral  dimension  of  20mm.  The  designs  in  Fig.l0(c)-(e)  and  (h-j)  correspond  to 
NIR  probe  with  a  minimum  40mm  lateral  dimension  which  is  not  suitable  for  endo-rectal  use. 
The  geometries  in  Fig.  10(a) — (e)  represent  a  10mm  spacing  between  the  closest  optodes,  and 
the  geometries  in  Fig.  10(f) — (j)  require  a  5mm  spacing  between  the  closest  optodes.  The 
smaller  spacing  in  Fig.  10(f) — (j)  will  be  challenging  for  fabrication  considering  the  number 
of  fiber  channels  and  the  side-firing  configuration  if  the  probe  is  to  be  integrated  to  TRUS 
probe  unless  the  internal  structural  of  the  TRUS  probe  can  be  altered. 

4.1  Sensitivity  comparison 

The  sensitivities  of  all  the  10  configurations  of  Fig.  10  are  compared  in  Fig.  11.  The  mesh  in 
Fig.  3  with  homogeneous  optical  properties  {(ia-  0.01mm-1  and  jus  =  1.0mm-1)  is  used  for  the 

sensitivity  calculation.  The  sensitivity  in  Fig.  10  is  evaluated  at  the  lines  identical  to  those  in 
Fig.  4.  Only  the  absorption  sensitivity  is  evaluated. 

The  observations  made  from  Fig.  11  are:  (1)  increasing  the  spatial  dimension  of  source- 
detector  array  generally  improves  the  sensitivity;  (2)  increasing  the  number  of  source-detector 
pairs  generally  improves  the  sensitivity,  as  demonstrated  previously  [43];  (3)  interspersed 
source-detector  layout  may  have  slightly  wider  lateral  sensitivity  but  is  comparable  to  non- 
dispersed  source-detector  layout  for  other  imaging  views;  (4)  the  geometry  of  26ssdd  (the 
upper  thicker  line)  has  the  best  sensitivity  feature  among  the  10  geometries,  therefore  it  can  be 
used  as  a  standard  to  evaluate  the  simple  geometry  of  7sd  (the  lower  thicker  line). 

4.2  Comparison  between  the  7 sd  design  and  the  26s sdd  design 

The  7sd  design  represents  an  array  which  is  less  challenging  in  fabrication  and  more  practical 
for  endo-rectal  use.  The  26ssdd  geometry  is  impractical  for  endo-rectal  application,  difficult 
to  fabricate,  but  has  the  best  performance  among  the  designs  listed.  It  is  shown  in  Fig.  1 1  that 
the  sensitivity  of  7sd  design  is  approximately  lOdB  lower  than  that  of  26ssdd  in  the  specified 
longitudinal,  lateral,  and  depth  directions.  Quantitative  comparison  of  the  performance  is 
conducted  between  these  two  geometries  for  representative  target  variations.  The  optical 
properties  listed  in  Table  2  (set  values)  are  used  for  these  comparisons. 
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Fig.  11  Sensitivity  comparison:  (a)  longitudinal  direction,  (b)  lateral  direction,  (c)  depth 
direction.  The  upper  thicker  line  corresponds  to  the  configuration  (h)  in  Fig.  10,  and  the  lower 
thicker  line  corresponds  to  the  desired  configuration  (a)  in  Fig.  10. 


4.2.1  Reconstruction  accuracy  versus  target  longitudinal  location 

A  target  of  10mm  in  diameter  is  placed  at  the  middle-sagittal  plane  of  y=40mm,  z=26mm,  and 
varied  in  longitudinal  coordinates  from  x=25  to  55mm  with  a  step  of  5mm  (Fig.  12(a)).  The 
optical  properties  reconstructed  by  the  two  geometries  are  plotted  versus  the  true  values  in 
Fig.  12(b)  and  (c).  The  optical  properties  recovered  by  26ssdd  and  7sd  designs  are  close  to 
each  other  at  most  of  the  longitudinal  locations,  but  the  7sd  design  shows  a  larger  variation  in 
the  recovered  absorption  contrast  at  x=30mm  and  x=50mm  compared  to  other  positions.  This 
may  be  related  to  fewer  source-detector  pairs  that  contribute  to  the  target  detection  when  close 
to  the  boundary  or  the  existence  of  any  irregular  elements  in  the  mesh. 
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Fig.  12.  Comparison  of  two  geometries  for  a  target  varying  in  longitudinal  location  in  the 
middle- sagittal  plane:  (a)  illustration  of  the  target  location  change  [Unit:  mm'1];  (b)  comparison 
of  absorption  coefficient  reconstruction;  (c)  comparison  of  reduced  scattering  coefficient 
reconstruction 


4.2.2  Reconstruction  accuracy  versus  target  depth 


A  target  of  10mm  in  diameter  is  placed  at  the  middle-sagittal  plane  at  x=40mm,  y=40mm,  and 
the  depth  is  varied  from  z=25  to  40mm  at  a  step  of  2.5mm  (the  last  data  point  is  simulated  at 
z=39mm  as  at  40mm  the  target  is  out  of  the  prostate)  (Fig.  13(a)).  The  reconstructed  optical 
properties  are  plotted  in  Fig.  13(b)  and  (c).  The  26ssdd  configuration  outperforms  the  7sd  one 
again,  however  beyond  z=30mm,  both  designs  are  incapable  of  recovering  the  absorption 
coefficient  of  the  target  from  the  prostate  background.  This  depth  limitation  is  related  to  the 
maximum  span  of  the  NIR  array,  the  absorption  coefficient  of  the  prostate,  and  the  size  of  the 
target.  For  a  larger  target  with  a  diameter  of  14mm,  it  is  verified  that  the  target  may  be 
resolved  up  to  36mm  depth  from  the  NIR  array  in  comparison  to  30mm  for  a  target  of  10mm 
diameter.  A  potentially  smaller  absorption  coefficient  of  intact  prostate  may  also  increase  the 
depth  limit  of  target  detection  due  to  an  increase  in  sensitivity. 


(a) 


(b) 


(c) 


Fig.  13.  Comparison  of  two  geometries  for  a  target  varying  in  depth  in  the  middle- sagittal 
plane:  (a)  illustration  of  the  target  location  change  [Unit:  mm'1];  (b)  comparison  of  absorption 
coefficient  reconstruction;  (c)  comparison  of  reduced  scattering  coefficient  reconstruction 


4.2.3  Reconstruction  accuracy  versus  target  size 

A  target  is  placed  at  middle-sagittal  plane  of  x=40mm,  y=40mm  and  z=26mm,  and  the 
diameter  is  varied  from  4mm  to  14mm  with  a  step  of  1mm.  The  target  diameter  change  is 
illustrated  in  Fig.  14(a).  The  reconstructed  optical  properties  are  shown  in  Fig.  14(b)  and  (c). 
It  is  clear  that  the  larger  the  target,  the  better  the  accuracy  of  reconstruction.  The  26ssdd  can 
recover  the  absorption  contrast  of  the  target  when  the  diameter  is  greater  than  6mm  and  the 
7sd  can  recover  the  target  for  target  diameter  greater  than  8mm. 

These  comparisons  suggest  that  the  7sd  design  is  inferior  to  the  26ssdd  design,  especially 
for  the  reconstruction  of  absorption  properties.  However,  the  accuracy  of  reconstructing 
scattering  properties  by  the  7sd  design  is  close  to  that  of  26ssdd.  Considering  the  challenges  in 
trans-rectal  NIR  probing  for  coupling  with  TRUS,  it  is  fair  to  develop  and  test  the 
instrumentation  with  the  7sd  design. 
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Fig.  14.  Comparison  of  two  geometries  for  a  target  varying  in  size  in  the  middle-sagittal  plane: 
(a)  illustration  of  the  target  size  change  [Unit:  mm'1];  (b)  comparison  of  absorption  coefficient 
reconstruction;  (c)  comparison  of  reduced  scattering  coefficient  reconstruction 


4.3  Capability  of  recovering  two  targets  by  the  7sd  design  in  sagittal  plane 

Capability  of  differentiating  two  targets  is  of  particular  relevance  to  prostate  cancer  imaging 
owing  to  the  existence  of  secondary  or  multifocal  tumors  [44-46].  The  multiple  lesions  may 
fall  into  the  same  TRUS  field-of-view  (FOV),  or  one  falls  outside  the  TRUS  FOV.  For  the 
former  case,  the  US  prior  could  be  used  to  guide  NIR  reconstruction  of  both  targets.  For  the 
latter  case,  as  NIR  actually  performs  3-D  imaging  it  may  interrogate  the  out-of-plane  target 
and  help  redirecting  the  US  to  that  target.  In  this  section,  however,  we  investigate  only  the 
former  case  of  having  two  targets  in  the  same  sagittal  plane.  This  requires  implementing  the 
multi-target  location  information  in  the  last  step  of  the  hierarchical  spatial  prior  routine.  The 
following  simulations  are  conducted  with  the  7sd  probe  design  only. 

4.3.1  Reconstruction  of  two  targets  located  at  the  same  depth  in  sagittal  plane 

Figure  15  shows  two  lOmm-diameter  regions  being  added  to  the  prostate,  at  coordinates  (25, 
40,  26)  and  (55,  40,  26),  respectively.  In  Fig.  15(a)  only  one  region  has  optical  contrast,  and  in 
Fig.  15(b)  both  regions  have  optical  contrasts.  In  both  cases  the  optical  contrast  can  be 

reconstructed  with  good  accuracy,  as  is  shown  in  Table  4.  The  fia  of  the  target  with  true 

optical  contrast  is  reconstructed  within  ±20%  of  the  set  value  and  jUs  can  be  reconstructed 

within  ±23%  of  the  set  values.  The  target  with  no  optical  contrast  is  reconstructed  with  some 
artifacts,  nevertheless,  the  target  with  optical  contrast  can  be  easily  differentiated  from  the  one 
without. 


(a) 


(b) 


Fig.  15.  Two  suspicious  regions  at  the  same  depth  [Units:  mm  and  mm'1] 
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Table  4.  Comparison  of  reconstructed  optical  properties  (mm1)  in  Fig.  15 


M?(  mm'1) 

Fig. 

Regions 

Peri-prostate 

Rectum 

Prostate 

Target  1 

Target  2 

(a) 

Set  value 

0.002 

0.01 

0.06 

0.06 

0.02 

Reconstructed 

0.002 

0.0101 

0.0601 

0.0778 

0.0208 

(b) 

Set  value 

0.002 

0.01 

0.06 

0.02 

0.02 

Reconstructed 

0.002 

0.01 

0.0597 

0.0207 

0.024 

Vs  (mm'1) 

Fig. 

Regions 

Peri-prostate 

Rectum 

Prostate 

Target  1 

Target  2 

(a) 

Set  value 

0.8 

1.0 

1.27 

1.27 

1.6 

Reconstructed 

0.7995 

0.9935 

1.261 

1.2187 

1.3216 

(b) 

Set  value 

0.8 

1.0 

1.27 

1.6 

1.6 

Reconstructed 

0.8007 

0.9953 

1.2343 

1.2302 

1.2837 

4.3.2  Reconstruction  of  two  targets  located  at  different  depth  in  sagittal  plane 


Two  targets  of  10mm  diameter  are  added  in  the  prostate  at  coordinates  of  (25,  40,  28)  and  (55, 
40,  24),  respectively.  Figure  16  shows  the  reconstructed  images  for  the  case  of  both  targets 
having  negative  contrast  and  the  reconstructed  values  are  listed  in  Table  5. 


(b) 


(c) 


Fig.  16.  Two  targets  at  different  depths:  negative  contrast  cases  [Unit:  mm  and  mm'1] 
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Table  5.  Comparison  of  reconstructed  optical  properties  (mm1)  in  Fig.  16 


(mm1) 

Fig. 

Regions 

Peri-prostate 

Rectum 

Prostate 

Region  1 

Region  2 

(a) 

Set  value 

0.002 

0.01 

0.06 

0.06 

0.02 

Reconstructed 

0.002 

0.01 

0.0594 

0.0674 

0.0252 

(b) 

Set  value 

0.002 

0.01 

0.06 

0.02 

0.06 

Reconstructed 

0.002 

0.0101 

0.0597 

0.0592 

0.0607 

(c) 

Set  value 

0.002 

0.01 

0.06 

0.02 

0.02 

Reconstructed 

0.002 

0.01 

0.0596 

0.0402 

0.0276 

Vs  (mm'1) 

Fig. 

Regions 

Peri-prostate 

Rectum 

Prostate 

Region  1 

Region  2 

(a) 

Set  value 

0.8 

1.0 

1.27 

1.27 

1.6 

Reconstructed 

0.8009 

0.9977 

1.2538 

1.2167 

1.3286 

(b) 

Set  value 

0.8 

1.0 

1.27 

1.6 

1.27 

Reconstructed 

0.8015 

0.9931 

1.2102 

1.2179 

1.2083 

(c) 

Set  value 

0.8 

1.0 

1.27 

1.6 

1.6 

Reconstructed 

0.8015 

0.9956 

1.2177 

1.2024 

1.2962 

When  the  target  is  at  a  depth  of  24mm,  the  jua  and  ^  can  be  reconstructed  within  ±20% 

and  ±25%  of  the  set  values,  respectively.  However,  a  target  at  28mm  depth  cannot  be 
reconstructed.  This  is  due  to  the  prostate’s  high  absorption  coefficient  of  0.06mm'1.  When  the 
prostate  absoption  coefficient  is  reduced  to  0.006mm'1  which  will  provide  positive  optical 
contrast  in  the  two  target  regions,  both  tartgets  can  be  recovered  as  shown  in  Fig.  17  and 
Table  6.  The  jua  of  the  target  can  be  reconstructed  within  ±5%  of  the  set  value,  while  the  ^ 
is  still  reconstructed  within  ±23%  of  the  expected  values. 


Fig.  17.  Two  targets  at  different  depths:  positive  contrast  cases  [Unit:  mm  and  mm'1] 
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Table  6.  Comparison  of  reconstructed  optical  properties  (mm'1)  in  Fig.  17 


mm'1) 

Fig. 

Regions 

Peri-prostate 

Rectum 

Prostate 

Region  1 

Region  2 

(a) 

Set  value 

0.002 

0.01 

0.006 

0.006 

0.02 

Reconstructed 

0.002 

0.0101 

0.006 

0.0083 

0.0191 

(b) 

Set  value 

0.002 

0.01 

0.006 

0.02 

0.006 

Reconstructed 

0.002 

0.0101 

0.006 

0.0116 

0.0058 

(c) 

Set  value 

0.002 

0.01 

0.006 

0.02 

0.02 

Reconstructed 

0.002 

0.0099 

0.006 

0.0208 

0.0199 

Vs  (min'1) 

Fig. 

Regions 

Peri-prostate 

Rectum 

Prostate 

Region  1 

Region  2 

(a) 

Set  value 

0.8 

1.0 

1.27 

1.27 

1.6 

Reconstructed 

0.7993 

0.9934 

1.2637 

1.2689 

1.6722 

(b) 

Set  value 

0.8 

1.0 

1.27 

1.6 

1.27 

Reconstructed 

0.8008 

0.993 

1.259 

1.2402 

1.3586 

(c) 

Set  value 

0.8 

1.0 

1.27 

1.6 

1.6 

Reconstructed 

0.8049 

0.9979 

1.2576 

1.2538 

1.424 

5.  Discussions 

The  scope  of  this  work  is  to  investigate  the  feasibility  of  trans-rectal  NIR  tomography  of  the 
prostate  in  the  context  of  concurrent  imaging  with  sagittal  TRUS  using  combined  endo-rectal 
NIR/US  probe.  Recently  there  have  been  considerable  interests  on  trans-rectal  NIR 
tomography  to  augment  existing  imaging  modalities  which  cannot  be  validated  without  the 
development  of  an  endo-rectal  NIR  applicator.  With  an  endo-rectal  applicator,  trans-rectal 
NIR  tomography  of  the  prostate  can  be  performed  stand-alone.  However,  without  a  position 
correlation  with  a  real-time  anatomic  imaging  modality  such  as  TRUS,  the  images  obtained 
by  trans-rectal  NIR  would  be  difficult  to  interpret.  Combining  the  NIR  applicator  with  TRUS 
is  a  viable  solution  for  accurate  positioning  of  the  NIR  probe  which  would  further  enable 
using  TRUS  anatomy  to  guide  the  image  reconstruction  of  trans-rectal  NIR  tomography.  A 
variety  of  designs  of  NIR  array  for  coupling  with  TRUS  are  possible,  however,  the  NIR  probe 
dimension  and  the  number  of  NIR  channels  on  the  probe  are  quite  limited. 

The  utilization  of  a  hierarchical  spatial  prior  is  under  the  condition  that  the  anatomic 
information  of  the  prostate  tumor  can  be  extracted  explicitly  from  TRUS.  The  prostate 
boundary  can  be  well-delineated  in  TRUS,  and  so  does  a  strongly  hypo-echoic  region 
indicating  a  suspicious  lesion.  This  is  when  the  NIR  imaging  may  help  determine  if  the 
suspicious  lesion  is  malignant  or  benign  based  on  optical  contrasts  [31].  However,  since  as 
many  as  40%  of  the  tumors  may  be  seen  as  iso-echoic  on  TRUS,  the  utility  or  accuracy  of  this 
hierarchical  imaging  approach  is  hindered  when  TRUS  images  do  not  specify  a  suspicious 
region,  or  when  it  is  difficult  to  define  the  spatial  extent  of  a  suspicion  region  in  TRUS.  Under 
these  circumstances,  the  third  step  of  the  proposed  hierarchical  image  reconstruction  routine 
may  be  performed  by  element-based  reconstruction  within  the  prostate  instead  of  region-based 
reconstruction  for  the  prostate.  Such  approach  is  proven  effective  based  on  our  initial 
investigations,  but  the  accuracy  and  robustness  may  be  affected  by  the  depth-dependent 
sensitivity  of  the  endo-rectal  NIR  probe  and  the  relatively  small  number  of  source/detector 
channels  that  may  be  engineered  on  the  NIR  probe.  More  dedicated  investigations  could  be 
conducted  when  the  trans-rectal  NIR/US  approach  is  experimentally  demonstrated.  Prostate 
trans-rectal  optical  imaging  is  a  relatively  new  area,  thus  it  is  imperative  to  focus  the  initial 
approach  of  trans-rectal  NIR/US  on  characterizing  lesions  identifiable  on  TRUS. 
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Characterizing  lesions  marginally  suspicious  to  TRUS  or  iso-echoic  on  TRUS  may  become 
less-challenging  when  the  knowledge  on  trans-rectal  NIR  tomography  of  lesions  most 
suspicious  on  TRUS  is  available. 

The  simulations  studies  presented  here  are  largely  based  on  setting  the  absorption 
properties  of  prostate  at  a  high  level  of  0.06mm"1.  This  is  significantly  larger  than  that  of 
breast  tissue,  and  it  is  this  parameter  that  dominates  the  detection  depth  of  trans-rectal  NIR 
tomography  for  the  NIR  array  being  discussed.  Recent  studies  indicate  that  improved 
measurement  condition,  such  as  suppressing  the  bleeding  interference,  may  lead  to  a  lower 
absorption  value  of  the  prostate  [22].  If  the  prostate  is  measured  at  its  intact  or  native  states  by 
modalities  such  as  trans-rectal  optical  tomography,  an  even  lower  absorption  coefficient  of  the 
prostate  may  be  obtained.  A  lower  absorption  coefficient  will  allow  the  same  NIR  array  to 
detect  deeper  targets.  In  fact,  a  TRUS-coupled  trans-rectal  NIR  tomography  may  not  only  help 
characterize  lesions  suspicious  to  TRUS,  but  also  help  quantify  the  optical  properties  of  intact 
prostate  that  are  unavailable  so  far. 

6.  Conclusions 

The  feasibility  of  imaging  the  prostate  using  a  TRUS-coupled  NIR  applicator  is  investigated 
by  simulation.  A  hierarchical  iteration  algorithm  is  first  developed  in  order  to  incorporate  the 
TRUS  spatial  prior  more  reliably  into  the  trans-rectal  NIR  image  reconstruction.  This 
hierarchical  reconstruction  method  uses  a  cascaded  initial-guess  approach  to  mitigate  the  local 
minimum  problem  common  to  NIR  tomography  reconstruction.  It  is  shown  that  trans-rectal 
optical  tomography  based  on  this  method  is  reliable.  This  hierarchical  reconstruction  method 
is  then  utilized  to  evaluate  a  number  of  designs  of  NIR  applicator  that  may  be  integrated  with 
a  sagittal  TRUS  transducer.  A  configuration  of  the  endo-rectal  NIR  applicator  is  proposed, 
that  contains  single  line  of  optode  on  each  lateral  side  of  the  sagittal  TRUS  transducer,  with 
20mm  lateral  separation  between  the  two  line  arrays  and  10  mm  longitudinal  spacing  among 
the  total  7  channels  on  each  line-array.  The  performance  of  this  simple  NIR  array  design  is 
evaluated  for  the  imaging  of  single  tumor  target  in  prostate  by  comparing  with  a  much  more 
complicated  design  that  is  impractical  for  endo-rectal  application.  The  simple  NIR  array 
design  is  also  evaluated  for  imaging  of  two  targets  in  the  prostate.  Results  suggest  that  trans- 
rectal  imaging  of  the  prostate  is  feasible  by  coupling  this  simple  NIR  array  with  TRUS. 

The  following  Part-II  paper  presents  the  instrumentation  of  a  TRUS-coupled  endo-rectal 
NIR  array  and  demonstrates  trans-rectal  optical  tomography  of  the  prostate  by  the  combined 
endo-rectal  NIR/US  applicator.  The  endo-rectal  NIR  array  has  incorporated  the  design 
suggested  by  Fig.  10(a).  Concurrent  trans-rectal  imaging  is  acquired  in  the  same  sagittal  plane 
by  both  US  and  NIR  optical  tomography.  The  real-time  TRUS  is  used  for  accurate  positioning 
of  the  endo-rectal  NIR  applicator  and  for  guiding  NIR  image  reconstruction  with  the  spatial  a 
priori  information.  Tests  on  phantoms  and  tissues  using  the  combined  trans-rectal  NIR/US 
imager  demonstrate  that  optical  contrast  may  be  recovered  by  endo-rectal  NIR  imaging  only 
but  with  improved  accuracy  when  the  TRUS  spatial  prior  is  incorporated.  Trans-rectal 
imaging  of  a  healthy  canine  prostate  in  situ  administered  with  tissue  contrast  validates  the 
endo-rectal  utility  of  the  NIR/US  probe  as  well  as  the  hierarchical  method  for  TRUS  guided 
trans-rectal  NIR  image  reconstruction. 
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Abstract:  We  demonstrate  trans-rectal  optical  tomography  of  the  prostate 
using  an  endo-rectal  near-infrared  (NIR)  applicator  integrated  with  a  trans- 
rectal  ultrasound  (TRUS)  probe.  The  endo-rectal  NIR  applicator 
incorporated  a  design  presented  in  our  previously  reported  work.  A 
continuous-wave  NIR  optical  tomography  system  is  combined  with  a 
commercial  US  scanner  to  form  the  dual-modality  imager.  Sagittal  trans- 
rectal  imaging  is  performed  concurrently  by  endo-rectal  NIR  and  TRUS. 
The  TRUS  ensures  accurate  positioning  of  the  NIR  applicator  as  well  as 
guides  NIR  image  reconstruction  using  the  spatial  prior  of  the  target.  The 
use  of  a  condom,  which  is  standard  for  TRUS,  is  found  to  have  minimal 
effect  on  trans-rectal  NIR  imaging.  Tests  on  avian  tissues  validates  that  NIR 
imaging  can  recover  the  absorption  contrast  of  a  target,  and  its  accuracy  is 
improved  when  the  TRUS  spatial  prior  is  incorporated.  Trans-rectal  NIR/US 
imaging  of  a  healthy  canine  prostate  in  situ  is  reported. 

©2008  Optical  Society  of  America 

OCIS  codes:  (170.3880)  Medical  and  biological  imaging;  (170.6960)  Tomography;  (170.7230) 
Urology;  (170.1610)  Clinical  applications. 
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1.  Introduction 

This  is  a  continuation  of  previously  reported  work  [1]  whose  objective  is  to  demonstrate  the 
feasibility  of  trans-rectal  ultrasound  (TRUS)  coupled  trans-rectal  near-infrared  (NIR) 
tomography  of  the  prostate. 

Prostate  cancer  screening  is  performed  by  measurement  of  serum  prostate-specific  antigen 
(PSA),  digital  rectal  examination  (DRE),  and  a  combination  of  these  tests  [2].  When  the 
suspicion  of  prostate  cancer  is  raised  by  abnormal  PSA  and/or  DRE,  the  diagnosis  is  made  by 
biopsy  performed  under  US  guidance  (most  often  by  TRUS).  Biopsy  is  also  used  to  confirm 
neoplastic  lesions  and  to  determine  their  clinical  significance  for  treatment  planning  [3].  The 
current  standard  of  prostate  biopsy  is  to  routinely  use  10  to  12  cores  of  tissue  obtained 
throughout  the  prostate  for  the  initial  assessment  [4].  It  should  be  noted  that  the  majority  of 
biopsies  are  found  to  be  negative,  and  in  men  with  persistent  suspicion  of  prostate  cancer  after 
several  negative  biopsies,  more  extensive  protocols  (>12  cores)  up  to  saturation  biopsy  (24 
cores)  represent  a  necessary  diagnostic  procedure  [5].  TRUS-guided  prostate  biopsy  is 
performed  following  a  “systematic  sampling”  strategy  because  less  than  60%  of  neoplastic 
lesions  appear  hypoechoic  on  TRUS  while  most  of  the  remaining  neoplastic  lesions  appear 
isoechoic  [6].  Also,  TRUS  does  not  reliably  differentiate  neoplastic  from  benign  tumors  [7]. 

We  have  suggested  [1]  that  the  accuracy  of  prostate  biopsy  may  be  improved  if  the  TRUS 
imaging  could  be  augmented  with  a  functional  or  “surrogate”  marker  of  a  prostate  tumor. 
Based  on  decades  of  research  on  cancer  imaging  [8-14]  and  prostate  measurements  [15-22], 
near-infrared  (NIR)  tomography,  being  non-ionizing  and  minimally-invasive  similar  to  TRUS, 
has  the  potential  of  providing  such  functional  or  “surrogate”  markers  of  prostate  tumors.  NIR 
optical  tomography,  if  carried  out  trans-rectally,  may  improve  the  specificity  of  TRUS 
imaging  such  that  prostate  biopsies  could  be  directed  to  the  most  suspicious  lesions. 

In  the  first  paper  of  this  two  part  series  [1]  we  discussed  the  challenges  of  trans-rectal 
optical  tomography,  particularly  the  fabrication  of  an  endo-rectal  NIR  applicator  for 
integrating  with  TRUS.  Specifically,  since  a  TRUS-coupled  endo-rectal  NIR  applicator 
requires  deploying  many  optodes  in  a  very  limited  space.  We  suggested  [1]  that  arranging 
optodes  longitudinally  is  feasible  and  the  configuration  would  allow  interrogating  deep 
prostate  tissue  in  a  sagittal  imaging  geometry.  Coupling  endo-rectal  NIR  to  TRUS  obtains 
functional  information  otherwise  unavailable  from  TRUS  alone.  Coupling  endo-rectal  NIR  to 
TRUS  also  renders  accurate  and  real-time  anatomic  information  from  TRUS  with  which  to 
correlate  NIR.  The  structural  information  obtained  from  TRUS  may  further  provide  the 
needed  spatial  prior  [14,  23-26]  to  improve  the  accuracy  of  NIR  image  reconstruction.  In  that 
initial  paper  we  further  suggested  [1]  an  array  configuration  of  the  endo-rectal  NIR  applicator 
for  direct  integration  with  a  TRUS  transducer. 

The  work  presented  here  details  the  development  of  an  endo-rectal  NIR  probe  built  on  our 
previous  work  [1]  and  demonstrates  TRUS-coupled  trans-rectal  NIR  tomography  of  the 
prostate.  It  validates  that  trans-rectal  NIR  tomography  helps  characterize  a  target  identified  by 
TRUS,  and  the  accuracy  of  quantitative  imaging  of  a  target  by  trans-rectal  NIR  tomography  is 
improved  when  the  TRUS  spatial  prior  is  incorporated.  TRUS-coupled  trans-rectal  optical 
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tomography  of  a  canine  prostate  in  situ  is  performed,  indicating  the  utility  of  this  approach  for 
in  vivo  imaging. 

2.  Instrumentation 

2.1  Development  of  a  TRUS-coupled  NIR  applicator  for  trans-rectal  optical  tomography 

The  integrated  sagittal-imaging  trans-rectal  NIR/US  applicator  consists  of  a  custom-built  NIR 
probe  and  a  commercial  bi-plane  TRUS  transducer,  as  shown  in  Fig.  1.  The  bi-plane  TRUS 
probe  is  equipped  with  a  proximal  7.5MHz  sagittal-imaging  transducer  and  a  distal  5MHz 
transverse-imaging  transducer.  The  sagittal  TRUS  transducer  occupies  a  60mmxl0mm 
window.  The  diameters  of  the  sagittal  and  transverse  imaging  sections  of  the  TRUS  probe  are 
18mm  and  20mm,  respectively.  Adapting  to  the  un-even  TRUS  cross-section,  the  NIR 
applicator  is  fabricated  to  a  cap-shape  and  attached  to  the  TRUS  probe  (Fig.  1(a)  -(e)).  The 
NIR  array  substrate  was  machined  from  a  black  polycarbonate  material  to  minimize  the 
surface  reflection.  This  substrate  was  then  connected  to  an  aluminum  bracket  and  securely 
fastened  to  the  TRUS  handle  using  a  bottom  clamp.  The  rectangular  TRUS  handle  (Fig.  1(f) 
-(g))  ensured  aligning  the  NIR  applicator  to  the  TRUS  transducer.  A  slot  of  60mmxl0mm 
was  opened  up  in  the  NIR  applicator  to  expose  the  sagittal  TRUS  transducer. 


Source 


(e)  Unit:  mm 


Fiber  conduit 

Handle 
clamp 

(9) 


(0 

Fig.  1.  The  combined  trans-rectal  NIR/US  probe:  (a)  top- view;  (b)  top- view  dimension;  (c) 
front- view;  (d)  side-view;  (e)  side-view  dimension;  (f)  side  and  top  views  of  the  NIR/US 
alignment  method;  (g)  rear- view  of  the  NIR/US  alignment  method. 


The  geometry  of  the  NIR  array  follows  the  design  we  previously  described  [1].  The  NIR 
probe  consists  of  two  linear-arrays,  one  for  the  source  and  the  other  for  the  detector,  separated 
by  20mm  and  placed  on  each  side  of  the  sagittal  TRUS  transducer.  Each  linear-array  consists 
of  7  channels  spaced  10mm  apart  and  covering  60mm  in  length.  The  60mm  long  NIR  array  is 
aligned  precisely  with  the  60mm  window  of  the  sagittal  TRUS  transducer.  The  sagittal  TRUS 
transducer  actually  has  a  longitudinal  field-of-view  (FOV)  of  50mm.  Therefore  the  NIR 
channels  1  to  6  on  each  linear-array  match  with  the  length  of  the  sagittal  TRUS  transducer, 
and  the  NIR  channel  7  of  each  linear-array  is  displaced  longitudinally  10mm  from  the  sagittal 
TRUS  transducer. 

A  metal-coated  600pm-core  diameter  fiber  (Oxford  Electronics)  is  chosen  for  its  thin 
cladding  and  mechanical  strength.  The  7  fibers  of  each  linear-array  are  packaged  into  one 
groove  of  ~4mmx4mm  in  cross-section  formed  in  the  black  substrate.  Bending  the  fiber  inside 
the  small  groove  for  side-firing  at  the  probe  surface  is  impractical.  Instead  micro-optics 
components  are  used  for  deflecting  the  light  side-ways.  As  shown  in  Fig.  2,  each  source 
channel  includes  2  gradient-index  (GRIN)  lenses  and  1  prism  attached  to  the  fiber  while  each 
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detector  channel  has  1  GRIN  lens  and  1  prism  attached  to  the  fiber.  The  GRIN  lens  (Newport 
Corporation)  has  a  pitch  of  0.25,  a  diameter  of  1mm,  a  length  of  2.61mm,  and  a  numerical 
aperture  of  0.46.  The  prism  is  a  coated  1mm  right  angle  micro-prism  (Tower  Optics).  Each 
fiber  is  polished  and  epoxied  to  a  prism  and  a  GRIN  lens  is  attached  to  the  other  side  of  the 
micro-prism  for  illumination  and  detection  at  the  probe  surface.  Each  source  channel  has  one 
GRIN  lens  attached  to  the  proximal  end  of  the  fiber  for  coupling  the  emission  of  a 
superluminescent  diode  (SLD)  using  a  spread-spectral-encoding  configuration  [27].  It  is 
shown  by  ZEMAX  (ZEMAX  Development  Co.)  simulation  that  for  collimated  incident  beam 
using  a  GRIN  lens  and  a  prism  gives  38%  more  coupling  than  without  the  GRIN  lens  for  a 
micro-prism  of  80%  reflection  at  NIR  band  that  is  typical  for  enhanced  aluminum  coating.  A 
coupling  efficiency  improvement  of  10%~15%  is  observed  experimentally  from  the 
completed  fiber  channels.  The  overall  coupling  efficiencies  of  the  14  home-assembled  fiber 
channels  are  shown  in  Table  1.  The  coupling  efficiency  is  at  about  50%,  which  could  be 
improved  if  the  assembly  can  be  made  more  precise. 


Metal  ranting 
multi-mode  fiber 


2.61 


GRIN  £ 
Lens 


-frr 


Metal  coating 

To  multi-mode  fiber 

spectra  meter _ ^ 


GRIN 

Lens  ™ 


- ^ 


GRIN  Lens  x — lh6(eorc)  Coated  Prism 

(a)  source  channel 


F4 


^ — 0.6  (core)  Coated  Prism 

fh)  detector  channel 


Fig.  2.  Micro-optical  configurations:  (a)  source  channel;  (b)  detector  channel.  Unit:  mm 


Table  1.  Measured  coupling  efficiency  of  the  source/detector  fiber  channels 


Source  channels 

si 

s2 

s3 

s4 

s5 

s6 

s7 

Coupling  efficiency 

46% 

48% 

49% 

49% 

51% 

43% 

49% 

Detector  channels 

dl 

d2 

d3 

d4 

d5 

d6 

d7 

Coupling  efficiency 

54% 

48% 

49% 

51% 

53% 

52% 

50% 

2.2  Development  of  a  continuous-wave  NIR  imager  for  coupling  with  US 

The  combined  trans-rectal  NIR/US  imager  is  schematically  illustrated  in  Fig.  3(a)  and  the 
photograph  is  given  in  Fig.  3(b).  The  US  scanner  is  an  ALOKA  SSD-900V  portable  machine. 
The  US  images  are  transferred  to  the  main  computer  of  the  combined  imager  by  a  PCI  image 
acquisition  card  (National  Instruments  PCI- 1405).  The  NIR  imager  uses  a  custom-designed 
superluminescent  diode  (SLD)  (Superlumdiodes  Inc.)  that  is  pigtailed  to  a  multi-mode  fiber 
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Grating 


(a)  (b) 

Fig.  3.  The  combined  trans-rectal  NIR/US  system:  (a)  System  diagram;  (b)  Photo  of  the  system 

on  a  custom  designed  three-layer  cart. 

and  delivers  lOOmW  of  840nm  NIR  light  with  14.2nm  FWHM  bandwidth.  The  SLD  output 
beam  is  dispersed  by  a  1200  groves/mm  grating  and  collimated  unto  linearly  aligned  7  fibers 
connecting  to  the  source  channels  on  NIR  applicator.  NIR  light  with  slightly  different 
wavelengths  are  coupled  to  the  7  fibers  to  form  a  spread-spectral-encoding  of  the  source 
channels  [27].  The  remitted  lights  collected  by  the  7  detection  channels  are  coupled  to  a 
spectrometer  (Acton  Research).  The  signal  corresponding  to  the  individual  source  channels 
are  discriminated  horizontally  by  the  spectrometer.  The  signals  corresponding  to  the 
individual  detector  channels  are  differentiated  vertically  based  on  the  position  of  the  detection 
channels  on  the  spectrometer  entrance  slit.  A  16-bit  intensified  CCD  camera  (Princeton 
Instruments)  acquires  a  complete  set  of  NIR  imaging  data.  The  exposure  time  for  one  frame  of 
data  is  in  the  range  of  100s  milliseconds,  depending  on  the  medium  being  imaged.  The  NIR 
system  resides  on  a  custom-built  cart  that  also  houses  the  US  scanner. 

The  CCD  has  a  maximum  dynamic  range  of  48dB,  which  may  not  be  sufficient  to 
accommodate  the  full  dynamic  range  of  the  signals  when  a  medium  is  highly  absorptive,  as 
the  minimum  and  maximum  source-detector  distances  of  the  trans-rectal  NIR  array  are  20mm 
and  63mm,  respectively.  Fortunately,  the  layout  of  the  dual-line  NIR  array  leads  to  lesser 
signal  dynamic  range  for  the  light  illuminated  from  the  sources  at  the  middle  portion  of  the 
array  than  the  ones  at  the  edges  of  the  array.  Reducing  the  total  dynamic  range  of  the  signal  is 
possible  by  proper  use  of  the  Gaussian  spectrum  of  the  SLD  source.  The  actual  light  couphng 
configuration  is  given  in  Fig.  4(a),  where  the  stronger  spectral  components  are  coupled  to  the 
peripheral  NIR  channels  (such  as  1  and  7),  and  the  weaker  spectral  components  are  coupled  to 
the  middle  NIR  channels  (such  as  3-5).  Comparing  with  the  configuration  of  couphng  the 
uniform  light  to  ah  7  source  channels,  the  strategy  in  Fig.  4(a)  offers  a  15dB  reduction  of  the 
overall  signal  dynamic  range  for  a  medium  of  0.0023mm1  absorption  coefficient  and  1.0mm'1 
reduced  scattering  coefficient  (equivalent  to  a  1%  Intralipid  solution). 
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Fig.  4.  Source  coupling  sequence  with  the  Gaussian-spectrum  (a),  and  the  received  signal 
dynamic  range  compared  with  the  case  when  source  channels  have  even  intensities  (b). 


3.  NIR  image  reconstruction  based  on  continuous-wave  data 


The  steady-state  diffusion  equation  for  a  photon  density  U  at  position  r  can  be  stated  as  [28] 

V  •  D(r)VU (r)  -  pa  (r)  •  U{r)  =  -S(r)  (1) 


where  jua  is  the  absorption  coefficient,  D  =  \$(jUc  ,+A)]"1  is  the  diffusion  coefficient  with  jus 

being  the  reduced  scattering  coefficient,  and  S  is  a  source  term.  For  a  collimated  source  and  a 
detector  at  a  semi-infinite  boundary,  the  diffuse  reflectance  may  be  described  by  [29] 


U(p)  = 


(2) 


where  p  is  the  source-detector  distance,  and  zb  is  a  length  term  determined  by  the  refractive 
index  mismatch  on  the  boundary  [30,  31]. 

The  raw  data  acquired  by  CCD  is  show  in  Fig.  5(a).  If  the  SLD  spectral  components  were 
coupled  orderly  to  the  source  channels  from  1  to  7,  the  higher  intensity  signals  would  have 
been  located  along  a  diagonal  line  in  the  CCD  acquired  image  [27,  32].  The  modified  source 
coupling  configuration  described  in  Fig.  4  leads  to  the  pattern  of  diagonal-shifted  high- 
intensity  signals  as  shown  in  Fig.  5(a). 

The  signal  non-uniformity  among  all  source-detector  pairs  is  calibrated  using  the  1% 
bulk  Intralipid  solution  based  on  the  linear  relationship  between  ln[/?2f/(/?)J  and  p  derived 

from  Eq.  (2),  as  shown  in  Fig.  5(b).  One  calibrated  data  set  corresponding  to  the  complete  7x7 
source-detector  pairs  is  displayed  in  (c)  in  comparison  to  the  non-calibrated  one. 


Fig.  5.  Data  Calibration:  (a)  raw  CCD  data;  (b)  linear  fitting  of  the  measurement  data  based  on 
the  semi-infinite  model;  (c)  calibrated  data  for  the  signal  corresponding  to  the  49  source- 
detector  pairs. 
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Recent  studies  have  demonstrated  that  absorption  and  reduced  scattering  coefficients  can 
be  reconstructed  quantitatively  from  steady-state  measurements,  by  updating  the  D  and  fla 

distributions  to  minimize  a  weighted  sum  of  the  squared  difference  between  computed  and 
measured  data  [28].  The  method  is  equivalent  to  performing  the  same  minimization  by 
mapping  the  DC  signal  to  frequency-domain  [33],  which  we  utilize  in  the  reconstruction. 

The  frequency-domain  diffusion  equation  [34] 

V  •  D(r  )V£/(r,  co)  -  [ \ia  (r )  +  —  j  -U(r,  CO)  =  -S(F,  co)  (3) 

where  CO  is  the  source  modulation  frequency,  is  used  to  map  the  Equ.  (2)  using  finite-element 
discretization  by  setting  CO  to  a  small  value  of  0.1  Hz.  The  optical  properties  are  updated  by 
using 

xk+i  =xk  +a-[JT  (xk)J(xk)  +  Myx  JT  (xk)Av(xk)  (4) 

where  *  is  the  parameters  to  be  optimized,  J  is  the  Jacobian  derived  from  Equ.  (3),  Av  is  the 
forward  projection  error,  A  is  a  penalty  or  regularization  term.  A  small  damping  factor  OC  in 
the  range  of  (0,  1)  is  introduced  to  stabilize  the  convergence  [1]. 

Both  absorption  and  reduced  scattering  coefficients  can  be  reconstructed  from  the  DC 
measurement,  but  it  is  noted  that  it  appears  to  be  generally  more  difficult  for  the  scattering 
image  to  be  reconstructed  than  the  absorption  image  [28].  In  this  work  we  focus  on  recovering 
targets  having  absorption  contrast  to  demonstrate  the  feasibility  of  TRUS-coupled  endo-rectal 
NIR  optical  tomography. 

The  geometry  used  for  image  reconstruction  is  illustrated  in  Fig.  6.  The  NIR  image 
reconstruction  uses  a  3 -dimensional  mesh  representing  80x40x60mm3.  The  TRUS  sagittal 
imaging  is  performed  at  the  mid-plane  of  the  NIR  mesh.  The  TRUS  image  is  used  to  develop 
a  mesh  with  target  spatial  information.  The  NIR  image  is  reconstructed  in  3 -dimension,  and 
displayed  at  the  mid-sagittal  plane  to  correlate  with  TRUS  image.  The  image  reconstruction 
typically  takes  ~10  minutes  on  a  3.0GHz  Pentium(R)  4  PC  for  10  iterations. 


Fig.  6.  NIR  Imaging  geometry:  (a)  NIR  arrays  are  parallel  and  symmetric  to  the  TRUS  sagittal 
plane;  (b)  the  positions  of  NIR  channels  with  respect  to  the  TRUS  image;  (c)  mesh  for  NIR 
image  reconstruction  generated  by  use  of  TRUS  image;  (d)  absorption  image  reconstructed 
from  a  simulation  data  for  an  absorbing  target. 


4.  Performance  evaluation  for  the  TRUS-coupled  trans-rectal  NIR  tomography 

4. 1  Imaging  single  target 

The  performance  of  phantom  imaging  by  the  endo-rectal  NIR  applicator  was  evaluated 
without  the  use  of  TRUS  prior.  A  1%  bulk  Intralipid  solution  was  the  background  medium.  A 
cylinder-shape  solid  phantom  15mm  in  diameter  and  25mm  long  (Fig.  7(a))  was  the  target. 
The  solid  phantom  was  fabricated  from  a  bulk  material  that  was  provided  and  calibrated  by 
the  NIR  imaging  laboratory  of  Dartmouth  College.  The  solid  phantom  has  an  absorption 
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coefficient  of  0.0056mm"1  and  a  reduced  scattering  coefficient  of  1.03mm"1,  in  comparison  to 
0.0023mm"1  and  1.0mm"1  of  the  background.  The  testing  setup  is  depicted  in  Fig.  7(b). 


(a)  (b) 

Fig.  7.  Solid  tissue  phantom  (a)  and  experimental  setup  (b) 


Figure  8  lists  the  images  reconstructed  without  the  TRUS  prior  by  using  a  mesh  of 
homogenous  density  throughout  the  entire  NIR  imaging  volume.  When  this  mesh  is  used  for 
reconstruction  based  on  only  the  NIR  information,  the  optical  properties  are  certainly  updated 
element-by-element  in  order  to  recover  the  heterogeneity  being  imaged.  The  solid  phantom  is 
placed  at  the  mid-sagittal  plane.  The  NIR  images  are  displayed  with  a  FOV  of  80x30mm2. 
The  left  most  optode  is  located  at  10mm  right  to  the  left  edge  of  the  image,  and  the  right-most 
optode  is  10mm  left  to  the  right  edge  of  the  image.  The  (a),  (b),  and  (c)  correspond  to  the 
target  depth  of  17.5mm  at  longitudinal  locations  of  20,  40  and  60mm  (counted  from  left  edge), 
respectively.  The  (d),  (e)  and  (f)  correspond  to  the  phantom  depth  of  12.5mm  at  longitudinal 
locations  of  20,  40  and  60mm,  respectively.  The  target  is  identified  clearly  against 
satisfactorily  recovered  background;  however,  the  absorption  contrast  of  the  target  is 
significantly  underestimated. 

The  same  data  sets  were  reconstructed  using  the  target  spatial  information  obtained  from 
TRUS.  The  TRUS  image,  similar  to  the  one  in  Fig.  6(b),  has  artifacts  around  and  shows  only 
the  lower  half  of  the  cylinder  owing  to  the  shadow  effect.  This  image/artifact  pattern  is 
specific  to  the  solid  cylinder  target  that  reflects  much  of  the  US  signal  on  its  surface.  This  type 
of  artifact  may  not  be  representative  for  tissue  imaging  in  situ.  The  artifact  is  thereby  ignored 
when  incorporating  the  spatial  information  of  the  target,  by  generating  a  mesh  having  a 
homogenous  background  region  and  a  circular  target  region  as  shown  in  Fig.  6(c).  As  there 
are  only  two  regions  to  reconstruct,  the  hierarchical  iteration  routine  introduced  in  [1] 
involves  only  2-steps.  The  results  are  given  in  Table  2  and  Figure  9.  The  absorption 
coefficient  of  the  background  medium  is  reconstructed  at  0.002mm"1  for  all  images.  The 
absorption  coefficient  of  the  target  is  recovered  to  within  14%  of  the  true  value. 
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Fig.  8.  N1R  imaging  of  a  solid  tissue  phantom  having  an  absorption  coefficient  of  0.0056mm'1. 
The  reconstructed  background  absorption  coefficients  are  0.0019mm'1  in  (a),  (b),  (d),  (e),  and 
0.0020  mm'1  in  (c),  (f). 


Table  2.  NIR  image  reconstruction  guided  by  TRUS-  prior 


Center  depth 

12.5mm 

Longitudinal  location  (x) 

20mm 

40mm 

60mm 

True  Value  (mm1) 

0.0056 

0.0056 

0.0056 

Reconstructed  value  (mm1) 

0.0067 

0.0064 

0.0063 

Center  depth 

17.5mm 

Longitudinal  location  (x) 

20mm 

20mm 

20mm 

True  Value  (mm1) 

0.0056 

0.0056 

0.0056 

Reconstructed  value  (mm1) 

0.0064 

0.0063 

0.0061 

ROm iti  8 Mm  in  ROmm 


Fig.  9.  TRUS  guided  NIR  image  reconstruction  for  solid  phantom  in  homogenous  medium 

4.2  Imaging  two  targets 

The  capability  of  recovering  more  than  one  target  by  the  endo-rectal  NIR  probe  is  also 
examined.  Two  cylinder- shaped  solid  tissue  phantoms,  including  the  one  shown  in  Fig.  7(a), 
are  used  as  the  target.  Both  targets  are  15mm  in  diameter  and  25mm  in  length.  The  newly 
added  solid  phantom  has  an  absorption  coefficient  of  0.0064mm"1  and  a  reduced  scattering 
coefficient  of  0.997  mm"1  in  comparison  to  the  other  one  of  0.0056mm"1  and  1.03mm"1.  Figure 
10(a)  and  (b)  are  the  images  for  placing  the  two  targets  at  depths  of  17.5mm  and  22.5mm, 
respectively,  in  the  mid-sagittal  plane,  with  20mm  longitudinal  spacing.  When  the 
reconstruction  is  performed  without  the  target  spatial  information,  the  two  targets  can  barely 
be  differentiated  at  the  depth  of  17.5mm,  and  not  discriminated  at  the  depth  of  22.5mm.  The 
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absorption  contrasts  of  the  targets  are  also  significantly  underestimated  in  both  cases.  Figure 
10(c)  indicates  the  results  of  using  the  spatial  information  from  the  TRUS  image.  The  two 
targets  can  be  recovered  as  having  different  absorption  contrasts,  at  both  depths.  The 
reconstructed  absorption  coefficients  listed  in  Table  3  are  also  close  to  the  true  values. 


17.5mm  (target  depth) 


22.5mni  (target  depth) 

- , - 


0.0019 


0.0021  0.0023  mm' 


0.0019 


0.0021  0.0023  mm*' 


10  20  30  40  50  60 

Longitudinal  location  (mm) 

(a)  (b)  (c) 

Fig.  10.  Imaging  of  multiple  targets:  (a)  NIR-only  reconstruction  when  the  target  depth  is 
17.5mm,  the  targets  can  be  barely  separated;  (b)  NIR-only  reconstruction  when  the  target  depth 
is  22.5mm,  the  targets  cannot  be  separated;  (c)  region  based  reconstruction  when  the  target 
depth  is  22.5mm. 


Table  3.  Reconstructed  absorption  coefficients  of  the  two  targets  in  Fig.  12 


Target  1 

Target  2 

True  value  (mm*1) 

0.0064 

0.0056 

Reconstructed  value  (mm1) 

0.0067 

0.0058 

4.3  Effects  of  condom  on  NIR  tomography 

The  TRUS  probe  is  always  covered  with  a  regular  latex  condom  when  imaging  the  prostate.  A 
TRUS-coupled  NIR  applicator  has  to  be  applied  with  a  condom  if  used  in  the  clinic.  In  the 
previous  applications  of  NIR  tomography  imaging  there  was  no  need  of  applying  a  latex 
barrier  between  an  NIR  applicator  and  the  tissue  being  interrogated,  therefore  the  effect  of  a 
latex  condom  on  NIR  tomography  was  not  reported. 

The  test  being  carried  out  used  the  setup  and  solid  phantom  shown  in  Fig.  7,  by  applying 
US  gel  to  the  NIR/US  probe,  covering  the  NIR/US  with  a  condom,  then  pressing  the  condom 
to  eliminate  any  gas  between  the  condom  and  the  probe.  US  gel  is  not  applied  on  top  of  the 
condom  when  imaging  the  Intralipid,  but  is  applied  for  all  the  tissue  imaging  tests  presented  in 
Section  5.  The  absorption  coefficient  of  the  target  reconstructed  using  TRUS  information  is 
listed  in  Table  4.  The  overall  distribution  of  the  reconstructed  absorption  coefficient  of  the 
target  is  plotted  in  Fig.  11(a),  and  the  distribution  specific  to  target  position  is  plotted  in  Fig. 
11(b).  These  results,  which  are  consistent  with  those  presented  in  Section  4.1,  demonstrate 
that  the  condom  has  minimum  effect  on  NIR  imaging.  Nevertheless,  an  attention  of  the  light 
power  is  observed,  as  it  takes  10%  longer  exposure  time  for  CCD  to  integrate  the  same 
amount  of  signal  when  the  condom  is  applied. 
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Table  4.  Absorption  coefficient  reconstructed  with  the  use  of  condom  on  the  NIR  probe  (mm1) 


Depth  (mm) 

22.5 

25 

27.5 

30 

32.5 

35 

37.5 

True  value 

0.0056 

0.0056 

0.0056 

0.0056 

0.0056 

0.0056 

0.0056 

X=20mm 

0.0065 

0.0058 

0.0070 

0.0066 

0.0048 

0.0043 

0.014 

X=40mm 

0.0053 

0.0051 

0.0047 

0.0049 

0.0043 

0.0052 

0.0063 

X=60mm 

0.0059 

0.0045 

0.0066 

0.0052 

0.0054 

0.0095 

0.0065 

Longitudinal  location  (mm) 

(b) 

Fig.  11.  Test  of  condom  effect  on  the  NIR  imaging:  (a)  reconstructed  absorption  coefficients 
for  all  cases  listed  in  Table  4,  and  The  red  dashed  line  shows  the  true  value;  (b)  data  points  in 
(a)  specific  to  target  depths  and  longitudinal  locations. 


5.  Results  of  tissue  imaging  by  TRUS-coupled  trans-rectal  NIR  tomography 

5. 1  Imaging  solid  absorbing  object  embedded  in  avian  tissue 

The  condom-covered  NIR/US  probe  was  enclosed  within  thick  layers  of  chicken  breast  tissue 
( jua  =  0.006 mm~\  jlis  =  0.757mm-1  [35],  and  a  black  object  (10mm  diameter  xlOmm  length)  was 

embedded  as  the  target.  The  photographs  in  Fig.  12(a)  and  (b)  show  the  setup  of  tissue  sample 
and  the  absorbing  target  being  embedded.  The  tissue  without  the  absorbing  target  is  also 
imaged  as  a  baseline  test.  The  baseline  images  by  US  and  NIR  are  shown  in  Fig.  12(c),  where 
the  NIR  image  is  reconstructed  using  a  homogenous  mesh  as  there  was  no  indication  of  a 
target  on  US.  Figure  12(d)  to  (f)  correspond  to  target  being  embedded  at  different  longitudinal 
locations  at  slightly  different  depths.  The  embedded  object  is  clearly  visible  in  the  US  images 
as  anechoic  shadows.  The  NIR  imaging  without  prior  information  can  clearly  recover  the 
object,  but  with  inconsistent  and  potentially  much  under-estimated  absorption  coefficients. 
When  the  target  location  and  size  information  are  used  to  guide  the  NIR  image  reconstruction, 
the  target  is  recovered  with  consistent  absorption  coefficients  indicating  a  strongly  absorbing 
object. 

5.2  Internal  imaging  of  avian  tissues 

As  shown  in  Fig.  13(a)  the  empty  abdomen  of  a  whole  chicken  was  filled  with  chicken  breast 
tissue  and  a  piece  of  chicken  liver  was  embedded  within  the  breast  tissue.  The  embedded  liver 
shows  up  as  the  hypo-echoic  region  circled  in  Fig.  13(b).  The  NIR  image  reconstructed 
without  any  target  information  is  given  in  Fig.  13(c),  in  which  the  absorptive  mass  does 
correlate  longitudinally  with  the  liver  tissue.  Fig.  13(d)  shows  NIR  image  reconstructed  by 
treating  the  liver  mass  as  a  sphere.  Fig.  13(e)  shows  the  NIR  image  reconstructed  by 
excluding  the  non-tissue  region  from  the  background  of  Fig.  13(b).  Fig.  13(f)  shows  the  NIR 
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image  reconstructed  by  taking  into  account  the  irregular  shape  of  the  liver  mass.  Among  the  3 
spatial  prior  implementations,  the  one  in  Fig.  13(f)  is  the  most  accurate  wherein  the 
absorption  coefficient  of  the  chicken  liver  mass  is  also  the  closest  to  the  realistic  value  [35]. 


Fig.  12.  Avian  tissue  imaging:  (a)  front  view;  (b)  top  view;  (c)  no  inclusion;  (d)  inclusion  is  at 
20mm;  (e)  35mm;  (f)  58mm,  longitudinally,  respectively. 


5.3  TRUS-coupled  trans -rectal  optical  tomography  of  canine  prostate  in  situ 

TRUS-coupled  trans-rectal  NIR  imaging  of  the  prostate  is  conducted  on  a  canine  cadaver.  The 
prostate  was  exposed  and  approximately  0.33ml  of  homogenized  foal  liver  was  injected 
ventral-dorsally,  paramedian  in  the  left  lobe  of  the  prostate.  The  prostate  was  then  enclosed  in 
thick  layers  of  peri-prostate  tissues.  On  the  TRUS  image  of  Fig.  14(a)  the  injected  liver  tissue 
is  visible  by  the  mass  proximal  to  the  center  of  the  prostate  and  the  vertical  hyper-echoic  strip 
at  the  ventral  side  of  the  prostate.  The  large  hypo-echoic  region  at  the  upper  half  was  due  to 
air.  The  mesh  generated  in  Fig.  14(b)  has  excluded  the  air  based  on  the  US  image.  The  rectum 
wall  was  not  outlined  in  the  mesh  because  of  its  close  proximity  and  near  apposition  to  the  US 
transducer  due  to  the  small  size  of  this  canine  cadaver.  The  rectum  layer  may  be  included 
when  imaging  a  larger  subject.  The  finalized  mesh  given  in  Fig.  14(b)  shows  nested-domains 
where  there  is  a  relatively  small  prostate  and  a  much  smaller  target  in  the  prostate.  The  NIR 
image  in  Fig.  14(c)  is  reconstructed  by  applying  the  3-step  hierarchical  reconstruction  method 
introduced  in  the  previous  paper  [1].  A  highly  absorptive  mass  is  clearly  recovered  out  of  the 
injected  liver  tissue.  The  absorption  coefficient  of  the  foal  liver  tissue  in  Fig.  14(c)  is  lower 
than  that  of  the  avian  liver  tissue  in  Fig.  13(f);  nevertheless  both  values  are  at  the  order  of 
0.1mm"1,  indicating  high  absorption  by  both  tissues.  The  reconstructed  absorption  coefficients 
of  the  prostate  and  the  peripheral  tissue  are  given  in  Table  5.  It  is  very  interesting  to  observe 
that  the  numbers  agree  with  the  values  suggested  by  the  literature  [16-18,  35]. 
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Fig.  13.  Internal  imaging  of  avian  tissue  embedded  with  a  piece  of  liver,  (a)  Whole  chicken 
sample;  (b)  US  image  of  the  embedded  liver  tissue;  (c)  NIR  image  reconstructed  without  US 
prior,  (d)  NIR  image  reconstructed  assuming  a  circular  target;  (e)  NIR  image  reconstructed  by 
adding  the  boundary  profile  of  the  background  tissue;  (f)  NIR  image  reconstructed  by  adding 
the  actual  contour  of  the  target. 


(a)  (b)  (c) 

Fig.  14.  Trans-rectal  NIR  /  US  imaging  of  in  vitro  canine  prostate  with  exogenous  tissue  contrast. 


Table  5  Reconstructed  absorption  value  of  in  vitro  canine  prostate  imaging 


Region 

Background 

Prostate 

Injected  liver 

Reconstructed  value  (mm1) 

0.0062 

0.0309 

0.1301 

6.  Discussion 

In  the  NIR  imager  an  840nm  superliminescent  diode  with  ~14nm  bandwidth  is  used  to 
simultaneously  excite  all  source  channels.  The  bandwidth  would  certainly  introduce 
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wavelength-dependent  attenuation  among  source  channels  for  the  same  tissue  chromophore. 
The  reconstructed  absorption  coefficient  is  therefore  a  value  averaged  over  the  band  coupled 
into  the  source  channels.  Luckily  the  absorption  of  the  tissue  chromorphores  like  hemoglobin 
is  less  wavelength-dependent  in  the  close  vicinity  of  840nm.  If  the  source  coupling  method  is 
performed  at  an  additional  wavelength,  such  as  780nm,  to  quantify  the  oxygen  saturation  [36], 
it  may  be  necessary  to  compensate  the  wavelength-dependent  absorptions  among  source 
channels  in  image  reconstruction. 

It  is  shown  that  the  absorption  coefficient  can  be  quantitatively  reconstructed  by  steady- 
state  trans-rectal  NIR  measurements.  The  accuracy  of  reconstruction  for  a  highly 
heterogeneous  domain  can  be  improved  dramatically  by  use  of  the  TRUS  spatial  prior. 
Implementing  frequency-domain  detection  to  the  NIR  system  will  also  allow  more  accurate 
reconstruction  of  the  absorption  coefficient  owing  to  the  reliable  differentiation  of  it  from  the 
scattering  by  true  phase  information. 

Generating  the  target  mesh  from  the  TRUS  image  has  been  based  on  the  assumption  that 
the  entire  imaging  volume  and  the  target  are  symmetric  and  centered  at  the  TRUS  imaging 
plane.  This  is  perhaps  necessary  if  only  one  TRUS  image  is  utilized  to  form  a  3-dimensional 
mesh  for  NIR  image  reconstruction,  but  is  mostly  inaccurate.  This  is  particularly  problematic 
in  imaging  of  the  prostate  if  the  NIR/US  probe  is  to  be  directed  away  from  the  midline  of  the 
prostate.  The  effect  of  asymmetric-domain  may  be  corrected  if  multiple  US  images  are 
available  to  generate  a  mesh  representing  the  true  3 -dimensional  contents  of  the  imaging 
volume  more  faithfully. 

The  trans-rectal  examination  of  the  prostate  on  the  sagittal  view  is  likely  being  interfered 
with  by  the  nearby  bladder  and  pelvic  bone.  If  the  bladder  is  empty,  the  photon  propagation 
may  be  less  interrupted.  If  the  bladder  is  full,  it  may  become  a  relatively  transparent  domain 
for  photon  propagation.  In  this  case  the  bladder  likely  becomes  a  diffusion-void  domain  [37, 
38].  A  transport-equation  based  model  of  photon  propagation  can  be  implemented  [37,  39]  to 
mitigate  the  interference  of  void-domains.  The  interference  from  the  pelvic  bone  is  similar  in 
principle  to  that  of  the  chest  wall  in  breast  imaging  when  using  a  planner  NIR  array  [40].  In 
this  work  we  applied  the  Robin-type  boundary  condition  under  diffusion  equation  to  the  void- 
tissue  interfaces  under  the  diffusion  model.  This  is  a  rather  straightforward  approach  that  is 
proven  helpful  for  improving  the  reconstruction  accuracy,  but  more  studies  are  needed  to 
quantitatively  evaluate  improvement  by  this  approach  in  comparison  to  the  more  rigorous 
transport-based  model. 

The  probe  sensitivity  curves  in  this  paper  and  those  in  the  previous  paper  [1]  were  all 
derived  under  the  assumption  of  an  absorptive  NIR  probe  surface.  This  boundary  condition 
produces  a  peak  of  depth-sensitivity  close  to,  but  not  at,  the  probe  surface.  This  feature  may 
be  desirable  for  prostate  imaging,  since  the  measurement  will  be  less  sensitive  to  the  existence 
of  a  rectum  layer.  On  the  other  hand,  the  peak  also  indicates  that  the  measurement  is  less 
sensitive  to  a  lesion  deeper  in  the  prostate  than  one  closer  to  the  surface  of  the  prostate. 
Prostate  tumors  are  known  to  occur  more  often  in  the  peripheral  zone  which  stretches  abaxial 
from  the  middle-line  of  the  prostate.  For  deeper  target,  we  have  showed  that  with  TRUS  prior , 
the  accuracy  of  target  recovery  can  be  improved  successfully.  However,  more  dedicated 
simulation  and  experimental  studies  are  necessary  to  evaluate  the  likelihood  of  retrieving 
deeper  target  in  vivo. 

In  all  the  tests  presented  in  this  work,  the  target  of  interest  is  visible  or  at  least  sensitive  to 
US,  therefore  US  spatial  information  of  the  target  is  readily  rendered  to  guide  the  NIR  image 
reconstruction.  If  the  target  is  ambiguous  to  US,  it  may  still  be  recognizable  by  NIR-only 
reconstruction.  One  example  is  given  in  Fig.  15  wherein  a  piece  of  chicken  breast  tissue  was 
dyed  with  diluted  Indian  ink  and  then  placed  in  the  middle  of  the  natural  chicken  breast  tissue. 
The  approximate  boundary  of  the  dyed  tissue  is  marked  with  the  dash-circle  in  the  US  image. 
This  embedded  tissue  is  hardly  distinguishable  in  the  TRUS  image;  therefore  the  NIR  image 
reconstruction  was  performed  without  any  a  priori  information.  The  target  is  identified  by 
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NIR  even  though  the  depth  and  size  were  not  accurate.  This  result  may  show  positive 
indication  on  the  clinical  application  of  trans-rectal  NIR/US  measurement  providing  the  fact 
that  up  to  40%  of  prostate  tumors  are  not  sensitive  to  US  [3,  6,  7].  Besides  using  NIR 
information  to  characterize  US-sensitive  lesions,  if  the  NIR  imaging  is  able  to  discover  a 
lesion  that  is  otherwise  non-suspicious  on  US,  the  overall  sensitivity  and  specificity  of  TRUS- 
based  prostate  imaging  may  be  improved. 


(a)  TRUS  Imaging  (b)  Homogenous  mesh  (c)  MR  Imaging  panim'1) 


Fig.  15.  Imaging  of  a  tissue  mass  that  is  ambiguous  to  TRUS  but  sensitive  to  NIR. 

7.  Conclusions 

We  implemented  the  NIR  array  design  suggested  in  our  previous  work  for  direct  integration 
of  the  NIR  array  with  a  TRUS  transducer.  The  combined  NIR/US  probe  &  system  enabled 
concurrent  acquisition  of  trans-rectal  NIR  tomography  and  TRUS  images  in  the  same  sagittal 
plane.  Although  the  NIR  imager  is  constructed  in  CW  mode,  accurate  quantitative 
reconstruction  of  the  absorption  coefficient  is  feasible  with  the  TRUS  spatial  a  priori 
information.  The  use  of  a  condom  is  found  to  have  minimum  affect  on  NIR  tomography 
measurement,  indicating  that  endo-rectal  NIR  may  be  applied  concurrently  with  TRUS  in  a 
clinic  setting  with  no  alteration  of  the  standard  procedures.  An  absorptive  target  may  also  be 
recovered  by  trans-rectal  NIR  only,  but  the  incorporation  of  TRUS  a  priori  information 
allows  trans-rectal  NIR  tomography  to  recover  an  absorption  target  more  accurately,  as 
demonstrated  by  imaging  of  the  alien  tissue  in  the  canine  prostate.  The  previous  Part-I  and 
this  Part-II  papers  together  have  validated  prostate  imaging  as  feasible  by  TRUS-coupled 
trans-rectal  NIR  tomography.  Work  is  on-going  toward  in  vivo  measurement  of  optical 
properties  in  intact  prostate  as  well  as  development  of  prostate  cancer  models  in  the  canine  to 
validate  that  augmenting  TRUS  morphology  information  with  trans-rectally  acquired  NIR 
contrast  may  improve  the  overall  accuracy  of  TRUS  prostate  imaging. 
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ABSTRACT 


A  dual-band  trans-rectal  optical  tomography  system  is  constructed  based  on  an  endo-rectal  near-infrared/ultrasound 
applicator  that  has  been  developed  previously  in  our  laboratory.  The  endo-rectal  NIR/US  applicator  consists  of  a 
commercial  bi-plane  ultrasound  and  a  NIR  probe  attached  to  the  sagittal  ultrasound  transducer.  The  NIR  probe  consists 
of  7  illumination  &  7  detection  channels  that  are  distributed  in  parallel  to  and  aside  the  sagittal  TRUS  transducer.  The 
emissions  from  a  780nm  and  an  830nm  light  sources  are  combined  and  delivered  sequentially  to  the  7  NIR  source 
channels  of  the  endo-rectal  NIR/US  probe.  The  7  NIR  detection  channels  are  coupled  to  a  spectrometer  for  separation  of 
the  signals  at  two  wavelengths  illuminated  from  single  source  channel.  The  dual-band  signals  from  all  source  channels 
are  acquired  sequentially  by  a  CCD  camera  synchronized  with  the  source  switching.  The  acquisition  of  dual-band  trans- 
rectal  optical  tomography  data  is  accompanied  by  position-correlated  concurrent  trans-rectal  ultrasound  imaging.  The 
reconstruction  of  a  target  at  dual-wavelength  illumination  is  guided  by  a  priori  spatial  information  provided  by  the 
sagittal  trans-rectal  ultrasound.  Liquid  phantoms  with  different  hemoglobin  concentration  and  oxygen  saturation  are  used 
to  test  the  feasibility  of  dual-band  trans-rectal  optical  tomography. 

Keywords:  Near- infrared  tomography,  oxygen  saturation,  trans-rectal  ultrasound 


1.  INTRODUCTION 

The  utility  of  non-invasive  imaging  of  the  tissue  oxygenation  by  multi-band  near-infrared  (NIR)  light  has  been 
recognized  for  decades.  The  distinctly  different  spectra  of  hemoglobin  absorption  of  the  NIR  light  when  being 
oxygenated  and  deoxygenated  have  enabled  finding  lesions  of  angiogenesis  or  altered  oxygenation  indicating  malignant 
changes  [1-2].  Studies  on  prostate  cancer  have  found  a  correlation  between  the  cancer  and  increased  micro-vessel 
density  [3].  Changes  of  the  local  oxygen  saturation  were  also  observed  during  the  prostate  photodynamic  therapy  under 
interstitial  NIR  measurement  [4].  Tissue  oxygenation  status  is  also  an  important  prognostic  indicator  for  androgen- 
deprivation  therapy  of  newly  diagnosed  metastatic  prostate  cancer  [5].  There  is  increasing  evidence  that  the  ability  of 
quantifying  the  hemoglobin  concentration  and  oxygen  saturation  in  prostate  is  important  for  the  diagnosis,  prognosis, 
and  treatment  monitoring. 

Encouraged  by  the  applications  of  multi-band  NIR  light  in  breast  cancer  imaging  [6-7],  researchers  recently 
start  to  develop  NIR  tomography  techniques  for  trans-rectal  imaging  of  the  prostate  [8-10].  In  this  work  we  report  the 
progress  on  the  instrumentation  of  dual-band  trans-rectal  NIR  imaging.  The  dual-band  trans-rectal  NIR  tomography  is  to 
be  combined  with  trans-rectal  ultrasound  (TRUS)  to  take  advantage  of  real-time  positioning  information  and  further  the 
anatomic  information  of  US  to  guide  the  NIR  image  reconstruction. 

2.  METHODS  AND  MATERIALS 

A  total  of  three  endo-rectal  imaging  NIR  applicators  have  been  historically  developed  in  our  laboratory  to  investigate 
dual-band  trans-rectal  NIR  tomography  that  may  be  coupled  with  TRUS.  The  first  one  was  an  axial-imaging  probe  of 
13mm  in  diameter  [11].  The  second  one  was  an  axial-imaging  probe  of  20mm  in  diameter  [8]  which  is  comparable  in 
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size  to  the  trans-rectal  ultrasound  transducer.  The  third  probe  was  for  sagittal  imaging  and  could  couple  directly  to  a 
sagittal  TRUS  transducer  for  concurrent  NIR/US  imaging  [10]. 


2.1  The  13mm-diamter  axial-imaging  endo-rectal  NIR  probe 


Figure  1  illustrates  the  13mm-diamter  axial  imaging  NIR  tomography  probe.  It  consists  of  16  optical  channels  including 
8  for  sources  and  8  for  detectors.  Those  source  and  detector  channels  were  interlaced  evenly  on  the  circle.  The  diameter 
of  each  optical  fiber  is  600pm.  A  lOmm-diameter  cone  lens  was  used  for  deflecting  the  light  circumferentially  for  axial- 
imaging.  The  fibers  and  cone-lens  were  mounted  together  in  a  13mm  diameter  black  sleeve  and  sealed  with  transparent 
optical  epoxy.  A  total  of  1 6  sealed  transparent  apertures  were  made  at  the  axial-plane  of  the  probe  for  the  light-path. 
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Figure  1.  Geometry  of  the  13mm  axial-imaging  NIR  tomography  probe 


Figure  2.  Spread-spectral-encoding  based  dual-wavelength  NIR  imaging  system 


For  this  probe  the  dual- wavelength  NIR  imaging  was  configured  in  spread- spectral-encoding  mode.  Shown  in 
Fig.  2,  two  broadband  (~14nm  bandwidth)  superluminescent  diode  (SLD)  light  sources  (780nm  and  830nm)  pigtailed  to 
single  mode  fibers  were  used  in  the  system.  Each  SLD  delivers  14mW  of  power.  The  SLD  outputs  were  collimated  at 
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different  angles  to  a  1200  groves/mm  grating  and  the  dispersed  beams  of  both  bands  overlap.  The  dual-wavelength 
dispersed  beam  was  collimated  by  a  planar-convex  lens  with  a  focal  length  of  300mm,  and  then  coupled  onto  8  linearly 
aligned  source  fibers.  Each  source  channel  coupled  ~1.5nm  bandwidth  out  of  each  SLD  bands  for  illumination  to  the 
medium,  and  the  dual-band  light  was  collected  by  detection  fibers  coupled  to  a  spectrometer  where  the  dual-band  signals 
were  differentiated  by  a  300  grooves/mm  grating  to  two  groups,  each  containing  a  complete  set  of  measurements  for  all 
source  and  detector  channels  for  each  wavelength  band.  The  complete  set  of  the  dual-band  data,  shown  in  Fig.  3  which 
was  acquired  by  a  16-bit  CCD  camera,  was  calibrated  based  on  photon  diffusion  equation  and  underwent  image- 
reconstruction  using  the  NIRFAST  algorithm  developed  by  Dartmouth  NIR  imaging  group.  The  use  of  photon-dissuion 
equation  for  this  13mm  probe  may  not  be  accurate  owing  to  the  smaller  size  of  the  probe  with  respect  to  the  mean 
scattering  distance;  however,  the  reconstructed  dual-band  images  of  the  hemoglobin  and  oxygenation  were  qualitatively 
correct  as  a  result  of  likely  canceling  out  of  the  inaccuracy  from  both  bands. 
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Figure  3.  CCD  image  and  raw  data  of  13mm-probe 


2.2  The  20mm  axial-imaging  endo-rectal  NIR  probe 

Figure  4  illustrates  the  configuration  of  the  20mm  axial-imaging  endo-rectal  NIR  probe.  Similar  to  the  previous  13mm- 
imaging  probe,  there  are  8  source  and  8  detector  optical  channels  spaced  circularly.  Each  optical  channel  was 
constructed  with  a  1.0mm  core  diameter  fiber  and  aligned  to  a  45°  rod  lens  of  2mm  diameter  to  deflect  the  beam 
transversely  for  side-firing.  A  2mm  drum  lens  was  attached  for  illumination  and  beam  focusing.  The  length  of  the  probe 
is  7”  which  was  comparable  to  a  standard  TRUS  probe. 
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Figure  4.  Diagram  of  20mm  endoscope  NIR  tomography  transverse  imaging  probe 

This  20mm  diameter  probe  was  connected  to  a  time-multiplexing  light  delivery  setup  for  dual-band  NIR 
tomography  as  shown  in  Fig.  5.  The  -lOOmW  outputs  from  two  laser  diodes  (FDs)  at  785nm  and  850nm  were  coupled 
by  a  bifurcated  fiber  to  a  fiber  switch  that  was  home-made  using  a  linear  motorized  translation  stage  (Zaber 
Technologies).  The  switching  of  the  light  upon  8  source  channels  was  synchronized  with  the  CCD  acquisition  of  the 
corresponding  signals  received  by  all  detectors. 
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Figure  5.  Dual  wavelength  NIR  tomography  system  for  20mm-probe 

Each  complete  set  of  the  data  consists  of  8  frames  corresponding  to  each  source  channel  being  turned  on  by  the 
fiber  switch.  Fig.  6  is  the  one  frame  data  when  the  source  3  is  switched  on.  Because  of  the  longer  distance  between  the 
source  and  detector,  the  total  64  data  sets  follow  a  better  linear  fitting  curve  than  the  previous  data  sets  for  the  1 3mm 
probe  in  Fig.  3.  The  use  of  fiber  switching  also  minimized  the  cross-talk  among  the  adjacent  source  channels.  The 
weakest  signal  however  may  submerge  in  the  background  noise  owing  to  the  limited  dynamic  range  of  the  CCD. 
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Figure  6.  CCD  image  and  raw  data  of  20mm  probe 


2.3  Combined  trans-rectal  NIR  and  ultrasound  imaging  probe 

A  sagittal  trans-rectal  NIR  imaging  probe  was  developed  for  coupling  with  the  TRUS  and  to  improve  the  imaging  depth 
as  in  sagittal  view  the  optodes  can  occupy  larger  array  for  deeper  interrogation.  The  integrated  sagittal-imaging  trans- 
rectal  NIR/US  applicator  consists  of  a  custom-built  NIR  probe  and  a  commercial  sagittal  TRUS  transducer  operating  at 
7.5MHz,  as  shown  in  Fig.  7.  The  sagittal  TRUS  transducer  occupies  a  space  of  60mmxl0mm.  The  NIR  applicator  was 
fabricated  from  a  black  polycarbonate  material  and  built  to  a  cap-shape  to  attach  to  the  TRUS  probe.  A  slot  of 
60mmxl0mm  was  opened  up  in  the  NIR  applicator  to  expose  the  sagittal  TRUS  transducer. 
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Figure  7.  Combined  trans-rectal  NIR/US  imaging  probe  geometry  (unit:  mm) 

The  source  array  and  the  detector  array,  separated  20mm  laterally,  each  has  7  channels  spaced  at  10mm.  A 
metal-coated  600pm-core  diameter  fiber  (Oxford  Electronics)  was  chosen  for  each  optical  channel.  As  shown  in  Fig.  7, 
each  source  channel  includes  2  gradient-index  (GRIN)  lenses  and  1  prism  attached  to  the  fiber  while  each  detector 
channel  has  1  GRIN  lens  and  1  prism  attached  to  the  fiber.  The  GRIN  lens  (Newport  Corporation)  has  a  pitch  of  0.25,  a 
diameter  of  1mm,  a  length  of  2.61mm,  and  a  numerical  aperture  of  0.46.  The  prism  is  a  coated  1mm  right  angle  micro¬ 
prism  (Tower  Optics).  Each  fiber  was  polished  and  epoxied  to  a  prism  and  a  GRIN  lens  was  attached  to  the  other  side  of 
the  micro-prism  for  illumination  and  detection  at  the  probe  surface.  Each  source  channel  has  one  GRIN  lens  attached  to 
the  proximal  end  of  the  fiber  for  coupling  the  emission  from  the  light  source. 
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The  combined  NIR/US  system  for  dual-band  trans-rectal  mode  imaging  is  shown  schematically  in  Fig.  8.  The 
US  scanner  was  an  Aloka  SSD-900V  portable  machine.  The  US  images  were  transferred  to  the  main  computer  of  the 
combined  imager  by  a  PCI  image  acquisition  card  (National  Instruments  PCI- 1405).  The  source  switching  is  by  the 
time-multiplexing  setup  described  in  Fig.  5.  This  system  takes  7  frames  to  acquire  a  complete  set  of  data.  One  example 
of  the  single-frame  image  for  the  fiber  switch  being  turned  to  source  channel  4  is  depicted  in  Fig.  9.  The  signals  have 
less  dynamic  range  than  do  in  Fig.  6  because  of  the  configurations  of  the  optodes. 
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Figure  9.  CCD  image  and  raw  data  for  trans-rectal  NIR/US  imaging  probe 


2.4  Setup  for  measurements  of  hemoglobin  concentration  and  oxygen  saturation. 

A  setup  shown  in  Fig.  10  was  built  for  measuring  the  blood  hemoglobin  concentration  and  oxygen  saturation.  Oxygen 
and  Nitrogen  gases  were  introduced  to  the  fresh  sheep  blood  (300mL)  in  a  tank  to  administer  the  level  of  blood  oxygen 
saturation  (St02)  which  was  monitored  by  an  oxi-meter.  A  pump  delivered  the  blood  from  the  tank  to  a  cylindrical 
container  placed  2.5mm  away  from  the  endo-rectal  NIR  probe.  Test  started  from  the  lowest  St02  level  in  the  blood  and 
stopped  when  it  became  saturated. 


Figure  10.  Diagram  of  blood  circulation  system  for  measurements  of  St02 
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3.  PRELIMINARY  RESTULS 


3.1  Dual-band  absorption  measurements  by  the  13mm  axial-imaging  endo-rectal  NIR  probe. 


Calibration  of  the  system  at  both  wavelengths  was  carried  out  using  homogenous  medium  with  different  concentrations 
of  diluted  India  ink  (concentration  0.005%~0.5%).  A  container  was  made  into  cylinder  shape  (10mm  inner  diameter) 
from  a  solid  phantom  (0.0056mm"1  of  absorption  and  1.03mm'1  of  reduced  scattering).  The  container  was  fixed  close  to 
the  probe  (2.5mm  away  from  container  surface  to  probe  surface)  in  a  1%  intralipid  background  (0.0023mm"1  in 
absorption  and  1.0mm"1  in  reduced  scattering).  Another  2ml  1%  intralipid  was  poured  into  the  cylindrical  container.  The 
container  wall  was  thinner  than  1mm.  The  total  volume  of  this  cylindrical  container  was  3ml.  The  change  of  the 
absorption  of  adding  one  drop  (0.02ml)  of  the  diluted  ink  (0.5%)  into  the  intralipid  medium  was  previously  calibrated  by 
use  of  transmission  measurements  and  Beer’s  law.  One  frame  data  was  taken  after  adding  each  drop  of  ink  up  to  10 
drops.  The  reconstructed  absorption  of  the  diluted  ink  in  the  container  was  shown  in  Fig.  11  where  a  linear  relationship 
is  found  comparing  with  the  true  values. 
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Figure  11.  True  absorption  coefficient  and  reconstructed  absorption  coefficient  of  diluted  ink 


3.2  Dual-band  measurements  of  blood  oxygenation  by  the  13mm  axial-imaging  endo-rectal  NIR  probe. 

By  using  the  imaging  system  configured  in  Fig.  10  and  the  calibration  results  in  Fig.  11,  the  reconstructed  absorption 
coefficients  of  the  sheep  blood  when  the  oxygen  saturation  level  changed  were  illustrated  in  Fig.  12.  Due  to  the 
fluctuation  of  oxygen  and  nitrogen  flowing  in  the  blood  and  some  other  aspects  including  temperature  control,  leakage 
of  gas  from  the  blood  tank,  the  blood  could  not  be  fully  oxygenated.  However,  the  tendencies  of  decreasing  and 
increasing  absorption  values  with  respect  to  the  increasing  oxygenation  at  the  two  wavelength  bands  were  correct. 
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Figure  12.  Variation  of  the  hemoglobin  absorption  coefficient  along  with  the  St02  changes 
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The  reconstructed  total  hemoglobin  concentration  of  the  sheep  blood  circulated  to  the  cylinder  container  aside 
the  13mm  probe  is  illustrated  in  Fig.  13.  Constant  level  of  the  total  hemoglobin  concentration  was  correctly  imaged.  The 
reconstructed  oxygen  saturation  of  the  sheep  blood  circulated  to  the  cylinder  container  is  illustrated  in  Fig.  14.  The 
images  show  correct  trend  of  increased  St02  levels. 
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Figure  13.  Reconstructed  map  (left)  and  the  level  (right)  of  the  total  hemoglobin  concentration  vs.  oxygenation  change. 


Figure  14.  Reconstructed  map  (left)  and  the  level  (right)  of  the  oxygen  saturation  vs.  the  true  oxygenation  change 


3.3  Dual-band  absorption  measurements  by  the  20mm  axial — imaging  endo-rectal  NIR  probe. 

In  the  previous  measurement  by  use  of  the  13mm  probe,  the  2.5mm  separation  of  the  target  from  the  probe  surface  was 
too  small  to  be  differentiated  by  the  NIR  tomography.  This  is  not  only  related  to  the  relatively  low  resolution  of  the  NIR 
tomography,  but  also  to  the  non-uniform  sensitivity  of  the  axial  imaging  probe  along  the  radial  depth  [12]  as  well  as  the 
use  of  photon-diffusion  equation  without  accounting  for  the  circular  curvature.  Further  test  on  the  20mm  axial-imaging 
probe  confirmed  that  the  depth  of  a  target  is  severely  compromised  at  this  axial-imaging  geometry  without  proper 
analytic  treatment  of  the  boundary  curvature,  as  indicated  in  Fig.  15.  A  solid  cylinder  shape  phantom  (10mm  in  diameter 
and  16mm  in  length)  was  put  into  a  homogeneous  background  of  1%  solution  intralipid.  The  phantom  was  put  10mm 
away  (center)  from  the  20mm  probe  surface.  The  solid  phantom  has  the  optical  properties  of  fia= 0.0064mm"1,  jj!s 

=0.9 1mm"1.  The  high  absorption  regions  reconstructed  at  two  wavelengths  are  consistent,  but  at  a  position  close  to  the 
probe  surface.  Improving  the  depth  localization  of  the  axial-imaging  NIR  tomography  may  be  rendered  by  incorporating 
spatial  information  of  the  target.  Yet,  the  curvature  of  the  axial-imaging  geometry  may  limit  the  depth  of  target  being 
detected. 
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Figure  15.  Stand-alone  optical  reconstruction  of  transverse  images 


3.4  Dual-band  absorption  measurements  by  the  sagittal-imaging  concurrent  NIR/US  probe. 


The  dual-band  measurements  of  absorption  based  on  the  sagittal-imaging  NIR  probe  are  shown  in  Fig.  16.  A  solid 
phantom  target  (jUa  =0. 0056mm"1,  jj!s  =1.03mm_1)  was  put  in  the  middle-sagittal  plane  of  the  probe  at  a  depth  of  10mm 

and  15mm,  respectively,  in  a  homogeneous  background  (1%  solution  intralipid).  Both  targets  appear  at  the  same  depth 
of  10mm  by  NIR  only  reconstruction.  This  target  depth  corresponds  to  the  maximum  depth  sensitivity  of  the  probe.  The 
absorption  value  of  the  target  is  also  underestimated.  Figure  17  shows  the  dual-band  reconstruction  after  incorporating 
the  spatial  information  of  the  target  available  from  the  TRUS.  The  target  spatial  information  is  used  to  develop  a  mesh 
representing  a  target  at  the  correct  location  and  the  background.  The  ^r/or-guided  NIR  reconstruction  at  both  bands 
gives  more  accurate  recovery  of  the  absorption  coefficient  of  the  target. 
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Figure  16.  Stand-alone  optical  reconstruction  for  sagittal  images 


Figure  17.  Sagittal  images  of  combined  trans-rectal  NIR/US  probe  in  different  depth 
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4.  SUMMARY 


We  described  the  development  of  three  NIR  imaging  probes  intended  for  dual-band  endo-rectal  imaging  of  the  prostate 
and  ultimately  coupling  with  TRUS  for  ultrasound-guided  NIR  image  reconstruction.  We  have  demonstrated  that 
oxygenation  changes  of  blood  can  be  quantified  by  the  axial-imaging  probe  being  developed.  Work  is  ongoing  to 
validate  the  dual-band  imaging  of  the  oxygenation  in  blood  samples  and  with  the  sagittal-imaging  NIR  probe,  as  well  as 
conduct  n  vivo  testing  with  cancer  models  developed  in  canine  prostate. 
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ABSTRACT 

An  approach  of  hierarchically  implementing  the  spatial  prior  information  in  trans-rectal  optical 
tomography  is  introduced.  Trans-rectal  optical  imaging  of  the  prostate  deals  with  photon  propagation 
through  the  rectum  wall,  the  peri-prostate  tissue  and  the  prostate.  Reconstructing  a  lesion  in  the  prostate 
is  challenging  due  to  the  structural  complexity  as  well  as  the  optical  heterogeneity.  Incorporating  spatial 
“hard”  a  priori  information  available  from  complementary  imaging  modalities  such  as  trans-rectal 
ultrasound  could  in  principle  improve  the  accuracy  of  trans-rectal  optical  tomography  reconstruction. 
However,  the  reconstruction  is  potentially  subject  to  the  local-minimum  sensitivity  if  the  values  of  all 
regional  optical  properties  are  to  be  initialized  simultaneously.  We  propose  a  hierarchical  spatial  prior 
approach  for  trans-rectal  optical  tomography  reconstmction.  Instead  of  assigning  the  initial  values  to  all 
sub-regions  at  once,  a  region  is  initially  assumed  homogenous,  and  the  reconstructed  optical  properties 
are  used  as  the  initial  guess  for  the  region  as  a  background  when  a  sub-region  is  included  in  the  next  step. 
This  approach  translates  to  a  3 -step  iteration  routine  whereby  the  first  step  reconstmcts  the  entire  imaging 
volume  as  a  single  region,  the  second  step  uses  these  results  as  the  initial  guess  of  peri-prostate  tissue  to 
reconstmct  the  prostate  and  the  rectum  wall,  and  the  third  step  assigns  the  updated  results  as  the  initial 
values  of  3  existing  regions  to  reconstmct  a  lesion  inside  the  prostate.  This  approach,  validated  by 
simulation  and  applied  to  experimental  measurements,  is  more  reliable  in  global  convergence,  robust  in 
imaging  of  single  or  multiple  targets,  and  accurate  for  the  recovering  of  optical  properties. 

Keywords:  Prostate,  optical  tomography,  trans-rectal,  image  reconstmction,  near-infrared. 

1.  INTRODUCTION 

Among  the  men  in  USA,  prostate  cancer  is  the  most  common  cancer  that  is  only  second  to  skin  cancer[1].  Prostate 
cancer  screening  is  performed  by  measurement  of  semm  prostate-specific  antigen  (PSA),  digital  rectal  examination 
(DRE),  and  a  combination  of  these  tests [2].  When  the  suspicion  of  prostate  cancer  is  raised  by  abnormal  PSA  and/or 
DRE,  the  diagnosis  is  made  by  biopsy,  most  often  performed  under  trans-rectal  ultrasound  (TRUS)  guidance,  to 
confirm  the  neoplastic  lesions  and  to  determine  their  clinical  significance  for  treatment  planning[3,4].  It  is  known  that 
most  times  prostate  cancer  presents  as  being  multi-focaf5,  6\  thereby  tissue  sampling  at  multiple  sites  within  the 
prostate  gland  is  necessary. 

The  most  common  appearance  for  prostate  cancer  on  TRUS  is  a  hypoechoic  lesion  in  the  peripheral  zone. 
With  PSA  based  screening  enabling  earlier  cancer  detection,  fewer  overt  abnormal  sonographic  findings  are  being 
detected  at  the  time  of  TRUS  guided  biopsy[7].  The  sonographic  finding  of  the  classic  hypoechoic  peripheral  zone 
lesion  has  a  sensitivity  of  cancer  detection  of  85.5%,  specificity  of  28.4%,  positive  predictive  value  of  29%,  negative 
predictive  value  of  85.2%  and  overall  accuracy  of  43%.  The  prevalence  of  isoechoic  or  nearly  invisible  prostate 
cancers  on  TRUS  ranges  from  25  to  42%.  To  date,  no  biologic  differences  have  been  noted  between  isoechoic  and 
hypoechoic  prostate  cancers[8,9].  The  notably  high  rates  of  iso-echoic  cancers  on  TRUS  imaging  leads  to  increasing 
the  number  of  biopsy  cores  to  improve  the  probability  of  detecting  the  prostate  cancer.  However,  the  current  prostate 
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biopsy  strategy  of  “systematic  sampling”  results  in  large  number  of  negative  diagnosis  performed  for  initial  biopsy 
evaluations. 

TRUS  is  widely  received  as  a  standard  imaging  modality  in  performing  the  prostate  biopsy  owing  to  its 
real-time  utility,  patient-and-physician  friendly  features  as  well  as  cost-effectiveness.  The  prostate  anatomy  and 
needle  track  are  visualized  accurately  with  TRUS  enabling  accurate  localization  of  the  tissue  sites  being  sampled. 
TRUS  will  maintain  as  the  primary  imaging  guidance  for  prostate  biopsy  in  the  future,  thereby  the  quest  of  reducing 
the  numbers  of  unnecessary  biopsies  demands  the  specificity  of  TRUS  be  improved.  Doppler  ultrasonography, 
contrast-enhanced  gray-scale  ultrasonography,  2nd-harmonic  imaging,  etc.  are  among  the  approaches  being 
investigated  to  improve  the  outcome  of  TRUS  imaging[10].  A  more  effective  approach  of  improving  the  outcome  of 
TRUS  prostate  imaging  is,  perhaps,  to  augment  the  TRUS  with  a  complimentary  modality  being  more  specific  to  the 
tumor  malignancy. 

It  has  been  recognized  during  the  past  decades  that  near-infrared  (NIR)  tomography,  being  non-ionizing 
and  non-or-minimally-invasive  similar  to  TRUS,  has  the  potential  of  providing  a  functional  or  “surrogate”  markers 
of  prostate  tumors.  Near-infrared  measurements  of  attenuation  through  tissue  have  demonstrated  significant  contrast 
gradients  between  blood  and  parenchymal  tissue  that  is  otherwise  difficult  to  obtain[11"16].  The  alteration  of 
vascularity  or  the  hemoglobin  content  in  the  tumor  provides,  often  very  high,  intrinsic  optical  contrast  between  the 
tumor  and  benign  tissues  which  has  been  well  demonstrated  in  breast  cancer  imaging.  When  multi-spectral  detection 
is  engaged,  NIR  imaging  is  also  capable  of  directly  quantifying  the  chromorphore  concentrations  important  for 
determining  the  local  tissue  malignancy[11"16].  In  prostate,  studies  have  shown  vascular  density  gradient  in  malignant 
versus  benign  tissue  specimens[17],  and  different  water  concentrations  in  cancerous  and  benign  tissues  in  v/fro[18] 
Invasive  NIR  measurements  of  prostate  have  been  conducted  for  experimental  prostate  tumors[19]  and  human 
prostate[20, 21].  Surface  measurements  of  implanted  prostate  tumor  have  also  been  reported^2, 23].  All  these  studies 
demonstrate  the  potential  of  using  NIR  to  detect  and  characterize  prostate  cancer.  NIR  diffuse  optical  measurement, 
performed  interstitially,  is  also  becoming  an  important  tool  for  monitoring  of  the  photodynamic  therapy  in  prostate 
[20, 21, 24]  prostate  njr  imaging  via  trans-urethral  probing  had  been  analyzed  and  tested[25].  Recently,  trans-rectal 
prostate  NIR  imaging  has  been  investigated  in  simulation  in  the  context  of  assisting  MRI  for  treatment  decision[26]. 
NIR  optical  tomography,  when  carried  out  trans-rectally  and  performed  concurrently  with  TRUS,  will  improve  the 
specificity  of  TRUS  imaging. 

The  implementation  of  trans-rectal  NIR  tomography  for  imaging  the  prostate,  similar  to  other  NIR 
tomography  applications,  will  be  benefited  by  knowing  the  structural  properties  of  tissue  volume  being  imaged.  The 
tissue  volume  interrogated  by  trans-rectal  NIR  imaging  constitutes  a  nested-domain  including  a  thin  layer  of  rectum 
wall,  a  large  volume  of  peri-prostate  tissue,  a  relatively  absorbing  prostate,  and  the  lesion  within  the  prostate.  These 
nested  imaging  domains  may  be  further  complicated  by  the  pelvic  bone,  the  bladder  if  not  completely  voided  and 
other  peri-prostatic  anatomic  structures,  that  could  interfere  with  the  light  propagation  in  the  prostate.  Previous 
investigations  have  discussed  the  issue  of  recovering  the  shapes  and  optical  properties  of  regions  with  optical 
contrast  inside  non-nested  or  nested  domains [27"31],  where  the  shapes  of  the  regions-of- interest  (ROIs)  were  derived 
from  optical  information  when  no  spatial  prior  is  available  from  other  complementary  imaging  modalities.  These 
methods  have  shown  sufficient  robustness  in  recovering  the  shapes  and  optical  properties  of  the  ROIs,  yet  the 
problem  of  stability  and/or  slow  convergence  was  noticed  in  such  approaches  dealing  with  nested-domains.  In  trans- 
rectal  NIR  tomography  reconstruction  the  spatial  information  from  TRUS  may  be  implemented  by  assigning 
homogenous  optical  properties  within  each  ROIs  of  the  imaging  domain.  However,  the  convergence  and  the 
accuracy  of  reconstruction  for  trans-rectal  NIR  tomography  of  the  prostate  will  likely  be  more  dependent  upon  the 
initial  guess,  given  sufficient  accuracy  of  the  spatial  prior  information,  due  to  the  possible  multiple  combinations  of 
optical  properties  in  the  nested-structures.  The  image  reconstruction  in  trans-rectal  optical  tomography  is  further 
complicated  by  the  non-consensus  regarding  the  optical  contrast  that  the  prostate  tumor  has  with  respect  to  the 
benign  prostate  tissue.  All  these  features,  as  shown  in  this  paper,  will  complicate  and  thereby  reduce  the  accuracy  of 
the  reconstruction  in  trans-rectal  NIR  tomography  even  with  the  structural  information  from  TRUS. 

In  this  work,  we  propose  a  hierarchical  reconstruction  scheme  that  may  be  more  robust  and  accurate  for 
structural  prior  guided  NIR  tomography  reconstruction  of  the  prostate.  This  hierarchical  method  improves  both  the 
stability  of  the  convergence  and  the  convergence  to  the  global  minimum  for  the  iterative  inverse  solver.  The 
improvement  and  robustness  of  the  hierarchical  reconstruction  approach  are  demonstrated  via  simulations  and 
experimental  investigations. 
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2.  SIMULATION  CONFIGURATIONS 


2.1  Geometry  of  TRUS  coupled  trans-rectal  NIR  imaging  of  the  prostate 


We  have  designed  and  developed  a  trans-rectal  NIR  probe  that  is  directly  coupled  to  a  bi-plane  TRUS  probe  (Aloka 
UST-672-5/7.5).  Having  more  NIR  channels  in  the  trans-rectal  probe  would  certainly  render  more  advantageous 
imaging  features,  however,  the  probe  dimension  and  fabrication  challenges  posed  limitation  to  integrate  only  7 
source  and  7  detector  channels  to  the  NIR  applicator.  The  NIR  optodes  are  set  laterally  aside  the  sagittal  TRUS 
transducer  spanning  60mmxl0mm.  On  each  lateral  side  of  the  TRUS  transducer,  the  7  sources  or  detectors  are 
spaced  at  1  Omm  thereby  covering  a  total  length  of  60mm.  The  source  array  and  detector  array  are  laterally  spaced  at 
20mm.  The  geometry  of  the  trans-rectal  NIR  probe  and  TRUS  transducer  with  respect  to  the  prostate  imaging 
domain  is  schematically  shown  in  Fig.  1. 


Rectq 

Ultras 

;ound  Transducer — * 

(a) 


Fig.  1  Illustration  of  the  TRUS  coupled  trans-rectal  NIR  imaging  of  the  prostate:  (a)  sagittal  view,  (b) 
transverse  view.  The  middle  plane  of  the  sagittal  TRUS  is  marked  with  a  red  dash  line.  The  red  circles  denote 
NIR  sources  and  the  blue  crosses  denote  the  NIR  detectors. 


2.2  Model-based  trans-rectal  optical  tomography  reconstruction 

The  prostate  and  peripheral  tissues[18"21]  are  scattering-dominant  that  justifies  describing  the  light  propagation  in 
them  by  use  of  the  diffusion  approximation  to  the  radiative  transport  equation.  We  use  the  frequency- domain  photon 
diffusion  equation[32]: 

-V-£>(r)Vt/(r,  *»)  +  (//„  +—  )U(r  ,co)  =  q0{r,co)  (!) 

c 

where  U(r,a>)  is  the  photon  density  at  position  r  ,  qQ(r,co)  is  the  source  term,  CO  is  the  source  modulation 
frequency,  C  is  the  speed  of  light  in  medium,  jda  is  the  absorption  coefficient,  and  p>  =  [3 ^  )]_1  is  the  diffusion 
coefficient  with  being  the  reduced  or  transport  scattering  coefficient. 

Equ.  (1)  is  solved  under  the  Type-III  boundary  condition[32]  by  use  of  finite- element  method[33]  The 
imaging  volume  to  be  reconstructed  in  trans-rectal  optical  tomography  can  be  divided  approximately  to,  as 
previously  described,  4  domains  or  ROIs:  the  rectum  wall,  the  peri-prostate  tissue,  the  prostate,  and  the  prostate 
tumor.  When  concurrent  TRUS  image  is  available,  the  thickness  of  the  rectum  wall,  the  size/shape  of  the  prostate, 
and  the  approximate  spatial  extent  of  the  prostate  lesion  can  be  defined  to  guide  trans-rectal  optical  tomography 
reconstruction.  We  utilize  the  TRUS  structural  information  by  the  “hard”  a  priori  method[34"36].  The  “hard”  a  priori 
method  treats  each  ROI  as  homogenous,  thereby  in  our  reconstruction  8  parameters  are  to  be  recovered  as  pa  and 

ps  of  the  four  ROIs.  The  Jacobian  values  are  calculated  for  each  ROI  instead  of  for  each  node  in  every  ROI  and  has 
the  form  of: 
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where  I..  (i,j  =  1,2,. ..,7)  and  </>..  (i,j  =  1,2,. ..,7)  are  the  intensity  and  phase  terms  of  U(r,co ),  respectively.  In 

(2),  k=l,  2,  3,  4  denote  “rectum  wall”,  “peri-prostate  tissue”,  “prostate”,  and  “prostate  lesion”,  which  are  the  four 
ROIs,  respectively.  The  optical  properties  are  iteratively  recovered  by  use  of  the  well-known  Levenberg-Marquart 
(LM)  algorithm  wherein  the  ROI-specific  values  of  jua  and  D  are  updated  according  to 

**+1  =  Xk  +  a  ■  [JT (xk  )J (xk )  +  All'  JT  (xk  )Av(xk )  (3) 

where  x  is  the  array  of  parameters  to  be  optimized,  J  is  the  Jacobian  matrix  defined  by  Equ.  (2),  A  V  is  the  forward 
projection  error,  Z  is  a  penalty  or  regularization  term.  In  the  iteration  a  small  damping  factor  CC  in  the  range  of  (0, 
1)  is  used  to  stabilize  the  convergence[37]. 

2.3  TRUS  prior  assisted  finite-element  mesh  for  trans-rectal  NIR  tomography  reconstruction 

Utilizing  TRUS  structural  information  to  guide  trans-rectal  NIR  tomography  reconstruction  could  be  performed  by 
using  the  location  and  size  of  a  suspected  lesion  to  directly  generate  a  mesh  that  contains  two  ROIs  to  represent  the 
lesion  with  a  dense  mesh  and  the  background  tissue  with  a  coarse  mesh[35, 36].  For  trans-rectal  prostate  imaging  a 
mesh  is  needed  to  represent  regions  including  the  rectum  wall,  the  peri-prostate  tissue,  the  prostate,  and  the  tumor 
within  the  prostate. 

The  TRUS  image  (available  in  open  source)  was  imported  into  a  software  3ds-MAX  [Autodesk  Inc]  for 
primary  3-D  model  generation.  The  finalized  3-D  mesh  of  the  prostate  is  converted  to  COMSOL  Multiphysics 
[COMSOL  AB]  compatible  format  with  MeshToSolid  [Syncode  Inc].  The  absorption  and  reduced  scattering 
coefficients  of  rectum,  peri-prostate  tissue,  and  prostate  are  assumed  homogenous  in  each  individual  ROI  and 
assigned  based  on[26].  A  spherical  shape  is  adopted  for  the  prostate  tumor.  Figure  2  illustrates  one  completed  NIR 
FEM-mesh.  The  prostate  in  Fig.  2  has  a  walnut  shape  of  -50x40x30  mm3  in  the  maximum  extensions.  The  rectum 
wall  has  a  thickness  of  4mm  and  a  curvature  radius  of  80  mm.  The  completed  mesh  similar  to  that  shown  in  Fig.  2 
contains  approximately  4000  nodes  and  20000  elements. 

The  source-modulation  frequency  CO  in  Equ.  (1)  is  set  at  100MHz,  and  1%  Gaussian  noise  is  added  to  all 
forward  calculations  to  form  the  measurement  data.  In  all  the  figures  presented  thereafter,  the  unit  of  the  length  is 
millimeter,  X  is  the  longitudinal  coordinate  in  the  range  of  [0,  80],  y  is  the  lateral  coordinate  in  the  range  of  [0, 

80],  and  Z  is  the  depth  coordinate  in  the  range  of  [0,  80]. 


(a)  Perspective  view  (b)Axial  view  (b)Sagittal  view  (c)  Coronal  view 

Fig.  2  FEM  mesh  of  the  NIR  imaging  volume  generated  based  on  TRUS  image 

3.  A  HIERARCHICAL  SPATIAL  PRIOR  APPROACH  FOR 
NIR  TOMOGRAPHY  RECONSTRUCTION 
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3.1  Trans-rectal  NIR  imaging  reconstruction  without  a  priori  information 


The  performance  of  recovering  a  tumor  target  by  trans-rectal  NIR  tomography  without  structural  prior  is  first 
examined.  With  the  accurate  forward  measurement  being  calculated  by  the  TRUS-defmed  geometry  as  shown  in  Fig. 
2,  the  iterative  image  reconstruction  is  conducted  using  mesh  of  homogenous  element  density  throughout  the  entire 
volume.  The  optical  properties  of  a  10mm  diameter  tumor  target  is  assigned  as  pa  =  0.02mm-1  and  p  -\ .6mm-1, 

with  the  parameters  of  other  regions  the  same  as  those  in  Fig.  2  and  listed  in  Table  1  coming  later  in  the  text.  The 
reconstruction  results  are  given  in  Fig.  3  for  the  tumor  target  being  placed  at  left,  middle,  and  right  locations  inside 
the  prostate.  The  (a)  is  for  pa  and  (b)  is  for  p  respectively.  The  1 -dimensional  profile  across  the  tumor  and  parallel 

to  the  probe  surface  indicates  the  quantitative  comparison  between  the  set  and  reconstructed  values.  It  is  observed 
that  the  existence  of  a  heterogeneous  structure  in  the  prostate  region  is  rendered,  but  the  spatial  characterization  of 
the  heterogeneity  is  poor.  Further,  a  tumor  with  negative  absorption  contrast  is  difficult  to  be  recovered  when  the 
spatial  prior  is  not  implemented. 


(a)  (b) 

Fig.  3  NIR-only  reconstruction  for  target  in  varying  longitudinal  locations.  Optical  properties  are  plotted  to 
along  line  y=40mm,  z=26mm  (actual  depth  is  15.1mm  because  of  the  bottom  curvature  of  the  mesh)  for  both 
the  target  and  the  reconstructed  images  (a)  Recovery  of  pa  .  pa  values  are  intentionally  plotted  along  the 

line  going  through  the  center  of  the  reconstructed  blobs  (~  y=40mm,  z=26mm,  actual  depth  is  approximately 
15.1mm  because  of  the  bottom  curvature  of  the  mesh),  which  is  shallower  because  the  non-uniform 
sensitivity  has  dragged  the  reconstructed  target  closer  to  the  NIR  array,  (b)  Recovery  of  p  for  the  target 
locations  identical  to  those  in  (a). 


3.2  Spatial  prior  guided  trans-rectal  NIR  reconstruction - simultaneous  updating  of  all  ROIs 

The  trans-rectal  prostate  imaging  domain  constitutes  nested- structures  including  a  thin  layer  of  rectum  wall  in  the 
vicinity  of  NIR  array,  a  large  and  under-defined  volume  of  peri-prostate  tissue  owing  to  the  one-sided  NIR  array,  a 
relatively  absorbing  prostate,  and  the  lesion  within  the  prostate.  In  our  approach  of  coupling  trans-rectal  NIR  with 
TRUS,  the  shapes  of  the  ROIs  for  trans-rectal  optical  tomography  reconstruction  could  be  defined  directly  from 
TRUS.  The  simplest  approach  of  implementing  the  TRUS  a  priori  information  may  be  assigning  homogenous 
optical  properties  within  each  ROIs  of  the  imaging  domain.  However  the  convergence  and  the  accuracy  of 
reconstruction  will  still  be  dependent  on  the  initial  guess.  This  can  be  attributed  to  the  gradient  based  solver  for 
which  the  local  minimum  feature  would  likely  be  exaggerated  in  prostate  imaging  in  which  the  multiple 
combinations  of  optical  properties  of  the  nested-structures  in  the  imaging  volume  could  give  multiple  combinations 
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of  optical  properties  that  may  fit  with  the  measurements,  in  which  case  the  projection  error  of  the  reconstruction 
routine  could  converge  to  a  very  small  value  but  the  recovered  optical  properties  may  be  far  from  the  desired  ones, 
as  is  shown  in  Fig.  4(c).  The  local-minimum  problem  makes  the  reconstruction  sensitive  to  the  initial  guess  of 
optical  properties.  In  trans-rectal  optical  tomography  finding  the  global  minimum  for  the  iterations  to  converge  may 
become  particularly  intriguing. 

When  TRUS  prior  is  available,  a  conventional  way  of  utilizing  the  TRUS  spatial  information  would  be  the 
“hard”  a  priori  method  in  which  the  optical  properties  of  each  ROIs  are  set  homogenous  and  updated  simultaneously 
at  each  iteration.  However,  we  have  found  that  this  conventional  approach  may  not  lead  to  reliable  convergence  for 
prostate  imaging,  which  is  believed  due  to  the  local-minimum  problem.  One  example  is  given  in  Table  1  where  the 
calculation  is  taken  for  the  NIR  array  in  Fig.  1(c),  the  prostate  model  in  Fig.  2,  and  a  target  of  10mm  in  diameter 
located  at  a  coordinate  setting  of  (40,  40,  26)  that  is  15.1  mm  from  the  rectal  surface.  When  the  four  ROIs  including 
the  rectum,  the  peri-prostate  tissue,  the  prostate,  and  the  prostate  tumor  are  updated  simultaneously  from  the  same 
initial  guess  of  pa  =0.0 1mm"1  and  p  =1.0mm_1,  the  iteration  stops  after  1  round  due  to  the  negative  pa  value 

obtained  for  the  rectum  wall.  This  is  due  largely  to  the  effect  that  the  iteration  has  fallen  into  local  minimums  which 
were  the  minimum  values  locally  but  not  globally. 


0.8  1.0  1.2  1.4  1.60.8  1.0  V2  1.4  1.60.8  10  1.2  14  1.6 

(a)  (b)  (c) 


Fig.  4  Local-minimum  issue  in  reconstruction.  The  forward  calculation  is  based  on  (a),  (b)  Reconstruction 
results  with  the  initial  values  of  all  regions  as  ^^O.lmm'1  and  ^  =2mm'1,  the  iterations  stopped  with 

projection  error  of  123.0100.  (c)  Reconstruction  results  with  the  initial  values  of  all  regions  as  pa  = 0.002mm' 
1  and  p  =0.5mm_1,  the  final  projection  error  of  the  iterations  is  0.8329. 

Table  1  Results  of  simultaneously  updating  the  4  ROIs  from  the  same  initial  guess 


“  ^  (mm'1) 

^  (mm"1) 

Regions 

Surrounding 

Tissue 
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Wall 

Prostate 

Tumor 

Surrounding 

Tissue 

Rectum 

Wall 

Prostate 

Tumor 

Set  value 

0.002 

0.01 

0.06 

0.02 

0.8 

1 

1.27 

1.6 

Simultaneous 

Update 

0.1216 

-0.008 

0.026 

0.0215 

1.1482 

2.3602 

0.6173 

0.7073 

3.3  Spatial-prior  guided  trans-rectal  NIR  reconstruction - hierarchical  updating  of  the  ROIs 

The  local  minimum  problem  may  be  mitigated  by  a  hierarchical  spatial  prior  approach  that  may  allow  more  steady 
and  global  convergence  of  the  iteration.  The  fundamental  idea  of  this  method  is  to  first  reconstruct  the  global  optical 
properties  of  the  entire  volume,  then  to  reconstruct  the  optical  properties  of  prostate  and  rectum  wall,  and  last  to 
reconstruct  the  tumor  lesion  area.  The  2nd  and  3rd  steps  use  the  value  obtained  in  the  previous  step  as  the  initial  guess 
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of  that  specific  ROI.  Therefore  at  each  step,  the  perturbation  by  a  relatively  smaller  region  is  less  influential  and 
convergence  of  iteration  is  better  warranted.  The  detailed  steps  are  shown  in  Fig.  5  and  described  in  below: 

(a)  The  first  iterations  assume  an  entirely  homogenous  imaging  volume.  In  this  round  the  initial  projection 
errors  will  be  large  and  the  converging  process  is  most  likely  to  be  affected  by  the  global  minimum.  Therefore,  a 
single  set  of  optimum  values  of  jua  and  are  determined  with  LM  algorithm  and  will  be  used  as  the  initial  guess  in 

the  second  step. 

(b)  The  second  iterations  consider  three  regions  of  rectum  wall,  peri-prostate  tissue,  and  prostate  within  the 
imaging  volume.  The  iterations  of  optical  properties  of  these  three  ROIs  start  at  the  same  initial  guess  provided  in 
step  (a)  and  converge  at  different  values. 

(c)  The  values  obtained  from  step  (b)  are  used  as  the  initial  guess  for  the  same  three  ROIs  but  with  a  tumor 
added  to  the  prostate.  The  tumor  and  the  prostate  take  the  same  initial  values  resulted  from  the  previous  step.  Now 
each  of  the  four  ROIs  (rectum  wall,  peri-prostate  tissue,  prostate,  and  tumor)  converges  to  different  end  values. 

The  comparison  of  the  optical  properties  reconstructed  at  each  step  is  shown  in  Fig.  5(d).  The  change  of  the 
overall  projection  error  for  the  three  steps  is  plotted  in  Fig.  5(e).  Rapid  and  reliable  convergence  can  be  observed. 


Fig.  5.  The  3-step  hierarchical  reconstruction  method  (a)  Step  1 — one  ROI  for  the  entire  volume;  (b)  step  2 — 
three  ROIs  representing  rectum  wall,  peri-prostate  tissue,  prostate;  (c)  step  3 — four  ROIs  representing  rectum 
wall,  peri-prostate  tissue,  prostate,  tumor;  (d)  The  recovered  optical  properties  along  the  z  direction  at 
x=40mm,  y=50mm.  The  dot  line  denotes  the  step  1,  the  dash  line  the  step  2  and  the  solid  line  the  step  3;  (e) 
changes  of  the  projection  error,  where  the  dash  lines  separate  the  converging  of  the  three  steps  in  (a) — (c). 


4.  Validation  of  the  hierarchical  spatial  prior  method 


4.1  Simulation  geometry 

Recently,  Li  et  al.  reported  simulation  results  for  trans-rectal  optical  tomography  reconstruction  in  the  context  of 
using  MRI  anatomic  informational  Their  work  is  referred  as  “NIR/MRI”  in  the  following  text.  We  have  tested  the 
performance  of  our  hierarchical  spatial  prior  method  using  the  probe  geometry  (shown  in  Fig.  6)  and  the  optical 
properties  presented  in  the  NIR/MRI  work. 
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Fig.  6  The  imaging  geometry  used  in  the  NIR/MRI  paper[26].  Two  source  arrays  are  distributed  laterally  with 
20mm  distance.  Each  array  has  5  sources  with  10mm  gaps,  as  denoted  by  red  circles.  28  detectors  are  evenly 
distributed  within  a  15mm  by  30mm  area,  as  denoted  by  blue  crosses,  (a)  Sagittal  view  of  the  distribution  of 
the  optodes.  (b)  Transversal  view  of  the  distribution  of  the  optodes. 
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Fig.  7  Comparison  of  the  reconstruction  values:  (a)  &  (b)  the  tumor  has  a  negative  absorption  contrast  with 
respect  to  the  prostate,  (c)  &  (d)  the  tumor  has  a  positive  absorption  contrast  with  respect  to  the  prostate.  The 
“3 -step”  refers  to  the  hierarchical  image  reconstruction  method. 


With  the  1%  noise  added  to  the  forward  simulation  data  as  is  indicated  in  the  NIR/MRI  paper,  the 
performance  of  our  hierarchical  spatial  prior  method  with  respect  to  the  NIR/MRI  one  is  given  in  Fig.  7  (a)  &  (b)  for 
a  prostate  lesion  of  negative  absorption  contrast.  It  is  observed  that  our  hierarchical  method  (listed  as  “3 -step”  in  the 
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table),  as  expected,  slightly  outperforms  the  NIR/MRI  method  in  terms  of  the  accuracy  of  recovering  optical 
properties.  We  have  also  tested  recovering  a  target  of  positive  absorption  contrast  in  Fig.  7(c)  &  (d)  using  the 
NIR/MRI  probe  geometry  and  our  3 -step  method.  Since  no  results  were  given  in  the  NIR/MRI  work  for  target  of 
positive  absorption  contrast,  only  the  recovered  properties  based  on  our  method  are  presented  for  comparison  with 
the  set  values.  In  (c)  &  (d)  the  absorption  coefficient  of  prostate  is  set  much  lower  than  that  in  (a)  &  (b)  but  the 
tumor  optical  properties  are  kept  the  same  as  those  in  (a)  &  (b).  If  the  absorption  of  prostate  in  (c)  &  (d)  is  kept  the 
same  as  in  (a)  &  (b),  the  positive  absorption  target  is  hardly  reconstructed  at  the  set  location.  The  choice  of  lower 
prostate  absorption  is  for  testing  our  hierarchical  method,  but  it  does  not  necessarily  disagree  with  the  actual  value  of 
prostate  absorption.  As  a  matter  of  fact,  the  absorption  properties  of  prostate  has  shown  large  variation  in  literatures 
where  the  measurement  are  either  taken  from  in  vitro  tissue  or  may  be  interfered  by  bleeding  of  in  vivo  tissue  under 
invasive  measurement.  The  absorption  coefficient  of  non-exposed  prostate  measured  in  situ  &  in  vivo  is  unavailable 
so  far,  and  if  available,  is  likely  lower  than  the  values  reported  in  literatures  for  non-intact  prostate.  The  results  in 
Fig.  7  demonstrate  the  robustness  of  our  hierarchical  spatial  prior  method  in  reconstructing  prostate  lesion  with 
either  negative  or  positive  absorption  contrast. 

4.2  Experimental  validations  of  the  hierarchical  reconstruction 

Two  sets  of  experimental  data  acquired  with  continuous-wave  NIR  imaging  system  were  employed  to  evaluate  the 
hierarchical  reconstruction  scheme.  The  samples  were  administrated  mostly  with  absorption  contrast;  therefore 
results  of  only  absorption  reconstruction  are  displayed  in  the  following  figures. 

Case  I:  Internal  imaging  of  avian  tissues 

As  shown  in  Fig.  8(a),  the  empty  abdomen  of  a  whole  chicken  was  filled  with  chicken  breast  tissue  and  a  piece  of 
chicken  liver  was  embedded  within  the  breast  tissue.  The  embedded  liver  shows  up  as  the  hypo-echoic  region 
circled  in  Fig.  8(b).  A  two-ROI  FEM  mesh  is  generated  to  simulate  the  experiment  setup,  as  is  shown  in  Fig.  8(d). 

Since  there  are  only  two  ROI  in  this  experimental  case,  when  implementing  the  hierarchical  reconstruction 
scheme,  only  a  two-step  reconstruction  is  performed.  The  results  are  shown  in  Fig.  8(e)  &  (f),  respectively.  It  can  be 
observed  that  the  absorption  values  are  reconstructed  quite  close  to  the  true  values  of  a  liver  tissue  indicated  by 
literature  findings.  Comparing  the  NIR  images  reconstructed  without  structural  prior  in  Fig.  8(c)  and  reconstructed 
with  hierarchical  method  in  Fig.  8(f),  the  utility  of  the  hierarchical  reconstruction  method  is  demonstrated. 

Case  II:  Trans-rectal  optical  tomography  of  canine  prostate  in  situ 

The  second  experiment  is  conducted  on  a  canine  cadaver.  The  prostate  was  exposed  and  approximately  0.33ml  of 
homogenized  foal  liver  was  injected  ventral-dorsally,  paramedian  in  the  left  lobe  of  the  prostate.  The  prostate  was 
then  enclosed  in  thick  layers  of  peri-prostate  tissues.  On  the  TRUS  image  of  Fig.  9(a)  the  injected  liver  tissue  is 
visible  by  the  mass  proximal  to  the  center  of  the  prostate  and  the  vertical  hyper-echoic  strip  at  the  ventral  side  of  the 
prostate.  The  large  hypo-echoic  region  at  the  upper  half  was  due  to  absorption  by  thick  layers  of  muscle,  skin  and 
exposal  to  the  air.  The  mesh  is  according  generated  in  Fig.  9(b).  The  rectum  wall  was  not  outlined  in  the  mesh 
because  of  its  close  proximity  to  the  US  transducer  due  to  the  small  size  of  this  canine  cadaver.  The  finalized  mesh 
in  Fig.  9(b)  shows  nested-domains  including  the  prostate  and  a  target  region  in  the  prostate.  The  NIR  images  in  Fig. 
9(c),  (d),  and  (e)  show  the  absorption  distribution  reconstructed  at  each  step  of  the  hierarchical  reconstruction 
method.  A  highly  absorptive  mass  is  clearly  recovered  corresponding  to  the  injected  liver  tissue.  The  absorption 
coefficient  of  the  foal  liver  tissue  in  Fig.  9(a)  is  lower  than  that  of  the  avian  liver  tissue  in  Fig.  8(b);  nevertheless 
both  values  are  at  the  order  of  0.1mm"1,  indicating  high  absorption  by  both  tissues. 
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Fig.  8  Experiment  I:  Internal  imaging  of  avian  tissues,  (a)  Tissue  phantom,  (b)  Ultrasound  imaging,  (c)  NIR- 
only  reconstruction,  (d)  FEM  mesh,  (e)  Reconstruction  step  1:  background,  (f)  Reconstruction  step  2  with 
the  target  structure 


,  v  0.02  0.06  0.1  0.12  0.02  0.06  0.1  0.120.02  0.06  0.1  0.12 
(a)  (b)  (c)  (d)  (e) 

Fig.  9  Experiment  II:  Imaging  on  canine  cadaver,  (a)  TRUS  imaging  (b)  FEM  mesh,  (c),  (d),  (e)  the  3  steps  of 

hierarchical  reconstruction 

5.  DISCUSSION  AND  SUMMARY 

In  this  paper,  a  structural  prior  guided  hierarchical  NIR  tomography  reconstruction  scheme  aimed  at  improving  the 
reconstruction  accuracy  of  trans-rectal  optical  tomography  of  the  prostate  is  presented.  The  method  is  developed 
under  the  framework  of  utilizing  the  trans-rectal  ultrasound  information  as  the  structural  prior.  Nevertheless,  the 
hierarchical  method  could  be  extended  to  implementing  other  high-resolution  imaging  modalities  such  as  MRI  to 
guide  NIR  tomography  and  to  applying  NIR  tomography  in  other  complicated  tissue  domains.  The  accuracy  of  the 
hierarchical  reconstruction  method,  as  indicated  by  the  simulation  and  experimental  investigations,  relies  on  explicit 
demarcation  of  the  inner  structure  of  the  imaging  volume  to  form  the  region-resolved  FEM  mesh.  In  prostate  a  large 
amount  of  the  cancer  may  be  shown  iso-echoic  on  TRUS.  Thereby  the  utility  or  accuracy  of  this  hierarchical 
imaging  approach  is  hindered  when  TRUS  images  do  not  specify  a  suspicious  region.  Under  such  circumstance,  the 
third  step  of  recovering  the  tumor  lesion  in  the  prostate  may  proceed  by  reconstructing  the  optical  properties  on 
every  element  within  the  prostate  but  the  overall  accuracy  of  the  reconstruction  must  be  compromised. 

In  summary,  an  approach  of  hierarchically  implementing  the  spatial  prior  information  in  trans-rectal  optical 
tomography  is  introduced.  Trans-rectal  optical  imaging  of  the  prostate  is  challenging  due  to  the  structural 
complexity  as  well  as  the  optical  heterogeneity.  Incorporating  spatial  “hard”  a  priori  information  available  from 
complementary  imaging  modalities  such  as  trans-rectal  ultrasound  is  potentially  subject  to  the  local-minimum 
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sensitivity  if  the  values  of  all  regional  optical  properties  are  to  be  initialized  simultaneously.  A  hierarchical  spatial 
prior  approach  for  trans-rectal  optical  tomography  reconstruction  does  not  assign  the  initial  values  to  all  sub-regions 
at  once,  instead  a  region  is  initially  assumed  homogenous,  and  the  reconstructed  optical  properties  are  used  as  the 
initial  guess  for  the  region  as  a  background  when  a  sub-region  is  included  in  the  next  step.  This  approach  translates 
to  a  3 -step  iteration  routine  whereby  the  first  step  reconstructs  the  entire  imaging  volume  as  a  single  region,  the 
second  step  uses  these  results  as  the  initial  guess  of  peri-prostate  tissue  to  reconstruct  the  prostate  and  the  rectum 
wall,  and  the  third  step  assigns  the  updated  results  as  the  initial  values  of  3  existing  regions  to  reconstruct  a  lesion 
inside  the  prostate.  The  robustness  of  this  approach  is  demonstrated  in  both  simulations  and  experiments. 
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In  vivo  trans-rectal  ultrasound  coupled  trans-rectal  near-infrared 
optical  tomography  of  canine  prostate  bearing 
transmissible  venereal  tumor 
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Abstract:  In  vivo  trans-rectal  near-infrared  (NIR)  optical  tomography  is  conducted  on  a  tumor-bearing  canine  prostate 
with  the  assistance  of  trans-rectal  ultrasound  (TRUS).  The  canine  prostate  tumor  model  is  made  possible  by  a  unique 
round  cell  neoplasm  of  dogs,  transmissible  venereal  tumor  (TVT)  that  can  be  transferred  from  dog  to  dog  regardless  of 
histocompatibility.  A  characterized  TVT  cell  line  was  homogenized  and  passed  twice  in  subcutaneous  tissue  of 
NOD/SCID  mice.  Following  the  second  passage,  the  tumor  was  recovered,  homogenized  and  then  inoculated  by 
ultrasound  guidance  into  the  prostate  gland  of  a  healthy  dog.  The  dog  was  then  imaged  with  a  combined  trans-rectal  NIR 
and  TRUS  imager  using  an  integrated  trans-rectal  NIR/US  applicator.  The  image  was  taken  by  NIR  and  US  modalities 
concurrently,  both  in  sagittal  view.  The  trans-rectal  NIR  imager  is  a  continuous-wave  system  that  illuminates  7  source 
channels  sequentially  by  a  fiber  switch  to  deliver  sufficient  light  power  to  the  relatively  more  absorbing  prostate  tissue 
and  samples  7  detection  channels  simultaneously  by  a  gated  intensified  high-resolution  CCD  camera.  This  work  tests  the 
feasibility  of  detecting  prostate  tumor  by  trans-rectal  NIR  optical  tomography  and  the  benefit  of  augmenting  TRUS  with 
trans-rectal  NIR  imaging. 

Keywords:  Prostate;  near-infrared  optical  tomography;  trans-rectal;  ultrasound;  transmissible  venereal  tumor. 


1.  INTRODUCTION 

Prostate  cancer  is  prevalent  in  the  men  in  the  U.S.  and  in  the  industrialized  countries  [1].  Most  of  the  prostate  cancers  are 
diagnosed  at  earlier  stage  now  with  the  use  of  the  sensitive  prostate-specific  antigen  (PSA)  test  [2]  aided  by  the  digital 
rectal  exam  (DRE).  Prostate  biopsy,  guided  by  trans-rectal  ultrasound  (TRUS),  is  performed  when  the  pathological 
evidence  is  needed  to  confirm  or  rule  out  the  onset  of  carcinoma  in  prostate  [3]. 

The  real-time  utility  and  the  excellent  morphological  information  of  TRUS  make  it  ideal  for  visualizing  the 
needle  trajectory  during  the  biopsy.  TRUS  had  been  utilized  to  assess  the  malignant  changes  in  prostate  tissue  [4], 
unfortunately,  the  lack  of  specificity  has  limited  the  value  of  TRUS  in  diagnostic  imaging  of  the  prostate  cancer  and 
transformed  the  role  of  TRUS  to  becoming  primarily  an  imaging  tool  for  placing  the  biopsy  needle.  The  accuracy  of 
TRUS  imaging  may  be  improved  by  augmenting  it  with  a  functional  contrast,  thereby  reducing  unnecessary  biopsies  and 
more  accurately  targeting  the  malignant  areas  for  biopsy.  The  near-infrared  (NIR)  light  is  known  to  reveal  significant 
functional  contrasts  between  a  tumor  and  normal  tissue  originated  from  the  angiogenic  and  hypoxic  changes  of  the 
malignancy  [5-11].  However,  using  NIR  alone  is  less  viable  for  prostate  imaging  due  to  its  low  spatial  resolution  and 
non-real  time  image  reconstruction.  An  appealing  approach  is  to  combine  NIR  and  TRUS  to  take  advantage  of  the 
complementary  features  of  NIR  and  TRUS. 

In  vivo  imaging  of  NIR  tomography  of  the  prostate  to  augment  TRUS  is  confronted  with  difficulties  uncommon 
to  the  previous  application  of  NIR  tomography  in  cancer  imaging.  First  of  all,  the  prostate  lies  deep  inside  the  pelvic 
compartment  wherein  the  imaging  by  trans-rectal  probing  may  be  the  optimal  choice.  In  trans-rectal  probing  the  NIR 
light  will  be  attenuated  by  the  condom  (required  when  using  TRUS)  first,  then  the  rectum  wall,  and  the  peri-rectum 
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tissue  before  reaching  the  prostate.  As  the  light  reaches  the  prostate,  reflection  on  the  prostate  capsule  and  the  potentially 
strong  absorption  by  the  peripheral  vascular  structures  may  limit  the  amount  of  light  interrogating  the  deep  prostate 
tissues.  Secondly  there  is  considerable  challenge  of  fabricating  an  endo-rectal  NIR  applicator  which  should  contain  a 
number  of  source  and  detector  channels  if  spatially-resolved  imaging  of  large  tissue  volumes  of  the  prostate  is  intended. 
The  challenge  becomes  more  pronounced  when  the  NIR  is  to  be  combined  with  TRUS  [12]. 

Recently  we  have  constructed  a  trans-rectal  NIR/US  probe  and  validated  the  utility  of  the  combined  trans-rectal 
NIR/US  imaging  in  a  cadaver  canine  prostate  in  situ  [13].  In  this  paper,  we  demonstrate  the  in  vivo  utility  of  the  trans- 
rectal  NIR/US  imaging  by  use  of  a  canine  prostate  bearing  transmissible  venereal  tumor  (TVT).  TVT  has  an  anechoic 
response  that  makes  it  visible  on  the  TRUS.  The  TVT  is  known  to  have  a  higher  cell  density  that  tends  to  absorb  and 
potentially  scatter  more  NIR  light  than  does  the  normal  prostate  tissue.  In  this  work,  the  TVTs  in  the  prostate  and  peri- 
prostatatic  tissue  were  located  by  TRUS,  then  imaged  by  NIR  tomography  in  situ.  The  TVTs  are  shown  as  strong  hyper- 
absorptive  and  moderately  hyper-scattering  on  NIR  tomography  images  being  reconstructed  without  any  prior. 

2.  METHODS  AND  MATERIALS 

2.1  Sagittal-imaging  trans-rectal  NIR/US  system 

The  details  of  the  integrated  sagittal-imaging  trans-rectal  NIR/US  applicator  and  system  are  described  elsewhere  [13]. 
Briefly,  this  combined  imaging  applicator  consists  of  a  custom-built  NIR  probe  and  a  commercial  bi-plane  TRUS 
transducer.  The  7.5MHz  sagittal-imaging  transducer  of  the  TRUS  probe  is  used  to  couple  with  the  NIR.  The  NIR 
applicator  was  fabricated  to  a  cap-shape  and  attached  to  the  TRUS  probe  on  the  top  with  a  slot  opened  to  expose  the 
sagittal  US  transducer.  The  NIR  array  substrate  was  fabricated  from  a  black  polycarbonate  material  to  satisfy  the 
absorptive  boundary  condition.  The  NIR  probe  consists  of  two  linear  arrays,  one  for  the  source  and  the  other  for  the 
detector,  separated  by  20mm  and  placed  on  each  side  of  the  sagittal  TRUS  transducer  of  60mm  in  length.  Each  linear- 
array  consists  of  7  channels  spaced  10mm  apart  and  covering  a  total  length  of  60mm.  Each  NIR  channel  consists  of  a 
metal-coated  fiber  epoxied  to  a  gradient-index  lens  and  a  prism  to  deflect  the  light  to  and  from  the  probe  surface. 

The  geometric  relation  between  the  imaging  probe  and  the  prostate  is  depicted  in  Fig.  1  where  the  optical 
channels  are  seen  arranged  cranial-caudally  from  1  to  7.  The  NIR  light  illuminates  through  the  condom,  the  rectum,  and 
propagates  into  the  prostate. 


Fig.  1  Illustration  of  the  TRUS  coupled  trans-rectal  NIR  imaging  of  the  prostate.  The  TRUS  middle-sagittal  plane 
is  marked  by  a  red  dash-line.  The  red  circles  denote  NIR  sources  and  the  blue  crosses  denote  the  NIR  detectors. 

The  TRUS  probe  is  connected  to  an  Aloka-3500  scanner.  The  US  images  were  transferred  to  the  main  computer 
of  the  combined  NIR/US  imager  by  a  PCI  image  acquisition  card  (National  Instruments  PCI- 1405).  The  NIR  imager 
uses  a  custom-designed  superluminescent  diode  (SLD)  (Superlumdiodes  Inc.)  that  is  pigtailed  to  a  multi-mode  fiber  and 
delivers  lOOmW  of  840nm  NIR  light  with  14.2nm  FWHM  bandwidth.  The  SLD  output  beam  is  focused  and  delivered 
sequentially  onto  7  source  fibers  of  the  NIR  applicator  by  a  translating  fiber  multiplexer.  The  entire  NIR  bandwidth  is 
coupled  to  each  fiber  at  a  time.  The  remitted  NIR  light  is  collected  by  7  detection  fibers  coupled  to  a  spectrometer 
(Acton  Research).  A  16-bit  intensified  CCD  camera  (Princeton  Instruments)  acquires  the  7  optical  signals  corresponding 
to  one  source  simultaneously,  and  a  complete  set  of  data  for  all  7  source  channels  is  taken  within  5  seconds. 

2.2  Animal  model 

Dogs  have  been  previously  used  for  prostate  tumor  model  studies  because  of  the  similarity  between  canine  and  human 
prostate  glands  [14-16].  For  this  study,  the  prostate  of  a  mature  dog  was  injected  with  canine  TVT  cells  with  the  purpose 
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of  developing  an  in  vivo  prostate  tumor  model  for  NIR/US  imaging.  The  unique  canine  TVT  is  a  round  cell  tumor  of 
dogs  that  mainly  affects  the  external  genitalia  and  can  be  transmitted  from  animal  to  animal  during  copulation, 
regardless  of  histocompatibility.  The  neoplastic  cells  are  generally  thought  to  be  of  histiocytic  origin  [17].  TVT  cells  can 
be  propagated  in  immunocompromised  (SCID)  mice  and  transferred  to  different  tissues  (such  as  the  prostate  gland)  of 
the  dog  to  result  in  a  neoplastic  mass  effect  useful  for  imaging  studies  [16]. 

The  experiments  were  conducted  in  the  Center  for  Veterinary  Health  Sciences  at  Oklahoma  State  University 
under  a  protocol  approved  by  the  university's  Institutional  Animal  Care  and  Use  Committee.  The  protocol  was  also 
approved  and  underwent  an  on-site  inspection  by  the  US  Army  Medical  Research  and  Material  Command.  For  these 
studies,  an  existing  TVT  cell  line  was  obtained  as  cryopreserved  tissue  from  MD  Anderson  Cancer  Center  (Houston, 
TX).  The  tumor  tissue  was  quick-thawed,  homogenized  and  inoculated  into  the  subcutis  of  a  NOD/SCID  mouse.  After 
approximately  12  weeks,  the  tumor  reached  an  appropriate  volume  (4. 0-4. 5cm3)  for  recovery  and  were  processed  and 
inoculated  into  a  second  NOD/SCID  mouse.  Following  growth  in  the  second  mouse  (8  weeks),  the  neoplastic  cells  were 
recovered  and  homogenized  for  inoculation  into  the  canine  prostate  gland. 

A  12  kg  sexually  intact  adult  Beagle  dog  estimated  to  be  approximately  4  years  of  age  was  housed  single  in  a 
inside  run  and  given  free  access  to  water  and  food.  For  TVT  cell  injection,  the  dog  was  anesthetized  using  an 
intravenous  injection  of  propofol  (8mg/kg)  followed  by  intubation  and  halothane/oxygen  inhalation  anesthetic 
maintenance.  The  animal  was  placed  in  left  lateral  recumbency  for  bowel  preparation  and  physical  examination  of  the 
prostate.  The  perineum  was  aseptically  prepared  for  injection  and  under  TRUS  visualization,  a  6-inch  16  gauge 
hypodermic  needle  was  inserted  transperineally  into  the  right  lobe  of  the  prostate  where  3  cc  of  the  prepared  TVT  cells 
were  injected.  During  retraction  of  the  injection  needle,  it  was  considered  unavoidable  that  TVT  cells  would  leak  from 
the  prostate  injection  site  and  be  "seeded"  along  the  needle  insertion  tract. 

The  dog  was  examined  weekly  by  TRUS  using  a  condom-covered  combined  NIR/US  probe.  For  all  TRUS  and 
NIR  follow-on  imaging  sessions,  the  dog  underwent  general  anesthesia  using  a  similar  protocol  used  for  TVT  cell 
injection.  During  the  first  two  weeks  following  TVT  cell  injection  there  was  no  evidence  of  tumor  growth.  The  next 
TRUS  examination  occurred  5  weeks  post- injection  when  hypoechoic  masses  were  observed  both  in  the  prostatic 
parenchyma  and  peri-prostatically  around  the  right  lobe  of  the  prostate.  The  dog  underwent  weekly  monitoring  for  2 
more  weeks,  and  was  then  humanly  euthanized  8  weeks  after  the  initial  TVT  cell  injection.  A  complete  necropsy  was 
performed  with  the  prostate  and  peri-prostatic  structures  submitted  for  histological  examination. 

3.  SYSTEM  PERFORMENCE 

The  combined  NIR/US  probe  &  system  enables  concurrent  acquisition  of  trans-rectal  NIR  tomography  and  TRUS 
images  on  the  same  sagittal  plane.  Our  previous  studies  have  validated  that  incorporating  TRUS  a  priori  information 
allows  trans-rectal  NIR  tomography  to  recover  an  absorption  target  accurately.  We  have  also  demonstrated  that  trans- 
rectal  NIR  imaging  could  recover  an  absorptive  target  with  no  spatial  prior.  In  this  paper,  the  TRUS  is  used  to  locate  a 
suspected  TVT  target,  and  NIR  image  reconstruction  is  performed  without  prior  information  in  order  to  examine  if  the 
TVT  shows  inherent  NIR  contrast.  As  our  trans-rectal  NIR  probe  allows  3 -dimensional  imaging,  the  system  performance 
is  evaluated  for  recovering  targets  located  on  or  off  the  sagittal  plane  when  using  only  the  NIR  information. 

3.1  Reconstruction  of  a  target  on  the  mid-sagittal  plane - Simulation 

The  performance  of  reconstructing  absorptive  targets  in  the  sagital-plane  using  continuous-wave  measurement  has  been 
investigated  previously  [13].  As  the  optical  scattering  properties  of  the  prostate  cancer  and  the  benign  tissue  may  be 
different,  simulation  is  performed  to  investigate  reconstruction  of  the  reduced  scattering  property  of  a  target  using 
continuous-wave  measurements.  Figure  2  lists  the  results  simulated  for  a  single  target  having  both  absorption  and 
scattering  contrasts  over  the  background.  The  homogeneous  background  (for  all  the  simulations  conducted  in  this  work) 

is  set  at  jUa  (absorption  coefficient)  of  0.002mm"1  and  jJi's  (reduced  scattering  coefficient)  of  1.0mm'1.  The  jda  of  the 

target  is  set  at  0.004mm"1  and  0.02mm"1;  the  jJi\  of  the  target  is  set  at  2.0mm"1  and  10mm"1.  The  target  is  a  sphere  of 

10mm  in  diameter  and  located  along  the  middle-sagittal  plane  (the  TRUS  imaging  plane)  at  a  depth  of  10mm.  Figure 
2(a)  shows  the  results  for  the  target  at  a  depth  of  10mm  and  placed  20  mm  longitudinally  from  the  mid-line  of  the 
sagittal  plane.  Figure  2(b)  shows  the  results  when  the  target  is  at  the  mid-line  of  the  sagittal  plane.  The  sagittal  and 
coronal  locations  of  the  target  are  correctly  reconstructed.  Cross-coupling  between  the  absorption  and  scattering 
reconstructions  out  of  the  continuous-wave  data  is  not  observed,  however,  the  optical  properties  are  underestimated. 
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(a)  Target  shift  20mm  longitudinally  from  the  mid-point 


(b)  Target  at  the  mid-point 


Figure  2  Reconstruction  of  the  absorption  and  reduced  scattering  of  a  target  located  in  the  sagittal  plane  based  on 
continuous-wave  data.  For  each  set  of  4-image,  the  top  row  corresponds  to  the  sagittal  view,  and  the  bottom  row 

corresponds  to  the  coronal  view. 
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(b)  Experimental  results  for  a  cylinder  target 
Figure  4.  Recovering  of  a  target  at  different  lateral  positions. 

3.2  Reconstruction  of  a  target  off  the  mid-sagittal  plane - Simulation  and  experimental  results 

The  geometry  of  the  trans-rectal  NIR  applicator  allows  3 -dimensional  imaging.  The  reconstruction  of  a  target  located  off 
the  middle-sagittal  plane  is  illustrated  in  Fig.  4.  Figure  4  (a)  shows  the  target  being  displaced  from  left-lateral  to  right- 


Proc.  ofSPIE  Vol.  7174  71741U-5 


lateral  at  5mm  per  step.  The  target  is  a  sphere  of  10mm  in  diameter  with  fla  of  0.0056mm"1  and  jil's  of  1.03mm'1.  The 

experiment  results  are  shown  as  Fig.  4(b)  for  a  target  of  cylinder  shape  (10mm  in  diameter  and  16mm  in  length)  having 
the  same  optical  properties  of  the  target  in  simulation  of  (a).  Both  the  simulation  and  the  experimental  results  indicated 
that  a  target  off  the  middle-sagittal  plane  may  be  recovered  by  NIR  only,  but  at  slightly  displaced  spatial  location. 
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(c)  Multiple  targets  with  same  absorption  but  different  scattering  coefficients 
Figure  1.  Reconstruction  of  multiple  targets 


3.3  Reconstruction  of  multiple  targets  on  the  mid-sagittal  plane - Simulation 

The  performance  of  differentiating  multiple  targets  on  the  mid-sagittal  plane  is  evaluated  by  simulation.  The  simulation 
is  performed  for  two  sphere  targets  of  10mm  in  diameter.  In  (a)  the  two  targets  have  the  same  fia  =  0.004mm"1  and 
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jUrs  =  2.0mm"1.  At  40mm  longitudinal  separation,  the  targets  can  be  separated  on  both  jUa  and  /u's  images.  At  closer 
separation,  overlapping  of  the  two  targets  is  seen  on  /ua  before  on  ju's .  In  (b),  the  two  targets  are  20mm  apart,  the  ju's 
of  both  is  2.0mm"1,  and  the  jUa  values  are  set  at  0.004mm'1  and  0.02mm"1  respectively.  In  the  jUa  image,  only  the 
higher  contrast  target  is  recovered.  While  in  the  jLl's  image,  the  lower  contrast  target  is  recovered  but  with  cross-coupling 
from  the  higher  contrast  target.  In  (c),  the  jda  is  0.004mm"1  for  both,  and  the  /u's  is  set  at  2.0mm"1  and  10mm"1  for  the 
two  targets.  The  two  targets  are  recovered  on  jLl's  image,  but  barely  separated  on  the  jUa  image.  The  above  simulation 

and  experimental  results  demonstrate  the  utility  of  reconstructing  both  /ua  and  ju's  based  on  continuous-wave 
measurement  even  though  the  recovery  of  the  target  is  not  satisfactory  for  all  conditions. 

4.  IN  VIVO  RESULTS 

The  results  of  in  vivo  trans-rectal  NIR  tomography  and  TRUS  of  the  canine  are  presented  in  Fig.  6  for  the  NIR 
absorption  images  and  Fig.  7  for  the  NIR  reduced  scattering  images.  All  NIR  images  are  displayed  with  consistent  color 

scales  for  both  jUa  and  jj!s  .  The  dimensions  of  US  and  NIR  images  are  60mm  (cranial-caudal) x3 0mm  (posterior- 

anterior)  for  sagittal,  40mm  (right  lateral  -left  lateral  )x30mm  (posterior-anterior)  for  axial,  and  60mm  (cranial- 
caudal)  x40mm  (right  lateral  -left  lateral)  for  coronal  views,  respectively. 

4.1  Hyper-absorptive  indication  of  the  TVT  on  NIR  image 

The  three  sagittal  TRUS  images  (Fig.  6)  were  taken  at  the  mid-line  of  the  prostate,  the  middle  portion  of  the  right  lobe  of 
the  prostate,  and  the  right  lateral  edge  of  the  prostate.  The  hypo-echoic  region  LI  indicates  an  abnormal  intra-prostatic 
mass.  The  large  hypo-echoic  region  L2  indicates  an  abnormal  mass  anterior  and  caudal  to  the  prostate  that  may  have 
connection  with  LI.  The  NT  on  the  right  lobe  US  image  denotes  the  needle  trajectory  for  introducing  the  TVT  cells. 
Longitudinal  hypo-echoic  regions  including  L3  are  seen  along  the  NT.  The  three  sagittal  NIR  images  correspond  to  the 
US  taken  at  the  middle  portion  of  the  right  lobe,  10mm  medial  to  it  and  10mm  lateral  to  it,  respectively.  Hyper- 
absorptive  regions  on  sagittal  NIR  images  correspond  to  LI,  L2,  and  L3  on  sagittal  US  images.  The  lOmm-medial  NIR 
image  shows  reduced  contrast  for  the  absorptive  masses,  and  the  lOmm-lateral  NIR  reveals  connected  strong  absorptive 
masses.  The  bladder  is  seen  as  slightly  hypo-echoic  at  the  2  more  medial  sagittal  US  images;  however,  no  hyper- 
absorptive  mass  is  presented  at  the  left-most  region  of  the  sagittal  NIR  images. 

The  three  axial  TRUS  images  were  taken  at  the  cranial  edge  of  the  prostate  crossing  LI,  the  caudal  edge  of  the 
prostate  crossing  L2,  and  the  peri-rectal  region  crossing  L3.  The  axial  US  images  show  a  small  hypo-echoic  intra¬ 
prostatic  mass  at  cranial  side  of  the  prostate,  the  bulging  of  the  right  lobe  and  extending  of  the  L2  over  the  prostate  mid¬ 
line  that  correlates  with  the  findings  on  mid-line  sagittal  US,  and  large  peri-rectal  hypo-echoic  mass  cranial  to  the 
perineum.  Of  the  five  axial  NIR  images,  the  longitudinal  positions  are  10mm  apart  from  cranial  to  caudal.  The  axial  NIR 
image  A2  should  contain  LI,  and  the  axial  NIR  image  A5  should  contain  L3. 

The  coronal  NIR  images  correspond  to  5mm,  10mm,  and  15mm  anterior  to  the  anterior  edge  of  the  rectal 
lumen.  The  hyper-absorptive  mass  indicative  of  LI  is  seen  medial  to  the  hyper-absorptive  masses  indicative  of  L2  and 
L3.  The  anatomies  of  the  hyper-absorptive  masses  on  the  NIR  images  in  3-views  agree  with  the  hypo-echoic  regions  on 
the  US  images  in  2-views. 

The  US  and  NIR  imaging  were  also  performed  at  the  middle-line  of  the  left  lobe  where  no  abnormal  features 
were  found  on  the  US,  and  globally  homogenous  and  low  absorption  was  seen  in  the  NIR  images. 

4.2  Hyper-scattering  indication  of  the  TVT  on  NIR  image 

The  hyper- absorptive  masses  in  Fig. 6  being  suspicious  of  TVTs  are  shown  in  Fig.  7  with  different  patterns  of  the 
reduced  scattering  contrast,  indicating  minimum  cross-coupling  between  the  reconstructions  of  jda  and  jj!s  from  the  in 

vivo  continuous-wave  measurement.  In  the  scattering  NIR  images,  L3  has  much  higher  contrast  than  do  the  other  masses 
corresponding  to  LI  and  L2.  The  contrast  elsewhere  is  also  slightly  more  homogenous  compared  with  that  in  the 
absorption  image. 
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Figure  2.  Trans-rectal  NIR  imaging  (absorption)  in  3-views  vs  TRUS  imaging  in  2-views 
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Figure  7.  Trans-rectal  NIR  imaging  (reduced  scattering)  in  3-views  vs  TRUS  imaging  in  2-views 
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4.3  Confirmation  of  the  TVT  growth 


The  mass  L3  indicated  by  the  US  and  NIR  correlates  with  a  large  peri-prostatic  nodule  along  the  needle  track  site 
determined  by  digital  rectal  examination  prior  to  imaging.  Necropsy  confirmed  that  L3  was  a  peri-rectal  nodule  along 
the  needle  track.  The  necropsy  also  confirmed  a  large  nodule,  corresponding  to  L2,  anterior  to  the  prostate  at  the  caudal 
side  of  the  prostate  that  extended  from  the  right-lateral  aspect  of  the  prostate  to  the  middle-line  of  the  prostate.  Exposing 
the  interior  of  the  prostate  also  confirmed  a  nodule  at  the  right  lobe  of  the  prostate  that  corresponds  to  the  mass  LI. 
Histological  examination  confirmed  all  nodules  to  be  TVT. 

5.  DISCUSSIONS  AND  CONCLUSIONS 

Among  the  hypo-echoic  masses  shown  on  the  US,  LI  was  intra-prostatic,  L2  was  largely  anterior  to  the  prostate,  and  L3 
was  caudal  to  the  prostate.  These  3  masses  are  however  shown  at  similar  depths  on  NIR  images.  This  is  attributed  to  the 
un-even  NIR  sensitivity  profile  along  the  posterior- anterior  direction  [12].  It  is  however  noted  that  the  existence  of  TVT 
nodules  have  been  recovered  by  use  of  only  the  NIR  information.  If  the  anatomic  location  of  the  TVT  available  from  the 
US  can  be  utilized  as  the  spatial  prior  for  NIR  image  reconstruction,  the  localization  of  the  TVT  nodules  on  NIR  images 
will  certainly  be  more  accurate. 

It  is  known  that  most  of  the  prostate  cancers  are  presented  as  multi-focal  [18,  19].  It  must  be  noted  that  the 
intraprostatic  TVT  tumors  were  initiated  and  developed  during  this  project  in  a  non-immunosuppressed  canine  model. 
The  subject  in  this  study  developed  multiple  TVT  nodules  intra-prostatically  and  peri-prostatically.  Although  the  sites  of 
the  TVT  tumors  were  not  totally  confined  to  the  prostate,  the  successful  imaging  of  multiple  TVT  tumor  nodules  by 
trans-rectal  NIR  tomography  implies  the  feasibility  of  recovering  multiple  intra-prostatic  TVT  tumors.  All  of  the  TVT 
tumor  nodules  developed  in  this  subject  were  shown  strongly  hypo-echoic  on  the  US,  which  correlates  with  the  strongly 
hyper-absorptive  and  moderately  hyper-scattering  findings  on  NIR  imaging.  The  bladder,  on  the  other  hand,  shows 
hypo-echoic  on  the  ultrasound;  however  it  is  not  shown  as  hyper-absorptive  or  hyper-scattering  on  the  NIR  images.  This 
confirms  that  trans-rectal  NIR  tomography  has  higher  specificity  of  imaging  the  regions  hypo-echoic  to  TRUS. 

It  has  been  well-demonstrated  that  cancers  in  organs  like  breast  are  presented  in  NIR  tomography  as  having 
increased  absorption  and  different  reduced  scattering  contrast  over  the  normal  tissues  when  imaged  intact,  in  situ  and  in 
vivo.  However,  there  was  previously  no  information  on  the  optical  contrast  of  prostate  cancer  compared  to  normal 
prostatic  or  peri-prostatic  tissue  taken  at  intact  prostate  in  situ  and  in  vivo.  Our  work  presents  a  new  paradigm  for 
prostate  imaging  that  shows  promises  of  improving  the  specificity  of  TRUS  imaging  by  augmenting  it  with  trans-rectal 
NIR  tomography. 

In  conclusion,  this  work  reported  in  vivo  imaging  of  TVT  tumors  in  the  canine  prostate  and  pelvic  canal  by 
trans-rectal  NIR  tomography  coupled  with  TRUS.  The  TVT  tumor  nodules  show  strong  NIR  absorption  and  moderate 
scattering  contrasts  over  the  normal  prostatic  and  peri-prostatic  tissues.  Correlation  of  the  TVT  locations  is  found 
between  trans-rectal  NIR  and  TRUS  results. 
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